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1.0 INTRODUCTION AND SUMMARY 


The AIRS study consists of two phases performed over a period of 
two years. The first phase of this study was devoted to define the 
functional requirements of AIRS, to define the testing required to 
demonstrate its performance, to provide an assessment of the technology 
and cost involves in such an undertaking, and to prepare a specification 
based on the above findings. This phase of the AIRS study were 
documented in Vol . II and the accompanying functional specification of 
the final report. 

The purpose of the Phase II study is to develop a computer model of 
the AIRS, as defined in the Phase I sutdy, to perform a computer 
simulation of the AIRS, and provide simulation results predicting the 
performance of the AIRS. The emerging TDAS requirements and techniques 
to enhance the AIRS capabilities were also incorporated into the AIRS 
baseline during this period. This volume of the final report serves to 
document the Phase II effort. The accompanying AIRS simulation software 
package has been installed at the Goddard CLASS computer. Instructions 
on using the simulation is contained in a CSS file called AIRS. 

1.1 Organization and Summary 

This volume is organized into 5 sections and 11 appendices. 

Section 2 discusses from a system operation point of view the autonomous 
and integrated aspects of the operation of the AIRS. The advantages of 
AIRS compared to the existing SSA receive chain equipment are 
highlighted. The three modes of AIRS operation are addressed in detail. 

The configurations of the AIRS are defined as a function of the 
operating modes and the user signal characteristics in Section 3. Each 
AIRS configuration selection is made up of three components: the 
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hardware, the software algorithms and the parameters to be used by these 
algorithms. These components and the rules for their selections are 
treated in this section. Notice that the first two components are 
modeled in the AIRS simulation software. The final parameter selection 
process is to be accomplished via the simulation software. 

A comparison between the AIRS and the wide dynamics demodulator 
(MOD) is provided in Section 4. They are compared in terms of their 
functional characteristics and implementation techniques. This serves 
to highlight some of the improvements that are designed into the AIRS. 

In addition, some unique features of AIRS are given. 

The final section of this report is devoted to describing the 
organization of the AIRS analytical/simulation software. Since the 
algorithms developed in the software simulation can be transferred to 
the AIRS microprocessor with minor modifications, a strict Monte Carlo 
technique is adopted whenever possible. In the cases where the Monte 
Carlo technique is impractical and/or overly time consuming, the 
software is supplemented with analysis. 

The modeling and analysis for simulating the performance of the PN 
subsystem is documented in detail in Appendix A and B. The acquisition 
performance is treated in Appendix A. The acquisition performance is 
determined based on the assumption that the code uncertainty region is 
uniformly distributed. In practice, a truncated 3-sigma Gaussian 
distribution is more appropriate since the actual received signal code 
epoch is more likely to be within the center of the predicted 
uncertainty region. If this is the case then the search should spend 
more time in the vicinity of the center of the uncertainty region. It 
is estimated that the acquisition time determined with the AIRS 
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simulation could be improved by a factor of 2 to 4. The analysis 
employed in Appendix A is not restricted to CCD matched filters. If a 
different implementation is used, e.g., using digital matched filter, 
the AIRS simulation is also appliable. However, when a digital matches 
filter implementation is used, the signal must first be sampled and 
hard-quantized. This introduces a 2 dB loss to the SNR. This factor 
must be accounted for when using the simulation. 

The tracking performance of the PN subsystem is treated ’,n appendix 


B. 

The frequency acquisition technique using a frequency-locked loo^ 
is documented in Appendix C. 

The portion of the AIRS simulation software concerning the PN and 
the frequency acquisition system is documented in Appendix D. A user's 
manual and sample runs are also given. 

The portion of the ARIS simulation regarding the phase-locked loop 
and the bit sync is documented in Appendix E. Sample runs are also 
given. 

The rest of the Appendices are a collection of miscellaneous 
items. Appendix F shows how autonomous Doppler compensation can be 
implemented. Appendix G addresses the technology aspects of employing 
CCDs for PN acquisiton. Appendix H shows how the AIRS AGC functions and 
how a signal strength indicator, accurate to G.2 dB, can be 
implemented. Appendix I shows that for all practical purposes, the 
Manchester symbol ambiguity can be resolved within 20U0 symbols. 

Finally, an expression is given in Appendix K to assess the allowabU* 
rms phase jitter in order to maintain a mean slip time of over 90 
minutes. 
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2.0 AUTONOMY AND INTEGRATED OPERATIONS FOR AIRS 


This section discusses from a system operational point of view the 
autonomous and integrated aspects of the Automated Integrated Receiver 
AIRS (AIRS). The advantages of AIRS compared to the existing SSA 
receive chain equipment are highlighted. Finally, the fhree different 
inodes of AIRS operations are discussed. 

2.1 Introduction 

The AIRS is designed to be a stand-alone receiver capable of 
autonomous, integrated operation for the S-band Single Access (SSA) 
return link of TDRSS. By autonomous operation, it is meant its 
capability to demodulate, detect and decode data, and to derive ranging 
and Doppler information from an IF input with a minimal dependence on an 
external executive computer such as the Automated Data Processor 
Equipment (ADPE) at White Sands Gr» >nd Terminal (WS r ' M . In order to 
achieve autonomous operation, the AIRS is capable of controlling its own 
receiver functions based on given a priori information as well as 
currently measured naracteristics of the input signal. The AIRS is 
also an integrated reco ; <er since the subsystems required for PN 
despreading, carrier and clock recovery, symbol detection and Viterbi 
decoding are designed, tailored, and interconnected to function -_s a 
single dedicated unit. 

Figure 2.1 depicts the differences between the autonomous 
integrated approach and the conventional centralized approach to a 
receive system design. Notice that in the conventional approach the 
subsystems are controlled by a central network controller which is 
external to the receive system. There are no communications and 
intersections among the subsystems except for che required signal 
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interfaces. In the autonomous integrated approach, the receive system 
accepts setup commands from the network controller initially. After 
initial acquisition, the receive system is normally detached from the 
controller except for responses to an interrupt command. This is 
indicated by the dashed line. Notice that the subsystems are integral 
parts of the receive system and are interconnected among themselves 
through the local controller. As we shall see, this autonomous 
integrated structure provides several advantages over the conventional 
approach. 

2.2 Advantages of the AIRS Concept 

Figure 2.2 highlights the differences between AIRS and the 
corresponding portion of the current SSA Receive Chain. In the current 
implementation, the units, i.e., Down £onverter/Ooppler Correlator (C/C) 
Low Rate Demodulation (LRD), Medium Rate Demodulator (MRD), Low and 
Medium Rate Bit Synchronizers (L/MRBS) ar.d the Viterbi Decoders (VD), 
are controlled directly by the ADPE through the Service Control Unit 
(SCU). The SCU merely acts as a common point of interface between the 
ADPE and the various units. There are two types of comnands from the 
ADPE: static setup comnands and dynamic Doppler commands. The one-time 

setup comnands are used to configure the parameters of the various 
units. The dynamic commands are provided continuously throughout a 
mission and they are the carrier Doppler estimates (used by the C/C) and 
the PN code chip rate Doppler estimates (used by the LRD). The 
individual unit reports its status to ADPE through the SCU. The 
following features of the current system are worth noting: 

t Doppler correction commands are required from the ADPE at all 
times for proper operation. 
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• There are no provisions for the ADPE to issue real-time chanyes 
to unit parameters to take advantage of the information p ovided 
by the status signals. The only ADPE reaction is to restart the 
system once an anomaly is detected. 

• The individual units function independently of one other. 

The AIKS also takes setup commands from and reports its status to 
ADPE as in the current implementation. However, the AIRS only uses the 
Doppler correction commands during acquisition. After acquisition, the 
AIRS generates its own Doppler estimates. Therefore, the .nterface 
between the ADPE and the AIRS can be removed after the link has been 
established except for configuration changes. During acquisition, the 
setup commands are similar to those issued by the ADPE. However, after 
acquisition, the AIRS can generate more accurate link characteristics by 
its own measurement than those provided by the a priori information. 
These measurements are fed to a decision algorithm (expert system) 
through a common data bus. The outcomes are used to instrut the local 
controller to modify the parameters of the functional module: to achieve 
the best desired performance. The key features of the AIRS are: 

• Detached Operation from ADPE 

• Real-Time Intelligent Operation 

• Integrated Operation 

• Self Diagnostics 

In what follows, we shall address the advantages provided by the AIRS 
implementation. 

2.2.1 Advantages of Autonomous Operation 

There are four major areas where the AIRS autonomous operation is 
advantageous. These are: 
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• Reduced Dependence on ADPE 
t Ease of Modification 

• Real Time Operation 

• Improved Performance 

2. 2. 1.1 Reduced ADPE Dependence 

Since the AIRS performs its own Doppler compensation once it has 
acquired, it can operate independent of the ADPE during the tracking 
operation. The ADPE can be freed to devote its attention to other 
processing needs of the ground station throughout the rest of the 
mission. Also, the receive chain operation will not be interrupted by 
possible disruptions in the ADPE. To accomodate emergency operations, 
the AIKS allows manual entry to simulate ADPE commands. 

2. 2. 1.2 Ease of Modification 

Currently, the setup commands from the ADPE contain both the data 
link characteristics (e.g., Data Group, Mode, expected signal power, 
data rates, etc.) as well as the settings for selectable unit parameters 
(e.g., search rate, IF bandwidths, loop bandwidths, etc.). These 
parameters are determined by the ADPE based on the estimated return link 
characteristics. A similar set of software for making that decision are 
incorporated into the AIRS control so that the only setup commands 
required frcm ADPE are those pertaining to data link characterization. 
This way, the AIRS configuration and parameters setup are completely 
ir -pendent of the ADPE. Any subsequent changes in the receiver 
structure as a result of new mission requirements or equipment 
modifications can be handled locally by the AIRS controller and will not 
impact t.ie ADPE. In other words, future modifications on the AIRS will 
not alter the AIRS/ADPE protocol and interface. Presently, major ADPE 
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software revision will be required since the changes must be made 
compatible with other activities of the ADPE. 

2. 2. 1.3 Real Time Operation 

The ADPE setup commands are geared to achieving the “best" 
performance based on the a priori information on the user link 
characteristics scheduled for mission support. The receive system 
parameters are then frozen throughout the mission. However, after the 
receive system has acquired the user signal, the receive system itself 
is capable of generating better estimates on the user signal 
characteristics than those based on a priori information. The best set 
of receive system parameters may not necessarily be the same as one 
dictated by the ADPE. In what follows, we shall dwell on some examples 
to illustrate this point. 

2. 2. 1.3.1 Signal Strength 

The coherent AGC in the receive system is capable of providing an 
accurate measurement on the received signal strength. This information 
can be used to optimize the carrier and bit sync tracking bandwidths. 

For example, if the signal strength is higher than the a priori 
prediction, then the loop bandwidths can be opened up to allow for 
tracking errors due to phase noise and dynamics. Similarly, the 
bandwidths can be narrowed if the converse is true. 

2. Z. 1.3. 2 Wide Dynamics Operation 

After acquisition, the receive system can measure the signal 
Doppler accurately. It does not have to rely on the ADPE for Doppler 
correction. This is particularly advantageous during a powered 
flight. Since the ADPE can only predict the burn occurrence to within 
+9 secs, the Doppler correction provided to the SSAR chain can be off by 
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a frequency rate up to 765 Hz/sec for +9 secs. The tracking loops in 
the current receive system must be able to track through this 
transient. Since the Doppler is internally compensated in the AIRS, the 
tracking loops will only see the frequency rate for a period of time 
proportional to the Doppler update rate which is on the order of a 
second or less. Since a tracking loop is more likely to slip and lose 
lock while operating under stress, the probability of loss of lock is 
proportional to the duration of the stress application. Therefore, the 
AIRS implementation is more tolerant to signal dynamics. 

2. 2. 1.3. 3 Reacquisition 

Occasionally, the communication link breaks down due to an anamoly 
during a support mission. Up until the signal is lost, the tracking 
receiver has the most current estimate on the carrier frequency, PN chip 
rate, and PN epoch uncertainty. These uncertainties are much smaller 
than the initial uncertainty predictions given by the ADPE. The AIRS 
can therefore use this information to help reduce the uncertainty range 
to be searched during reacquisition. This way, the reacquisition time 
can be reduced. In the current system (except to some degree in the 
WDD), this information is ignored and there is no distinction between 
acquisition and reacquisition. 

2. 2. 1.4 Performance Improvement 

The ability to monitor the link conditions allows the AIRS to 
better match the receiver parameters with the i .coming signal 
characteristics. This translates directly i.'to BER, tracking and 
ranging performance improvements. The AIR5, is al ^o capable of resolving 
certain cases of I,Q channel ambiguities ^.g. wien the data rates are 
close to each other) that the current s.vs em is not capable of doing. 
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The AIRS relies on its monitor signals to achieve this goal. However, 
these signals are either unavailable, inaccessible or impractical to be 
made available to the ADPE in the current system. 

Through the use of optimal receiver algorithms, the threshold 
performance of the AIRS is greatly enhanced. Operationally, this means 
that the AIRS can be used to support degraded operation involving 
malfunctioning users. 

2.2.2 Advantages of Integrated Operation 

The AIRS subsystems are designed and interconnected to achieve the 
best receive system function as a single entity. This provides several 
advantages over the current system. 

2.2.2. 1 Reduced Hardware and Interface Complexity 

Since the subsystems are not built as individual units, hardware 
that perform similar functions can be shared. For example, the analog- 
to-digital conversion function can be shared by both the carrier and the 
bit sync loops. Other obvious examples are power supplies and interface 
signal processors. In addition, since a single vendor is responsible 
for the AIRS, the interface requirements between the units (e.g. the LRD 
has to interface with the SCU and the L/MRBS) can be eliminated. 

2. 2. 2. 2 Ease of Maintenance and Seif Diagnostic Capability 

Reduced complexity translates directly into easier maintenance. In 
particular, because of the digital implementation approach chosen by the 
AIRS, maintenance is further simplified by reduced analog component 
counts (required only at the front end of the receive sysem before A/D 
conversion) and extensive use of digital components/numerical 
algorithms. Digital components eliminate the traditional analog 
problems associated with drifts, aging and temperature sensitivity which 
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requires periodic adjustments. In addition, the numerical algorithms 
are amendable to self diagnostic procedures which greatly simplify 
trouble-shooting. Because of the built-in intelligence of the AIRS, 
self diagnostics ca r be achieved in many levels. During operation, any 
anamoly that causes a signal dropout will initiate a self-diagnostic 
procedure to identify the origin of the failure mechanism, i.e., whether 
it was a result of the receiver failure or an anamoly in the link 
condition. The diagnostic procedure, whether during power-up or during 
operation, will be used to isolate fault to the board level. 

2. 2. 2. 3 Performance Improvements 

Integrated operations mean that the AIRS subsystem functions are 
coupled together. This is most noticeable for the various tracking 
functions (i.e., code, carrier and bit sync loops!. The AIRS employs a 
rather sophisticated data-aided loop (DAL) design for carrier tracking 
and bit sync. The advanages of a DAL implementation has been 
demonstrated to some extent by the WDD (for which a rather simple DAL 
implementation is employed). In particular, the inherent bit sync/data 
detector performance of the WDD has been demonstrated to be superior 
than the performance of the WDD used in conjunction with the L/MR bit 
sync. No doubt, the more sophisticated nature of the AIRS will further 
enhance this performance capability. 

The DAL is also a better design in terms of the carrier tracking 
and threshold perfomance. The AIRS will have a lower threshold and 
smaller rms phase jitter for a given CNR input. The AIRS DAL also 
eliminates undesirable lock points on the S-curve. Thus, it is capable 
of differentiating between the I and Q channel for unbalanced QPSK 
signals. This eliminates certain I.Q channel ambiguities associated 
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with the current system. In addition, since the bit sync and carrier 
loop acquire simultaneously instead of sequentially, acquisition time is 
'ced 


The AIRS employs a Doppler compensation loop to minimize the loop 
stresses. The Doppler loop is coupled with the DAL and improves the 
tracking and threshold performance. Furthermore, the PN loop is aided 
by the DAL to further minimize the PN loop stress. This improves the 
ranging perfoinance. 

2.3 AIRS Operation Modes 

The AIRS has th^ee modes of operation. The two autonomous modes 
(Modes 1 and 2) are very similar except for the receiver 
configurations. The other mode, the test mode, is a mode where the 
receiver is operated manually by a keyboard/display via an external 
interface. The receiver parameter selections are menu-driven to allow 
friendly user interface. All three modes can be controlled by a local 
controller. Mode 1 and Mode 2 are normally controlled by a remote 
controller/computer such as the ADPE. 

2.3.1 Mode 1 - Normal Mode 

The flow diagram for the autonomous modes (Modes 1 and 2) of 
operation is shown in Figure 2.3. The AIRS accepts the setup commands 
from the ADPE. The setup commands define the user signal 
characteristics as shown in Figure 2.4. Based on these commands, the 
receiver will be configured and parameters will be selected by the AIRS 
control unit to best-match the data characteristics. The receiver 
configuration and parameters are distinct for Mode 1 and Mode 2. For 
Mode 1, the flow diagram for acquisition is shown in Figure 2.5. The 
various loops are configured to take advantage of all the information of 
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Mode Select (1, 2, or 3) 

PN Code Tap Setting 
I and Q Data Symbol Rate 
I/Q Power Ratio 
Code Search Direction 
Epoch Delay 
Search Range 

Register A PN Code Kernal 
Register B PN Code Kernal 
DG-l/DG-2 Select 
BPSK/QPSK Mode 
Acquisition Start 
Code/Uncoded 
Combined/Uncoded 
Differential Format 
EIRP 

Figure 2.4. List of Setup Command Functions 
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the data characteristics of the incoming signal. For example, both the 
I and Q channel signals are processed for acquisition and tracking «.id 
the respective data rates are used to set the * ious bandpass filter 
bandwidths. (Currently, only the strong channel is used.) During 
acquisition, the ADPE also provides AIRS with the dyanmic commands which 
are the carrier Doppler and code Doppler corrections. Unce AIRS has 
acquired, it signals ADPE to discontinue supporting the Doppler 
compensation. From this point, the AIRS i c . comnletely independent of 
the ADPE and operates functionally as shown in Figure 2.6. However, the 
ADPE c n interrupt AIRS at any time to regain control the AIRS. 

During autonomous oepration, the AIRS continually monitors for any 
system anamoly such as loss of lock. If such an anamoly occurs, the AIS 
enters a reacquisition mode. In this mode, the AIRS tries to reacquire 
the signal in the most efficient manner by using all the signal 
information mointored by the AIRS prior to the anamoly. If the 
reacquisition is successful, the AIRS resumes autonomous oepration. If 
the AIRS ~.'Hnot acquire, after a preassigned amount of time determined 
by the search range uncertainties, the anamoly will be reported to ADPE 
for attention. The . i.TS then continues to reacquire .ntil the ADPE 
issues new setup command. 

2.3.2 Mode 2 - Felxible Data Format Mode 

The operational flow of this mode is very similar to iode 1 except 
for the receiver configurations during acquisition and tracking. The 
purpose of this mode if to allow the user to switch their baseband data 
charateristics (data rates, oata formats) without having to reacquire 
the PN code and the carrier. To achieve this goal, the carrier loop and 
the bit sync must be decoupled. In addition, the AIRS cannot use the 
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data rates information in the PN and the carrier loops. It can only 
assume that the highest data rate used is compatible with the user 
transmitter iIRP. Figure 2.7 shows the acquisition flow diagram for 
mode 2. Note that two separate bit syncs are used and they acquire 
independently . After acquisition the functional diagr. is shown in 
Figure 8. Note that in this case, the bit syncs can lose lock 
temporarily while the data rates are being switched. However, the PN 
and carrier loops will still be in lock. Therefore, no PN and carrier 
loop reacquisition are necessary. 

2.3.3 Mode 3 - Test Mode (Manual Mode) 

In the test mode, the receiver is configured via external 
interactive keyboard commands. The receiver parameters must also be 
selected externally. In this mode, the AIRS allows external access to 
various portions of the receiver. This also allows an operator to 
experimentally optimize the receiver performance by deviating from the 
autonomous control program of the receiver. This feature is helpful if 
the AIRS is to be further modified to accomodate new requirements after 
its development. 

3.0 AIRS CONFIGURATIONS 

Since the AIRS must accomodate a variety of data modulation schemes 
under three different modes of operation, it possesses a multitude of 
configurations. Surprisingly, the amount of hardware that needs to be 
reconfigured is relatively small. This can be attributed to the use of 
digital/microprocessor control techniques in the AIRS architecture. As 
will be demonstrated shortly, any change in data modulation can usually 
be accomodated via changing some nominal clock rates and parameters used 
in the software. 
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Figure 3.1 is the functional block diagram of the AIRS. The 
received IF siynal is split into two channels. These channels are first 
processed by the PN subsystems if they are spread. The I and Q channels 
are then coherently demodulated and A/0 converted to digital samples. 
These samples, at a relatively high rate are processed in the digital 
signal processing (OSP) subsystem to generate the detected bits or 
symbols, carrier frequency and phase errors, bit sync timing errors, as 
well as other monitoring signa’s. The data rates of the errors and 
monitoring signals at the output of the DSP subsystem are typically 
reduced to a few kHz so that they can be processed by the succeeding 
microprocessor software. Based on these errors and monitoring signals, 
the software residing in the microprocessor then makes fine adjustment 
on the control and clock inputs to the A/D and DSP subsystem so as to 
maintain coherent demodulation and synchronous data detection. The 
exact control applied to these systems is a function of the modulation 
format and the operating mode selected. 

Depending on the user data characteristics and the operating mode, 
the AIRS selects a particular configuration by arranging the hardware 
paths, the appropriate PSP blocks, the set of appropriate software 
algorithms, and a set of appropriate receiver parameters. 

The user data characteristics are determined from the interface 
data supplied by the service control unit (SCU). Figure 3.2 shows the 
data characteristics that can be derived from the original SCU data. 

The effective C/Nq can be estimated based on the minimum required E b /N 0 
and the data rate. The selection of the receiver mode ( IRM) can be made 
via external switches or as a new item supplied by the SCU data. 

The AIRS configuration table is shown in Figure 3.3. Depending on 
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DATA CHARACTERISTICS 

PH USE: 

iPM/ O ""= I ?:Q (DGl-1,2) 

1 = I ONLY (DG1-3) 

2 = NONE (DG2) 

SYMBOL FO RMAT: (=BIT FORMAT IF UNCODED) 

ISF/v B I -PHASE 

1 = NRZ 

2 = UNKNOWN 

MODULATION FORMAT: 

IMF/O = BF'SK (DG2) 

1 = 5QPSK (DG2) 

2 = UOF'SK (DG2) 

3 = BPSK*2 (DG1) 

4 = UNKNOWN (DG2) 

COD E RA TE: 

ICR/O = UNCODED 

1 = RATE 1/2 

2 = RATE 1/3 

3 = UNKNOWN 

RECEIVE R MODE SELECTION: 

IRM/U = NORMAL MODE 

1 = FLEXIBLE DATA FORMAT MODE 

2 = TEST MODE 

OT HER PARMETERS: 

I/Q POWER RATIO 
EFFECTIVE C/No 
I BIT RATE 
Q E<IT RATE 

DE RIVED INFO RMATION: 

1. SYMBOL RATE = BIT RATE / CODE RATE 

2. I CHANNEL EFFECTIVE C/No = EFFECTIVE C/N 0 * IQF'R/(1 + 

3. O CHANNEL EFFECTIVE C/No = EFFECTIVE C/N 0 / (1 + IQF’R) 

I OPR = I CHANNEL POWER / 0 CHANNEL POWER 


I QF'R) 


Figure 3.2. User Data Characteristics. 
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AIRS CONFIGURATIONS 


F'N CONFIGURATION 

2-CHANNEL PN TRACK. STRONG CHANNEL ACQ 
I -CHANNEL F'N TRACK AND ACQ 
E 'ASS F'N SYSTEM 

BIT SYNC CNOF I GURAT I ON 
B I -PHASE (NRZ 0 2*SYMB0L RATE) 

NRZ 

N/A 

FLL CONFIGURATION 
B I -PHASE (NRZ @ 2*SYMB0L RATE) 

NRZ 

B I -PHASE 

FLL CONFIGURATION 

IzQ-Qz I 

IzQ-OzI 

IzQ-OzI 

2 * (IzQ-QzI) 

IzQ-Qz I 

IF FLL CONFIGURATION 

DATA-AIDED 2-F'HASE COSTAS 

DATA-AIDED 4-F'HASE COSTAS 

DATA-AIDED 4-PHASE COSTAS 

2*DATA-AIDED 2-F'HASE COSTAS 

N/A 

N/A 

N/A 

N/A 

2*2— PHASE COSTAS 
4-PHASE COSTAS 


Figure 3.3. AIRS Configuration Table. 
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the data characteristics represented in Figure 3.3 by the various values 
of the identifiers IPN, ISF, IMP, and IRM (see Figure 3.2 for what they 
represent), the AIRS uses different configurations for the PN subsystem, 
bit syncs, the frequency-locked loop (FLL) and the phase-locked loop 


(PLL). 


3.1 Configuration Rules 

Before describing in detail each of the subsystems and laying out 
all the different configurations, it is instructive to understand the 
basic rules being applied. 

3.1.1 Normal Mode 

In the normal mode of operation, all the necessary information on 
the data modulation format are known. The AIRS will take advantage of 
this knowledge to configure the receiver to exploit the modulation 
structure. The tracking loops, i.e., PN tracking loop, carrier recovery 
loop and bit sync will operate on both the I and Q channel in a coupled 
fashion. For example, in the DG-1 mode, there are two independent data 
channels. However, since they both have the same carrier, the AIRS 
. carrier recovery loop processes the phase error samples from both 
channels. As another example, data-aided loops are used to obtain the 
optimum performance. 

3.1.2 Flexible Data Format Mode 

The flexible data format mode only affects the receiver 
configuration following PN despreading. Since the data modulation 
format cannot be assumed known at any given time, the carrier loops must 
be decoupled from the bit sync. 

3.2 AIRS Subsystem Hardware 

As shown in Figure 3.1, there are three distinct hardware 
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subsystems. They are the PN subsystem, the A/0 subsystem and the DSP 
subsystem. 

3.2.1 PN Subsystem 


The PN subsystem is used for DG-1 signal when the input signal is 
spread. Initial acquisition is accomplished by using the CCD matched 
filters when the input data rate is less than 50 Ksps.* The CCD 
technique speeds up the acquisition process considerably at low data 
rates by providing a parallel processing capability. At higher data 
rate, because of the low number of PN chips to data bit ratios, the 
advantages of the CCD matched filter diminishes (see Appendix A). In 
addition, less than 1 sec acquisition time can be accomplished using 
con/entional sequential detection schemes. Therefore the AIRS uses 
sequential detection for data rate higher than 50 Ksps. The use of 
sequential detection, however, does not involve additional hardware, 
since the correlators required for sequential detection are shared with 
the correlators required for tracking. 

AIRS uses only one PN channel for acquisition, since there is no 
simple way to take advantage of the signal power provided by the other 
PN channel to aid acquisition, with the same level of complexity in 
hardware. In case of DG1 mode 1 or 2, the PN signal with the higher 
power is used. For DG-2, the PN subsystem is bypassed. 

A dithered early-late gated loop structure is used for PN 
tracking. The dithered method minimizes the effects of hardware 
imbalances while providing a potential 3 dB improvement over the strict 
early-late gated loop. For DG1 mode 1 and 2, the dithering is done on 


^When Manchester coding is used, 25 Ksps. 
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both the I and Q channel so that error signals from both channels can be 
used. This provides a potential 3 dB improvement over a scheme which 
tracks only one channel if both the I-Q channels have the same power. 

For DG1 mode 3 where only the I channel is spread, the dithered loop 
assumes the conventional structure. The details of the PN subsystem 
operation are described in Appendices A and B. 

The CCD acquisition technique is inherently digital in nature so 
that the CCD acquisition unit interfaces with the AIRS microprocessor 
software directly. In order to tap the power of the AIRS software, the 
PN tracking loops and the portion of the PN acquisition system using 
sequential detection are implemented in a hybrid manner. The outputs of 
the correlators used in both systems are A/D converted. Loop filtering 
and sequential detection functions are then performed in software. 

Figure 3.4 shows the hardware arrangement for the tracking/sequential 
acquisition subsystem. The accumulator (ACM) is controlled by the 
dithering signal during tracking. For sequential acquisition, the ACM 
is controlled by the processor software. The digital data output to the 
AIRS processor will be on the order of 1 KHz or less. 

3.2.2 A/D Subsystem 

Tue A/D subsystem is shown in Figure 3.5. The synthesizer provides 
the local reference for coherent demodulation. The incoming data is 
filtered by the bandpass filter. The filter bandwidth is selected as a 
function of the incoming symbol rate. The bandwidth is selected so as 
to minimize the signal distortion on the filter output while maintaining 
the output signai-to-noise ratio to be as high as possible. The 
baseband signal at the mixer output is amplified to a suitable amplitude 
for A/D conversion by the gain control amplifier (GCA). The GCA serves 
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to close the coherent A6C loop. 

In the scheme depicted in Figure 3.5, two separate A/D converters 
are used to sample the inphase and quadrature channel. Another 
possibility is to use one A/D converter to sample directly at 4x the IF 
frequency of the incoming signal as done in the Wide Dynamics 
Demodulator (WDD). This eliminates the need for the two mixers. 

However, this technique is judged to be unsuitable for the AIRS 
application as discussed in Appendix I. 

The output of the A/D is routed to two accumulators (ACM). ACMj 
(and ACM-} in the lower branch) is the data accumulator because its dump 
clock T} is normally coherent with the symbol clock so that its output 
is an approximation to the optimum integrate and dump matched filter 
detector. Notice that clock for T 1 is higher than T $ . The midbit 
accumulator ACM 2 (ACM^ in tht lower branch) is normally dumped at the 
same rate at T^ however, its phase differs from T 2 by half a cycle. 
Tnerefore the accumulator output is approximately the integral of the 
first half of a symbol and the last half of a previous symbol. This is 
illustrated in Figure 3.6. However, the clocks for T} and T 2 can be 
related in a different way for some configurations. The A/D samples are 
also available directly as A} (B^. 

In the normal mode of operation, the A/D sampling clock, at an 
integer multiple of the data rate, is derived from the bit sync clock. 

In the flexible data format mode operation the sampling rate is set to a 
rate determined solely by the highest data rate that the signal power 
can support. This is done becasue the actual data rate can vary without 
the knowledge of AIRS. Tying up the sampling clock with the data rate 
does not provide any advantage. 
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(a) Data 




(b) A/D Sample Clock 



(c) T 1 Dump Clock 





(d) Dump Clock 


/N 


Fig. 3. 6. Typical Timing Relationships. (The dump clock 

transfers the ACM content and resets the ACM for the 
next summation. ) 
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3.2.3 DSP Subsystem 

The digital samples from the A/D subsystem must be further 
processed before they can be used as inputs to the microprocessor 
software. The DSP subsystem consists of several blocks, each of which 
performs a distinct function as shown in Figure 3.7. Depending on the 
receiver configuration, appropriate outputs of the DSP subsystem are 
routed to the software. Notice that in the flexible data format mode, 
there is a separate DSP block for bic sync and for carrier recovery. In 
the normal mode, the DAL DSP block is used. 

3.2.3. 1 Frequency Error Signal Generator 

The frequency error signal generator serves as the frequency error 
detector for the frequency-locked loop (FLL) and its implementation is 
shown in Figure 3.8. The dump rates for Tj and T 3 is equal to M x of 
the symbol rate where M is an integer between 5 and 15. The sampling 
rate is nominally set at 32 times the symbol rate until it reaches the 
sampling speed limit. At that point, it reduces to the largest integer 
multiple of the symbol rate that does not exceed the sampling speed of 
the A/D converter. 

The accumulator serves to filter as well as to slow down the data 
rate of the FLL error signal to match that of the software. Since the 
FLL bandwidth will be less than 10 Hz, the ACM can be dumped at a rate 
less than 1 KHz. 

Also shown is the signal to be processed by the software for lock 
indi.'tion. The lock indicator signal is proportional to cos awTj while 
the frequency error is proportional to sin aoJj when the FLL has 
acquired lock. Since AaJj is small compared to 1 radian by design 
during tracking, a comparison between the two signals can be used to 
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Fig. 3.8. FLL DSP. 
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indicate lock. 

3. 2. 3. 2 Data-Aided 2-Phase Costas Loop DSP 

Fiyure 3.9 shows the DSP block required for generating the DAL 
errors. Notice that T^, T 2 and T3 are all derived from the bit sync 
clock. Tj and T2 are the same and T-j is offset from them by half a 
symbol. The errors are output at a rate of roughly 10 x the loop 
bandwidths. The lock indicators consists of filtered versions of Aj, 
A 2 , and A3. When the system is locked, A 3 «» 0, A 2 « (cos <t>| , and 
A3 ® (sin (j»|. Bit sync lock can be indicated by comparing A^ and A 2 
while carrier lock can be indicated by comparing A 2 and A3. 

3. 2. 3. 3 Data-Aided 4-Phase Costas Loop DSP 

Data-aided 2-phase DSP from the I and Q channel are combined to 
form the 4-phase DSP as shown in Figure 3.10. For staggered QPSK. the I 
and Q sampling clocks are the same and the dump clocks from the I 
channel is offset from those of the Q channel by half a symbol. For 
general unbalanced QPSK (UQPSK) operation, the sampling clocks and dump 
clocks for the I and Q channels are unrelated. 


3. 2. 3. 4 Flexible Data Format Mode PLL DSP 

A novel approach is used for the flexible data format mode DSP for 
generating the phase error as shown in Figure 3.11. The Aq and Bq are 
dumped at 4 x the highest symbol rate allowed by the flexible data 
format mode. The (Aq,Bq) pair defines the phase angle of the signal 
vector. The KOM converts this phase angle to a phase error. The rate 
is slowed down by t, ' output accumulator for the AIRS software. 

Figure 3.12 shows the required programming for the ROM. The phase 
angle e, defined by x = Aq and y = Bq defines an address (x,y) of the 
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• DATA-AIDED 2-PHASE COSTAS LOOP DSP FROM I AND Q CHANNEL 



Q 



Q TIMING ERROR 
PHASE ERROR 


• SQPSK 

• ACM DUMP RATES 

Tj = T^ 3 (Symbol Rate)'* 

T 2 Offset by Half Symbol 

I&Q Dump Clocks Staggered by Half Symbol 


• UQPSK 

• ACM DUMP RATES 

Tj = T 3 = (I or Q Symbol Rate) - * 

T 2 Offset by Half Symbol 
I&Q Dump Clocks Unrelated 

Fig. 3.10. Data-Aided 4-Phase Costas Loop DSP. 







• ACM DUMP RATE 

• T = = (4 x Highest Symbol Rate)" 1 

Fig. 3.11. Flexible Data Format Mode PLL DSP. 
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ROM. The content of the ROM at that address is proportional to the 
entry in the table shown in the figure. For example, for BPSK 
modulation, the ROM content will be porportional to &-n if e is in 
quadrant 8. In general, if the modulation is not specified, then BPSK 
must be assumed for OGI and UPQSK must be assumed for UQPSK. It is 
interesting to note that a for UQPSK is selected as a compromise between 
a 1:1 and a 4:1 power split. When a = 45° this is optimal for balanced 
QPSK (1:1). When a= tan -1 (l/2) = 26.6°, this is optimal for (4:1) 
unbalanced QPSK. The penalty in not knowing the power split is an 
increase in jitter due to this 9° average offset. 

The signal for lock indicator derived from Aj is proportional to 
| cos 4 >| where <j> is the phase error. The signal derived from Ag is 
proportional to |sin <t>| . Lock indication can be obtained by comparing 
the two signals. 

3. 2. 3. 5 Flexible Data Format Mode Bit Sync DSP 

In this mode, the bit sync must operate independent of the PLL as 
reflected in Figure 3.13. Notice that the samples A^ come directly from 
the A/D converter. The integrate-and-dump is accomplished in this DSP 
blocks. The dump clocks are controlled by the bit sync synthesizer as 
done in a conventional DTTL. The detected data out of the full symbol 
ACM is also sent to the data detector DSP. 

The full sybmol magnitude is much larger than the mid symbol 
magnitude when the bit sync has acquired. They can be used for lock 
indication. 

3. 2. 3. 6 AGC DSP 

The processing block for AGC is shown in Figure 3.14. Every time 
the A/D sample registers an extreme value (e.g. 1 or 64 for a 6-bit A/D 
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converter), the counter outputs a one. This is accumulated for a number 
of samples so that when the ACM is dumped, its content is proportional 
to the out-of-range p r obability. The output rate to the software is 
approximately 10 x the AGC bandwidth. 

The mean of jA^ is proportional to the signal amplitude and is 
used in the software that computes signal strength. The technique is 
explored in Appendix H. 

3. 2. 3. 7 Data Detector DSP 

The purpose of the data detector DSP shown in Figure 3.15 is to 
provide soft (6 bit) and hard decision (1 or 0) detected symbols. In 
case of Manchester-coded symbols, the symbol is converted first to NRZ 
symbols. This conversion involves an ambiguity. The DSP block is 
capable of resolving this ambiguity. The technique used and analysis of 
its performance is documented in Appendix J. 

3.3 Software Algorithms 

There are four main functions of the microprocessor software. They 

are: 

• Loop Filtering 

• Synthesizer Control 

• Monitor 

• Logic 

The first two functions are required to close the carrier loop and the 
bit sync. The third function is implemented to monitor the status of 
the receiver. The last function is used to configure the receiver, 
select the appropriate parameters and control the sequence of various 
operational procedures such as acquisition or reacquisition. 

3 . 3 . \ Loop Filtering 
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The loop filtering algorithms takes an input digital 
sequence {x k } from the DSP subsystem and generates an output sequence by 
computing 




ax k a by k _j ♦ cy k _ 2 


The coefficients a, b, and c define the loop order (up to third order 
loop' _• the loop noise bandwidth. The selection of these coefficients 
is done during the setup phase based on the given user 
characteristics. They are modified if a better user characterization 
can be obtained, e.g., after the receiver has acquired so that there is 
a better estimate on the signal level. These coefficient are also 
changed when, for example, the PLL switches to a third-order loop from a 
second-order loop after the initial acquisition. When the DSP loop 
error signal samples are available from both I&Q channels, they are 
weighted by the respective power split and combined before they are 
processed by the filtering algorithm. This is done for QPSK (DG-2) 
signals. 

3.3.2 Synthesizer Control 

The filtered output y k is used to control the frequency (phase) of 
the appropriate synthesizers to close the carrier recovery and the bit 
sync loop. The filtered phase error controls the frequency of the local 
reference used to coherently demodulate the signal in the A/D 
subsystem. The filtered bit sync error is used to control the frequency 
of the bit sync synthesizer from which the T $ , T^, T^, and T^ clocks in 
the A/D subsystem and other clocks in the DSP subsystem are generated. 

3.3.3 Monitoring 
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Lock detector indicators and the signal strength indicator are 

lowpass filtered using a recursive difference equation similar to the j 

loop filtering operation described above. The appropriate filtered 

signals are compared to verify lock (see Appendix K). The signal 

strength indicator samples are averaged and scaled properly to provide | 

l 

the signal strength estimate. These operations are performed 

continusouly as long as the AIRS is in operation. | 

3.3.4 Logic 

I 

The logic portion of the AIRS software is probably the most • 

intricate of the software functions. It defines the acquisition, j 

reacquisition, self-test and diagnostics sequence. It sets the receiver 
configuration (hardware, software and parameter selections) as a j 

function of the set up commands. Currently, it is conceived to use a 
tree structure for the decision making. However, the AIRS concept is 
well adapted to be implemented as an expert system using artificial 
intelligence (AI) techniques. 

3.4 Parameter Selection 

The selection of the AIRS subsystem parameters such as carrier and 
bit sync loop bandwidths is Dased upon a tradeoff between a set of the 
drivers such as: 

• Acquisition time 

• Reacquisition time j 

• Mean slip time j 

• Phase noise j 

0 Dynamics 

0 BER 

The effects of these drivers have been considered in this report (see 


i 
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Vol . 1, Appendices F,L). Fortunately, the AIRS has been designed so 
that most of these drivers can be accomodated in a very simple manner. 

Typical carrier loop mechanizations require a compromise between 
the loop bandwidth for acquisition and tracking since the acquisition 
performance improves with a larger loop bandwidth while the tracking 
performance degrades. The AIRS gets around this problem by using a FLL 
to aid acquisition and a third-order loop to minimize the effect of the 
dynamics so that the PLL loop bandwidth can be designed based on mean 
slip time and phase noise alone. 

In what follows, we consider various parameters and suggest values 
to be used for the key subsystems of the AIRS. If the operation of AIRS 
is compatible with the WDD, the system parameters will be selected on 
that basis. However, the AIRS simulation will be used as a tool to 
optimize these selections when the AIRS is ready for a final design. 
3.4.1 Frequency-Locked Loop 

The PLL is first-order since a first-order loop has the best 
acquisition performance for this application. In order to acquire the 
±4.2 kHz frequency uncertainty within a short time (~ 1 sec) at the 
lowest acquisition signal level, a multiple bandwidth technique can be 
used. As an example, the FLL can initialize with a 6 Hz bandwidth, 
reduces to 3 Hz at the end of 0.2 sec and further reduces to .3 Hz at 
the end of 0.6 sec. 

At higher signal levels, the loop bandwidth can be switched 
sooner. The actual switching time to achieve optimum performance can be 
determined by the AIRS simulation. 

When the signal level is high (i.e., for symbol rate > 20 ksps), 
the use of FLL does not provide any advantage since a PLL can be 



- 49 - 



* 


'<jCinCt 


om 


designed to acquire rapidly and without any possibility of false-locking 
to data side-bands. The situation is considered in the PLL discussion 
to fol low. 

3.4.2 PLL 


The PLL is designed to first acquire phase with a second-order loop 
configuration and then to switch to a third-order loop for tracking. A 
selection of the bandwidth and loop order are shown in Table 3-1. The 
symbol rate is the sum of the I and Q channel. Notice that if the 
symbol rate is higher than 20 Ksps, a 5 kHz bandwidth loop is used to 
acquire frequency directly, without first aiding with a FLL. Again, the 
AIKS simulation can be used to optimize these parameters. 

The selection of the lower limit of the PLL is based on satisfying 
the worst case tracking signal to noise level of C/Nq = 28 dB-Hz. This 
will give a rms phase jitter of about 11° theoretically. This will also 
yield a slip time of about 90 minutes (see Appendix L). However, 
because of phase noise effects, the phase jitter will be slightly higher 
and the mean slip time will probably be on the order of 1 minute. 
Ideally, the phase jitter can be reduced by decreasing the loop 
bandwidth. In practice, oscillator phase noise of the system will begin 
to dominate at low loop bandwidths and introduce an additional component 
of the jitter which can make the total higher than the 11°. 

In Table 3-1, notice that the loop bandwidth selection is 
determined by the symbol rate. This, of course, assumes that the signal 
level is barely meeting the BER requirement. In many cases, the 
transmitter will be transmitting more power than required for certain 
bit rates. In this case, the loop bandwidth can be modified to provide 
better performance. The AIRS will take advantage of the information 
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Table 3-1. 

AIRS Carrier Loop 

Bandwidth. 

IAQ Symbol 

Rate (Ksps)* 

Loop Order 

Loop Bandwidth 


SR < 1 

1 < SR < 18.3 
SR > 18.3 


Third 
Third 
Thi rd 


23 Hz 

(.023) x SR 
420 Hz 


ACQUISITION 


I4Q Symbol Rate (Ksps)* 


Loop Order 


Lood Bandwidth 



SR < 

3 

Second 

90 Hz 

3 < SR < 

20 

Second 

(.029) x SR 

SR > 

20 

Second 

5 kHz 

NRZ only. 

Half for 

Manchester-coded symbol 



* 


'<jCinC^c 


om 


provided by the signal strength estimator for this purpose. 
3.4.3 Bit Sync 

The bit sync bandwidths can be selected as follows: 



Symbol Rate, 1/T S B L T S (*) . 

< 10 Ksps 1.5 ! 

1 

10 - 250 Ksps .78 | 

0.25 - 1 Msps .39 | 

1 - 6 Msps .2 | 

i 


I 

\ 

} 

I 

i 

{ 

! 

i 


dCinCo m 


- 52 - 



* 


'cjCinC^c 


om 


4.0 AIRS AND WDD COMPARISON 


To a certain degree, the wide dyanmics demodulator (WDD) is a 
predecessor of AIKS. This is exemplified by its digital implementation 
and its use of the data-aided loop for carrier and data clock 
recovery. In this section, the AIRS and WDD are compared in terms of 
their functional characteristics and implementation techniques. This 
serves to highlight some of the improvements that are designed into the 
AIRS. In addition, some unique features of AIRS are given. 

4.1 Functional Characteristics Comparison 

Figure 4.1 summarizes the differences between the AIRS and the 
WDD. The WDD is designed to replace the low rate demodulator (LRD) and 
therefore is basically a low data rate PN/BPSK demodulator. The bit 
sync and data detection capability is a byproduc' of its data-aided 
carrier recovery loop. The AIRS, on the other hand, is a complete SSA 
return receive system that replaces the LRD, the medium rate demodulator 
(MRD), and the low and medium rate bit syncs (LRBS.MRBS). Therefore, it 
is a dual channel PN/BPSK and UQPSK demodi ator and bit sync. It also 
interfaces directly with the Viterbi decoder (VD) to form an integrated 
receive system. 

The WDD is chartered to provide wide dynamics operation. The 
technique employed is to use a third-order loop for carrier recovery. 

The AIRS also employs a third-order loop. In addition, it performs its 
own Doppler compensation, so that in a typical operational scenario, the 
Doppler rate uncertainty can be reduced significantly. This further 
enhances its dynamics handling capability. 

The WDD reacquisition is basically the same as its acquisition 
except that the search range is reduced. The AIRS uses a multistep 
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• REPLACES LRD 


• REPLACES LRD, MRD, LRBS, MRBS 


• UP TO 300 Kbps • UP TO 3 Mbps 

• WIDE DYNAMICS CAPABILITY t IMPROVED WIDE DYNAMICS 

CAPABILITY 

• REACQUISITION MODE « IMPROVED REACQUISITION 

• SINGLE STEP • MULTISTEP 

• ACQUISITION PERFORMANCE • IMPROVED ACQUISITION 

• 15 SEC AT THRESHOLDS • 1 to 4 SECS 

ff FIXED LOOP BANDWIDTHS • VARIABLE LOOP BANDWIDTHS 


Figure 4.1. Functional Characteristics Comparison Between AIRS 


and WDD. 
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reacquisition strategy to adapt the searching algorithm to the diffusion 
of uncertainties as a function of the time from loss of lock. 

The acquisition performance for the WDD is comparable to (slightly 
worse than) the LRD. The AIRS uses a parallel searching technique 
implemented with CCD to speed up the acquisition time at low signal 
levels. 

Finally, the AIRS uses adaptive loop bandwidths to handle received 
signal level variations. 

4.2 Implementation Techniques Comparison 

A comparison of implementation techniques used by the AIRS and the 
WDD is shown in Table 4-1. 

4.3 Unique AIRS Features 

Features unique to AIRS are shown in Figure 4.2. These features 
are not found in the WDD. 
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WDD 


AIRS 


a PN ACQUISITION • 

• SINGLE-CHANNEL SEARCH 

• SEQUENTIAL DETECTION 
REQUIRES AGC CAL 

• FALSE LOCK PROTECTION 
WITH VARIABLE THRESHOLD 
TECHNIQUE DEGRADES 
THRESHOLD PERFOMANCE 

• PN TRACKING • 

• ONLY ONE CHANNEL IS 
TRACKED 

• MODIFIED E/L DESPREADER 

• FIXED BANDWIDTH 

• CARRIER ACQUISITION • 

• ONE CHANNEL 

• SWEEP ACQUISITION REQUIRES 
FALSE LOCK REJECTION 
ALGORITHM THAT CAN BE TIME 
CONSUMING 

• CARRIER TRACKING • 

• ONE CHANNEL 

a SECOND/THIRD ORDER 
LOOP 

a FIXED BANDWIDTH 


PN ACQUISITION 

• MULTI-CHANNEL SEARCH 

• FOR 4-CHANNEL SEARCH ONLY; 
NOT REQUIRED WITH CCD 
IMPLEMENTATION 

• MULTI-CHANNEL OPERATIONS 
ALLOWS FALSE LOCK 
VERIFICATION AFTER CODE HIT 

PN TRACKING 

• UP TO TWO CHANNELS ARE 
TRACKED (3 dB IMPROVEMENT) 

• DOUBLE TAU-DITHERED LOOP 

• VARIABLE BANDWPTH 

CARRIER ACQUISITION 

• BOTH CHANNELS 

• FLL IMPLEMENTATION CANNOT 
FALSE LOCK 


CARRIER TRACKING 

• BOTH CHANNELS 

• SECOND/THIRD ORDER LOOP 
AND DOPPLER COMPENSATION 

• VARIABLE BANDWIDTH 


• IF SAMPLING AT 5 MHz • 

• REACQUISITION (ONE STEP) • 

• FREEZES PN LOOP AND 
PERFORMS CARRIER REACQ 


BASEBAND SAMPLING 
REACQUISITION (MULTISTEP) 

• DISTINGUISHES BETWEEN PN, 
CARRIER AND BIT SYNC LOSS 
OF LOCK COMBINATIONS 

• WAIT, MINISEARCH & COLD START 


Table 4.1. Implementation Comparison of WDD and AIRS. 
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• INTEGRATED DEMOD/BIT SYNC/DECODER OPERATION 

• AUTONOMOUS STAND ALONE OPERATION 

• PROVISION FOR RFI MITIGATION 

• DOPPLER COMPENSATION 

• AUTONOMOUS (NORMAL) /FLEXIBLE DATA/TEST MODE 
t I/O CHANNEL AMBIGUITY RESOLUTION 

• VARIABLE LOOP BANDWIDTHS ALLOWS REAL TIME PHASE NOISE/ 
THERMAL NOISE TRADEOFFS 

Figure 4.2. Unique AIRS Features. 
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5.0 AIRS ANALYTICAL/SIMULATION SOFTWARE 


The primary purpose of the AIRS simulation software is to provide a 
computer emulation of the AIRS hardware and the AIRS software 
algorithms. The simulation is a tool for checking out and verifying the 
soundness of the technique used to implement the AIRS, and for exposing 
bugs and problems in the early stage of the design. The software tool 
is also used to evaluate and predict the performance of various AIRS 
subsystems. Finally, the software tool can be used to select the 
various system parameters whenever a tradeoff exists. 

Since .he algorithms developed in the software simulation can be 
transferred to the AIRS microprocessor software with minor modifications 
(perhaps to accomodate the microprocessor computation cycles), a strict 
Monte Carlo technique is adopted whenever possible. In the cases where 
the Monte Carlo technique is impractical and/or overly time consuming, 
the software is supplemented with analysis. 

Figure 5.1 shows a typical application of the AIRS software. In 
this case, the simulation is used to optimize a design parameter. 

5.1 AIRS Simulation Software Structure 

The AIRS simulation software package is designed to evaluate the 
performance of the AIRS in terms of acquisition, tracking, BER and 
reacquisition. The software is organized into the following modules: 

a. Interface module 

b. Configuration module 

c. Performance module 

The interface module is used to translate the output of a CLASS 
execution module, whose purpose is to .define the link conditions and 
user characteristics such as data group and modes, into a form usable by 
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the succeeding configuration module. The configuration module is used 
to define the setup of the AIRS as a function of the link and user data 
characteristics, under the criterions established by the selected mode 
of receiver operation. Finally, the performance module computes and 
displays the various desired performance measures. 

5.1.1 Interface Module 

The pertinent output from the CLASS module defines the data 
characteristics for the I and Q components as well as the received 
effective C/Ng. The interface module translates the CLASS output into 
an explicit description of the I and Q channels. They describe: 

a. Data format (NRZ/bi -phase) 
d. Data rate 

c. Code rate 

d. PN spreading 


e. BPSK/QPSK 


i. Power ratio 
g. Effective C/Ng 

A summary of the I and Q channel characteristics are displayed as part 
of the interface module output. 

5.1.2 Configuration Module 

The configuration module takes the output of the interface module, 
and based on the selected mode of receiver operation, configure the AIRS 
signal porcessing architecture. The output of the module is a summary 
of the AIRS configuration for: 

a. PN acquisition 

b. PN tracking 

c. Frequency acquisition 
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d. Carrier tracking 

e. Bit sync 

as well as the associated subsystem parameters. 

5.1.3 Performance Module 

The performance module is a collection of software packages that 
are used to predict ARIS subsystem performance. The module uses a 
combination of analysis and Monte Carlo simulation as shown in Figure 
5.2 .he performance packages that are analytical in nature are: 

a. PN acquisition package 

b. PN tracking package 

c. BER package 

d. Decoder package 

the performance packages that are Monte Carlo in nature are: 

a. Frequency acquisition package 

b. Carrier tracking package 

c. Bit sync tracking package 

d. Carrier/bit sync tracking package 

These pack, es can be used to support four main executive programs that 
determine AIRS sytem acquisition, reacquisition, tracking, and BER 
performance. The interconnections between the modules and the main 
programs a< e shown in Figure 5.3. The packages all employ a user- 
friendly interactive approach to help the user define the pertinent 
subsystem parameters. 

5. 1.3.1 Analytical Packages 

5. 1.3. 1.1 PN Acquisition Package 

The PN acquisition packer determines both tn^ mean time to acquire 
and the time to acquire wi:h .i k probability for three PN acauisition 
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Performance Module Archi tecture. 
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schemes. Two schemes are based on using CCD PNMFs and one scheme is 
based on multi-channel sequential search. 

5. 1.3. 1.2 PN Tracking Package 

This package determines the initial pull-in, tracking, and slip 
performance of the configured subsystem. The output performance 
measures are: pull-in time, rms PN chip jitter, and mean time to lose 

lock. 


5. 1.3. 1.3 BER Package 

This package determines the bit error rate (BER) for uncoded data 
or the symbol error rate (SER) for coded data. For coded data the BER 
package also computes the probability mass at each quantization level. 
There are two portions of the package. The first one deals with the 
transient BER, i.e., during acquisition and high dynamics operation. 

For the transient BER computation, the effects of the mean carrier phase 
error are incorporated. The second portion deals with steady state 
operation. For the steady state BER computation, the effects of phase 
jitter and timing jitter are considered. 

5. 1.3. 1.4 Decoder Package 

This package computes the decoded BER for convolutional ly encoded 
data. It is based on computing the computational cutoff rate of the 
quant i Ted symbol. 

5. 1.3. 2 Monte Carlo Packages 

The main purpose of the Monte Carlo packages is to duplicate the 
AIRS hardwai e/software structure, so that it can be checked for 
potential problems. As such, the Monte Carlo packages should not be 
used for extensive performance predictions as the computational time 
will be prohibitively long. 
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5. 1.3. 2.1 Frequency Acquisition Package 

This package duplicates the AIRS algorithms. The received signal 
is sampled and processed in the manner of a frequency-locked loop. The 
algorithm stops once the incoming signal carrier frequency can be 
matched to within 25 Hz. This is done by implementing a lock 
detector. The output of this package is the mean time to acquire 
frequency lock. 

5. 1.3. 2. 2 Carrier Tracking Package 

This package is used in the Flexible Data Format mode in which the 

carrier loop and the bit sync loop must be decoupled. The output of 

this package is the mean time to acquire phase lock, the mean phase 
offset (which will be a function of time during transient), the rms 
phase jitter. The mean time to lose lock determination is partly 
analysis and partly simulation. 

5. 1.3. 2. 3 Bit Sync Tracking Package 

This package is used in the Flexible Data Format mode in which the 

carrier loop and bit sync loop must be decoupled. The output of this 

package is the mean time to acquire timing, the mean timing error, the 
rms time jitter, and the mean time to lose lock. The mean time to lose 
lock determination is partly analysis and partly simulation. 

5. 1.3. 2. 4 Carrier/Bit Sync Tracking Package 

This package is used during normal operation in which a data aided 
tracking loop configuration is used. In this case, the bit sync and the 
carrier loop acquires simultaneously. The output is the mean phase and 
lock acquisition time, the mean phase and clock error, the rms phase and 
clock error, and the mean time to lose carrier and bit sy” lock. The 
mear time to lose lock determination is a combination o* ’ysis and 
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Monte Carlo simulation. 

5.2 An Application Example 

Figuare 5.4 shows an application example of using the AIRS software 
to determine system acquisition time. Notice that the analytical 
portion of the software computes the time spent during code search, code 
loop acquisition, and Viterbi decoder node sync. The time spent during 
frequency, phase and bit sync acquisition are determined via Monte Carlo 


simulation. 











APPENDIX A 


RAPID PN ACQUISITION SCHEMES FOR AIRS: 
TECHNIQUE, ANALYSIS AND PERFORMANCE SIMULATION 


This appendix is a revised and updated version of an earlier technical 
report, TR-0484-8211, prepared by Y. T. Su and C. M. Chie under this 


contract. 
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1 . Introduction and Summary 

1.1 Introduction and Objective 

An Autonomous Integrated Receive System (AIRS) has been proposed by 
LinCom to improve the current Tracking and Data Relay Satellite System 
(TDRSS), ard to use as a receiving system for the future Tracking and 
Data Acquisition System (TDAS). The object of this study is to examine 
in depth three possible code acquisition subsystems for the AIRS, namely, 

• serial search using charged-coupled device (CCD) pseudonoise 
matched filter (PNMF) 

• parallel search using PNMF 

• four-channel parallel search using sequential detect 'on, 

and to develope a software analysis tool for predicting the performances 
of these systems. 

The software package is to have many useful features, for example, 
it can 

o generate the operation characteristics for different detector 
cc nfigurations, parameterized by SNR and integration time 

• produce sensitivity curves (or data), e.g., detection (false 
alarm) probability vs. threshold, noise power, SNR, etc. 

• evaluate the system performance (with given design parameters) 
under different SNR environments 

• optimize system parameters. 

1 .2 Organization of the Report 

This report is organized as follows. 

Section 2 defines the general PN acquisition problem and classifies 
search algorithms. Section 3 describes and motivates the three schemes 
to be investigated. Section 4 identifies assumptions needed for 
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analysis and summarizes performances of various systems. Conclusions 
and some suggestions for further studies are contained in Section 4. 
Appendices A to G detail the performance analysis and software 
appl ications. 

An executive summary is contained below in section 1.3. 

1.3 Executive Summary of Progress and Results 

We have developed analytical models for the three acquisition 
subsystems mentioned in Section 1.1. Complete performance measures, 
e.g., mean acquisition time, time for 90% acquisition, and a 
corresponding software package (as one part of the AIRS demonstration 
program) which contains the features listed in Section \ .1 are 
obtained. In the course of analyzing system performance, we have 
developed analytical expressions, both exact and approximation. For 
detection and false alarm probabilities of each detector structure 
involved. Efficient numerical algorithms for computations are also 
found. To assess, according to merit, the three candidate systems, 
numerical optimizations have been performed. As expected, the CCD 
receivers outperform the sequential detection method. For the two 
receivers using CCD PNMF, the serial search is superior to the parallel 
search though this advantage dwindles as the reference C/Nq increases. 

The AIRS' current acquisition requirements are summarized in Table 
1-1. This study has found that these requirements can all be met by 
both CCD systems without difficulty, as exemplified in Table 1-2. This 
table is evaluated with the following parameters: 


Code Epoch Timing Uncertainty 
Reference C/Nq for One Channel 
Number of Doppler Bins 


±1900 chips 
36.5 dB-Hz 
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Table 1-1, 


Acquisition Time vs. C/N n and Mode Requirements 
for the AIRS 


ACQUISITION TIMES FOR % ACQUISITION 


C/N 0 dB-Hz 

5 SECONDS 

15 SECONDS 

45 SECONDS 

39.5 

N/A 

DG-1 Modes 1, 2A, 3 
DG-2 Modes 2A, 3 

DG-1 Mode 2B 
DG-2 Mode 2B 

49.5 

DG-1 Modes 1, 2A, 3 

DG-1 Mode 2B 

N/A 


Table 1-2. Acquisition Performance for Algorithms 
Using CCD PNMF 


-) 


Algorithm 

Parallel Search 

Serial Search 

Acquisition Time (sec.) 

3.95 

3.89 

Acquisition Probability 

0.914 

0.925 

Number of CCD PNMF's 
Per Channel 

8 

8 

Length of a PNMF (cells) 

255 

255 

Required Memory Size (bits) 

2 io 

? 10 

for Noncoherent Integration 

( - 8K) 

(- 8K) 
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0.4 dB 
0.2 dB 
3M chips/sec 
6 MHz 
6 MHz 
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Sampling Loss 
Doppler Loss 
Data Transition Loss 
PN Code Rate 
IF Bandwidth 

CCD Transport (Sampling) Clock Rate 
The inclusion ot the 4-channel sequential detection method serves 
to illustrate the superiority of the CCD receivers. For example, with 
the same se* of parameters shown above, the minimum mean acquisition 
times (in sec) for the CCD receivers are 2.47 (serial) and 3.95 
(parallel) respectively while that for 4-channel sequential detection is 
at least 6.65. that the minimum mean acquisition times for the 

C-D receivers art . . "erent from the acquisition times in Table 1-2. 

Trie reason is that the data obtained in Table i-2 is based upon the 
criterion of minimizing the time required for 90% acquisition. 

The above numerical results are shown in part to highlight che 
versatility of the software analysis tool developed by this study. Not 
only the goals set up in Sec. 1.1 are accomplished but several more 
functions are added. Among these additions are the capabilities to: 

• Obtain the optimal system parameters under the criterion of 
either minimizing the mean acquisition time or the time for 90% 
acquisition. 

• Estimate the correlation loss due to pre-despreading filtering, 
t Help to decide between performance and complexity (i.e., number 

of channels, length of CCD PNMF's, bin width, a f •< bit numbers). 
Details can be found in Appendix G. 
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2. Background 

The PN acquisition problem can be defined as to find a pair of 
Doppler and code phase estimate (u> d , t) for the incoming waveform: 


Y (t ) = /2S d(t) P n (t-r) cos[(u» c +a) d )t+(|)] +n(t) (2-1) 


such that 


|(i) d -ai d | <5u), (2-2a) 

and j x — x j -C6 T c . (2-2b) 

In the above expression, P n (t-x) is the received PN code with rate 
T c ~l, delayed by t with respect to an arbitrary time reference, d(t) is 
the data modulation, n(t) is additive white Gaussian noise with a one- 
sided power spectral density Nq. S, <d c and $ are the carrier's power, 
frequency and phase respectively and is the Doppler frequency. The 
requirement of simultaneous Doppler estimate u» d is due to t!v' fact that 
the Doppler offset must be small enough not to incur too much signal 
power degradation when correlation between received signal and locai 
estimate is performed. Hence, the choice of 6u> depends not only on the 
pull-in region of the carrier-phase tracking loop but on the Doppler 
loss tolerable in obtaining code phase estimate x. The parameter 6T C , 
on t 1 other hand, is typically chosen as T c or a fraction of T c , which 
is within the pull-in capacity of the code-tracking loop. 

Consider the rectangular time-frequency uncertainty region shown in 
Fig. 2-1, where a small 6o> x 6T C rectangular is defined as a cell. The 
PN acquisition problem is now equivalent to find out which cell 
satisfies the condition (2-2). Here, there may be more than one such a 
ceil because usually we test only one specific point within a cell; see 
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Fig. 2 -1. A Time- Frequency Uncertainty Region. 


AN ENLARGED CELL (OR A SECTION) SEARCH DIRECTION 
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also Fig. 2-1. Thus, acquisition algorithms are in fact searching 
strategies and can be classified into three categories: serial, 

parallel or compound search. This report discusses two algorithms which 
fall into the third category and will be called as SPS (serial-parallel- 
serial) and SSP (serial-serial-parallel) algorithms according to the wa^ 
signal is processed. The compound nature of these algorithms, as we 
shall see, is due to the usage of the so-called CCD PN match filters 
(PNMF) [1] and the low chip-signal energy to noise power density ratio 
(E c /Ng) environments. Algorithm similar to SPS has been discussed [2], 
though no rigorous definition and analysis are given. Pure serial or 
parallel search algorithms using either digital match filter [3] or SAW 
device [4] were also investigated recently. 

3. Acquisition Algorithms Using CCD PNMF 
3.1 An Analytical Model for CCD Detector 

An equivalent CCD detector is shown in Fig. 3-1. The received 
waveform, after frequency deconversion (mixing), is integrated for a 
fraction of the chip time and sampled at a rate R S =1/6T C (integrate-and- 
dump) in both I and Q arms. The sampled taseband I & Q signals are then 
fed into an I/Q pair of CCD PNMF's (correlation), squared in each arm 
and finally added to yield a tet statistic R(oTq) at time t=nTQ, 
where T Q - 6T.. 

Denote the sampled skoals at the input of I and Q channel PNMF's 
by yj(kTq) and yq(kTq) respectively. Then the oucputs of I/Q PNMF's are 

M 

R a (nT 0 ) = ^ P n (,cT 0 ) Y q C(<+n-M)T 0 ], M>n>M (3-1) 

y <=l y 

where g denotes I or Q (channel). 
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Obviously, with the capability of retaining thus reusing the 
received information, y(<Tq)'s, the CCO detector can produce the test 
statistics, R(nTg)'s, at a rate M times faster than conventional analog 
I-Q detector. However, this improvement factor M can not be made 
arbitrary large due mainly to the presence of the data modulation 
d(t). On the other hand, in small E C /N Q situations, long correlation 
period is necessary to get a robust statistic. This dilemma can be 
solved by resorting to noncoherent integrations (NI). NI, or post- 
detection integration, is a technique borrowed from pulse radar 
detection theory [6] referring to the case in which phase coherence 
between consecutive pulses is not maintained. In our case, NI is 
equivalent to the following summation operation: 


where 


"I 


R i = R x 2 (i) + R Q 2 ( i ) 


(3-2) 

(3-3) 


R(i) = l P n C(i-l)MT 0 + <T ] y [(i-l)uT 0 +<T 0 ] 
y <=1 y 

iMT n 

= / P n (t'h a (f + (i-lUT n )df 
(i-l)MT n 9 0 


(3-4) 


u is a constant larger than or equal to M and aTq = (u-M)Tq. 

3.2 Algorithm Using CCD PNMF 

We are now in a position to define the SPS and SSP algorithms. 
Basically, they are all two-stage tens. The difference lies in the 
first stage test only. The second stage test, or the verification mode 
is used to reduce the false alarm rate so as to avoid a costly false 
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code tracking attempt (usually modeled as a penalty time [7]). A 
single-dwell test (SDT) will be used as the second stage test for both 
algorithms. As we shall see, the acquisition time is insensitive to 
that consumed by a properly designed SDT for the second stage. 

(A) The SSP algorithm: 

(A.O) Load the PNMF with M samples of Yg(t) (initialization). 

(A.l) Fix a M-cell segment of the PN code for local reference. 

(A. 2) Let the incoming waveform samples run through the PNMF for 
MTq seconds and collect serially M test statistics R^-, (i-1) M 
+ 1 <; j < iM, corresponding to a section of M cells (Fig. 2-1) 
in the uncertainty range. These numbers are stored in a first 
in first out (FIFO) register for further (noncoherent) addition; 
See Fig. 3.2. 

(A. 3) Repeat (A.l) & (A. 2) for Nj mutually disjoint segments of the 
local PN code. The test statistics corresponding to the same 
code phase are added (i.e. NI) and at the end of the Nj th 
iteration, only the largest* and its phase are kept in two 
registers. 

(A. 4) Repeat the process (A.1)-(A.3) for all (disjoint) sections 
(or bins) in the uncertainty range. 

(A. 5) The code phase with the largest test statistic is to be 
tested by the second stage test and will be used as our 
estimate . (see eq. (1-2)) for further code phase tracking if 
passes; otherwise, go back to (A.O) and start all over. 

The SSP algorithm is based on the maximum likelihood estimate principle 


*Here, the largest is chosen from the M cells in the section under 
^••'nation and the previous largest one. 
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to use the code phase which yields the highest correlation value. Such 
an algorithm is optimal in the sense of minimizing the error probability 
[8]. Nevertheless, our criterion is to minimize the mean acquisition 
instead and hence the cost is time not error rate though they are 
related. The second algorithm, taking the above fact into 
consideration, is essentially the same as the well-known serial search 
algorithms [7] with enlarged cells, where an enlarged cell is what we 
call a section (bin) in (A. 2) above, i.e., M consecutive cel?s; see also 
Fig. 2-1. Therefore, 

(B) The SPS algorithm (Fig. 3.2): 

(B.O) = (A.O) 

(B.l) = (A.l) 

(B.2) = (A. 2) 

(B.3) Repeat (B.l) and (B.2) for Nj iterations. If the largest 
test statistic is larger than a threshold then go to the 2nd 
stage test, otherwise go to the next section and repeat this 
process (B.l) to (B.3). 

The purpose of this algorithm is to shorten the time needed to go to the 
second stage if the resulting statistics of a certain section indicate 
strongly that a good estimate is obtained. 

3.3 Four -Channel Parallel Search Using Sequential Detection 

This approach divides the uncertainty range into four equal length 
regions. Each region is searched by one channel. Once coarse 
synchronization is declared in one of the four channels the other three 
sequential detectors w'll stop searching and the system then enters the 
code tracking mode. Theoretically, the acquisition time will be reduced 
to only a quarter of that of a single channel case, assuming very low 
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false alarm rate, since the Hj-cell exists in only one region. Thus it 
suffices to analyze the one-channel sequential detection. An analytical 
model of the sequential detection used for PN acquisition is depicted in 
Fig. 3.3. The product of the local code and the incoming waveform is 
fitered by a bandpass filter which is wide enough to accomodate the data 
modulation and the Doppler uncertainty. A sequence of independent and 
identical distributed (i.i.d.) observations are then obtained by 
sampling the envelope detector output at a rate = 1/28 jp where B j p is 
the bandwidth of the BPF. These observations, , are further modified 
by first multiplied by a and then subtracted by b, where a,b are 
constants depending on E c /Nq. The cumulated sum of modified 
observations 


i .e.. 


X. $ aR.-b , 



» 


are fed into the decision logic. There are two decision logics that can 
be simply implemented. The first one, denoted by S^A,B), uses the 
decision rule: 

> A stop sampling and decide that is true (D^) 

e(A,B) continue by taking another sample (D) 

< B stop sampling and decide that is true (D Q ) 

while the second one, denoted by Sg ( B , N t ) , uses the following rule: 





with Sequential Detection. 
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> B for all 1 < k < 


S. < 8 for any 1 < k < N. 


(+f 


i 


otherwise 


S2 ( B , N t ) is usually a preferable choice because: 

(1) The truncation time is needed to set an upper bound for the 
per cell dwell time so as not to reduce the Doppler range 
over which a coarse synchronization can be detected. 

(2) The acquisition time is insensitive to the dwell time spent 
in Hj-cells. After all, the number of Hg-cells is hundreds 
or thousands times as large as that of Hj-cells. 

Except for the mean dwell time in a H^-cell, the two tests are 
practically the same as far as the performance is concerned. 

4. Performance Analysis and Numerical Results 
4. Summary of Performance Results 

In Appendix A, we show that the test statistics R(nTg) produced 
serially by the CCD correlator are independent if 6 = 1 and exhibit only 
very small correlations among neighboring R(nT 0 )'s if 6 < 1. 

Furthermore, if we define a H^-cell as one which satisfies (2-2) and a 
HQ-cell which does not, Hj, the hypothesis that the present cell under 
examination is u H^-cell and Hq its complement hypothesis. Then, as 
mentioned in the beginning section, there are more than one H^-cell in 
most cases. Moreover, the number of Hg-cells is much larger than that 
of H^-cells. To facilitate our analysis, we shall make the folic mg 
aesumptions: (1) all test statistics are independent no matter what 6 

is. (2) there is one and only one H^-cell in the uncertainty region, (3) 




I 

I 

I 

i. 


1 






the code self noise* is negligible, ^nd (4) M( in chins) >> 1. 

Under these assumptions, it can then be shown (Appendix B and D) 
that the detection and false alarm probabilities, Pg^, Pp^, for the 
first stage test are 


P 0 ; / PtyjHjJdy 

1 Y 

P FA. • W P M V dy 

1 Y 


(for SPS) 


(4-la) 


(4-lb) 


P D ■ J P(y| H !)[/ P(x|H 0 )dx] °dy 


(for SSP) 


PpA l = 1 - % 


(4-2a) 


(4-2b) 


respectively, 


where 


P(y | H l) = , (4-3a) 


p (y| H 0 ) = 


(4-3b) 


N = the number of noncoherent integrations, 

3 = E c ' N 0’ 

^ N- 1 ( * ) = mod if' ecl Bessel function of order N-l, 

Cg = the nubmer of Hg-cells. 

Detailed analyses and computational algorithms pertaining to P and 
Pp^ are to be found in Appendix B and D. As for the 2nd st^ge test, 
Aprvndix C contains the derivation of the operation characteristic, P 
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•s. Pp^. The mean acquisition times for both algorithms are (see 
Appendix F) 


T ACQ “ t D + (p D " I^ x D 1 +P FA 1 T D 2 +P FA' T r ‘Wt’ ^ SPS ^ 


and 


t acq ■ t d * p d 1[n ' c No +1 > * ^ 2 p fa/ t c + p fa- t p ]t c ‘ SSP > ( 4 - 5 > 


where 


P D = P D 1 * V 
P FA " P FA X * p FA 2 » 

T D ■ MT c> 

T Dj = Nt D* 

= dwell time for the second stage est, 
,p = Penalty time (in chips), 

N d = number of Doppler bins, 

N t = number of time bins. 


The analysis details of the sequential detection method are contained in 
Appendix E. 

4.2 Numerical Results 

Based upon the above analysis, the numerical behaviors of both CCD 
algorithms are demonst ' ;ed in Fig. d-1 through Fig. 4-5. An empirical 
optimization rule is adopted: Pp^ must be small enouyh to offset Tp. 

According to this rule, we set Pp^ = 10"^ such that (Pp^ • Tp) T . is of 

d 

the order of msec even with Tp as large as } n ‘ . The other system 

parameters: P D ^, Pq^, Pp/^> t Dj* t D ? * un ^ ess specified, are chosen to 

minimize T^.g. 

Fig. 4-1 shows the SSP performance as a function of N- and M 
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(chips) when E c /Nq = -30.5 dB, 5T C = .5T C and p ^ |t-t|/T c = .25. An 
optimal Nj do exist for each parameter set. The accumulated acquisition 
probability FacqU) is shown in Fig. 4-2, where we use the formula: 


F ACQ W = Pr £ T ACQ U ^ 


where 


A n 


T 0, 


[t/T A ] 

1 

I (I'Pd)" P D» 

(4-6) 

r=l 

I 

FA. + P FA*V T c’ 

1 

(4-7) 

1 


[x] = the largest integer less than x. 


Comparisons can be made: between the best case (p=0) and the worst caes 
(p=.25), among different M (coherent integration period). The choice of 
M's is of course dependent on the data rate. 

Performances for the SSP and SPS algorithms as a function of E c /Nq 
and Nj are described in Fig. 4-i and 4-4, which exhibit the importance 
of a judicious choice of Nj as the E c /Ng varies. This is because Nj has 
put an inherent lower bound on the time required to make a decision (see 
the flat portions of the curves corresponding to different Nj's in Fig. 
4-3 and 4-4). 

Comparison between the ^PS and SSP algorithms is made in Fig. 4-5 
where Nj has been optimized. It can be seen that the SPS algorithm is 
superior to the SSP algorithm though this superiority is dwindled as 
E c /N 0 gets higher. The SSP algorithm, however, is less sensitive to the 
variation of E c /Ng owing to the fact that a threshold for the first 
stage test is not necessary; see Fig. 4-7 and Fig. 4-8. 
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The relative insensitivity of to can be found in Table 4-1 
to 4-3 for both algorithms, where t(0.9) is defined as the time required 

to acquire with probability 0.9. Note that this number for the SPS is 

an upperbound; see Appendix F. 

The above numerical results are computed with the following 
parameters: 

N d (Doppler bin) 3 

1+Cq (cells per Doppler bin) 7600 

6T C (step size) 0.5 chip 

Sampling loss 2.5 dB ( p= .25 ) 


Doppler loss 
Data transition loss 
Code chip rate 


IF bandpass filter bandwidth 
Processing gain 


0.5 chip 
2.5 dB ( p= .25 ) 
0.4 dB 
0.2 dB 
3 MHz 

512 (chips) 

6 MHz 


30 dB 


Evidently, even with C/Nq = 36.5 dB-Hz (only one channel is used), i.e., 
E c /N 0 = -31.9 dB (sampling loss excluded), we still can acquire within 4 
sec with a probability equal to 0.9. The AIRS' acquisition requirement 
thus can be met without difficulty. For comparison purposes we include 
here the best sequential test performance (Fig. 4-6). The normalized 
mean acquisition in Fig. 4-6 can be converted into in seconds via 


T ACQ^ sec ^ = V 1+C o >V 


(4-8) 


5 . Conclusion and Further Sti ies 

The main results achieved by this study can be summarized as 
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Table 4-1.* System Performance Parameers for the SPS. 




C/N 0 (dB-Hz) 

36.5 

37.5 

38.5 

/ 

39.5 

40.5 

^ACQ ^ SeC ^ 

2.58 

1.65 

1.05 

0.722 

.494 

tu 2 (sec) 

1.4(-1) 

9.2(-2) 

6.2(-2) 

4.2(-2) 

3.53(-<?) 

N I 

700 

500 

400 

300 

200 

S 

.743 

.794 

.908 

.948 

.9o0 

P FAi 

3.64(-2) 

2.94(-2) 

3.23(-2) 

2.78(-2) 

4.bo.-2) 


.990 

.991 

.992 

.993 

.990 


f 

i 

i 

I 


i 

i 


Table 4-2.* 

System Performance 

Parameters 

for the SSP. 

C/Nq (dB-Hz) 

36.5 

37.5 

38.5 

39.5 

40.5 

T ACQ^ sec ^ 

3.88 

2.53 

1.65 

1.09 

.737 

^ (sec) 

8.3(-2) 

5.9(-2) 

3.9(-2) 

2.6(-2) 

2.33(-2) 

N I 

700 

500 

2 'JO 

200 

150 

S 

.712 

.778 

.719 

.728 

.811 

P FAi 

.288 

.221 

.281 

.272 

.189 

S 

.995 

.995 

.995 

.994 

.993 


*R C * 3 MHz, M = S12 (chips), BWjp * 6 MHz, processing gair« = 30 dB. 
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Table 4-3.* 

Performance Parameters 
Alyorithms when C/Nq = 

for the SPS and SSP 
36.5 dB-Hz. 

Algorithm 

SSP 

oPS 

T ACQ ^ sec ) 

4.27 

2.56 

T|, (sec.) 
u 2 

1.0"' -1) 

1 .67 ( -1 ) 

Ni 

1000 

950 

t(0.9) (sec. 

) 3.91 

3.89 

p o 

.915 

.927 

Pfa, 

8.19(-2) 

4.90(-2) 


I ! 


r.. air- -v 


(h. i i 
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follows. 

(I) Precise definitions of two algorithms for code acquisition 
subsystems using matched filter type correlators. 

(II) Rigorous performance analysis for both algorithms. This 
analysis can be easily applied to other algorithms, e.g., the 
fast local code algorithm described below or, the sequential 
test algorithm discussed in Appendix E. 

(III) The addition of the second stage test, -he improvement due 
to this addition is similar to that obtained by a two dwell 
system over a single dwell system. 

(IV) Numerical optimizations have been carried out for both 
algorithms. 

(V) Sensitivities of both algorithms are examined and performance 
comparisons are made. 

(VI) Development of a software analysis tool (Appendix G) with 
versatile capabilities. 

We have investigated thus far the behaviors of the SSP and SPS 
algorithms v. the region which features low E c /Nq, long code period and 
large or medium uncertainty ranges. This condition, i.e., the 
parameters chosen in Sec. 4, corresponds to the worst condition the AIRS 
is supposed to operate. Evidently, from a worst case consideration the 
SPS algorithm is a preferrable choice. 

As to the hardware aspect, it appears that a dual 256-tap CCD 
correlator will be aviiable this summer from Fairchild and Fo r d. The 
device is based on a mask developed by S. C. Munroe of the MIT Lincoln 
Laboratory [10]. It appears to be directly compatible with the AIRS 
acquisition requirements. Becasue of the size of the bulky substrate 
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required, SAW correlators are deemed to be unsuitable for the AIRS 
application [11], 

An improvement on the acquisition subsystem performance presented 
in Sec. 4 can be achieved by minor modification on our CCD algorithms. 
The riodifications are made in (A.l) and (A. 2) (i.e., (B.I) and (B .2) ) 
[9]. Instead of fixed local code reference segment during every 
coherent integration period, we let the local code run at a rate N+l 
times faster than the sampling rate. Hence, before a new sample 
entering the PNHF, N test statistics are produced. The N+l-th one 'S 
discarded because of the concurrent shift of the sampled signal and the 
local code. Judging from eq. (4-4), (4-5) and Tables (4-1), 

(4-2), T arn can be reduced by as much as N times so long as in /N » 
max(TQ,P FAi *x D2 ) . The achievable local code rate (N+1)/5T C is of course 
another technology issue. It appears now, for a sampling rate 1/5 T c = 

6M chips/sec, not much room (i.e., N) left to be improved. 

Recently, two systems similar to the SSP scheme have been 
reported. The CCD receiver described in [5] used a quantized NI 
method: every coherent integration output becomes either 1 or 0, 

depending on whether it passes a threshold or not. The other one 
proposed by Holm [7] used 64 bits digital correlators built by TRW along 
with the fast local code algorithm mentioned before. Because of the 
digital correlators, hard-limiting on the received waveform is 
necessary. Both systems are built to operate in a small uncertainty 
region and high (relatively) E c /Ng environment. We hope that the 
performance evaluations of these two systems can be accomplished in the 
future so that we can make assertions on the losses incurred by 
quantizations. 
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APPENDIX A 

AUTOCORRELATION FUNCTION OF THE CCD CORRELATOR OUTPUT 

We want to show, in this appendix, that the output samples of the 
CCD correlator are (almost) uncorrelated and independent if they are 
samples corresponding to the code phases in the Hq region. 

Because of the nonlinear operation in both I and Q channels, the 
uncor relatedness of the u g samples is difficult to prove directly from 
the statistic of R(t), the correlator output. However, if we can show 
that samples before the square-law device in both channels are 
independent then this independence can be carried over to the output, 
since the squaring operation is memory less. 

For analytical simplicity, we shall assume that (1)* there is no 
data modulation and (2) the sampling value at t = nTg, Tq = T c /m, in « I 
or 2. Then the output of the I-channel PWff is 

K 

Rj(nT 0 ) = cos(M>nT 0 +*) J {P n (kT 0 )[/SP n ((n-M+k)T 0 +r; 

+ n((n-M+k)T 0 +t)]AT 0 } (A.l) 


where ji| > T , n > M, S is the signal power and aTq is a constant 
related to t q . 

From assumption (2) and the uncorrelatedness between P n (kTQ)'s and 
n(kTg) 's, we have 


*As a matter of fact, this assumption is redundant unless |t| < T-, 
because a binary random data modulation sequence will not alter trie 
statistic of Rj or Rq if | x J 2 .T C ; see also (a. 4). 
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EtRjCnT^R^mTQ)^} 

= cos ( AajnT 0 +(j» )cos (AutfnT 0 + 4 » )aTq{ I [ SECP n ( kT Q ) P n («.Tq) P p ( ( n-M+k ) Tq+x ) 

k St 

* P n ( (m-M+i ) T q+t ) ] 

+ I I E[P n (kT 0 )P n (4T 0 )]ECn((n-M+k)T 0+ T)n((m-M+£)T 0 +x)]} 
k n 

- cos(AujnT 0 +<t»)cos(AcdmT 0 +«t.)R I (n,m) (A. 2) 

Hence, 

£{ R j ( nT Q ) R j (mT Q ) } ■ ^ cosCAufn-mjTQlR^n.m) (A. 3) 

Observing that 


E{n[r-MH)T 0 n]n[(m-M+*)T 0 +T]} 


and 


= 6 


n+k,m+Ji 


'1 

- 

,0 


if n+k = m+fc 
otherwise 


E[P n (kT 0 )P n UT 0 )] = 0 if | k-fc | > m 


(A.4) 


(A.5) 


we then conclude that the second double summation in the bracket becomes 
zero if I n-m | > in. 

In case that the PN sequence { P ( t ) } can be modeled as a random 


non 
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binary (±1) sequence, i.e., MT C « L(Period of fP n (t)}), then a similar 
conclusion for the first double summation can be reached. Furthermore, 
in the low E c /Ng envirorment, this term is negligible even if |n-m| 

< m. Looking back at (A.l), we can see that Rj(nTg) is actually a 
linear combination of Gaussian R.V.'s and is itself a Gaussian R.V. 

Independence of Rj(nTg) follows from the uncor relatedness. As mentioned 

2 2 
before, Rj(nTg) are then independent, so are Rg(nTg). The independence 

between I and Q channel is generally valid, so is that of 

Rj( nT 0 ) + Rg(r>T 0 ) $ R(nT 0 ). 

Notice that there still exists some "residual" correlations among 
m consecutive samples if in > 1, this residual correlation will be 
reduced by a factor of N, after N noncoherent integrations, where N » 1 
in all cases of interest. Therefore, for practicl purposes, we can 
neglect the correlations among these m (usually, m=2) consecutive 
samples. 
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APPENDIX B 


STATISTICS OF THE FIRST STAGE TEST 


We shall derive the detection (Pq^) and false alarm (Pp^) 
probabilities for the 1-st stage test of the SPS scheme, those of the 
SSP search will be discussed in Appendix D. Because 


S. 

J 


I 

l 

i=l 


R ij 


(B.l) 


where R. . is the sampled value obtained from the i-th coherent 

* J 

integration at the j-th cell, the statistics of R. . under both Hq and Hj 
must be obtained first. 

Fig. B.l shows an equivalent coherent integration circuit, where Rj 
and Rq can be written as 

Rl(HT c ) ■ x(*T c ) cos <|> + [x(MT c ) cos <j> + N'(MT c )] (B.2a) 
and 

R q ( mt c ) = x(MT c ) sin * + [x(MT c ) sin 4 + N‘(MT C )] (B.2b) 

where 

x(t) ft x(t) - E[x(t)] ft x(t) - 7(t) (B.3a) 

MT c 

x(t) ft (/5T ) _1 / PN(t+sT )PN(t+sT )dt (B.3b) 

c 0 c c 


and N'(MT r ), N'(MT r ) are i.i.d. Gaussian R.V.'s with zero mean and 
^ E c 

variance ^ -rp (note we have a normalization factor 1//5T in both 

c ”0 


correlators) . 

At low E c /Ng, the fluctuation part of x(t) is negligible, so for a 
fix 4 both Rj and Rq can be regarded as sum of a constant (possibly 
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Fig. B.l. An Equivalent Coherent I-Q System. 
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zero) and a Gaussian R.V. With this assumption, the density functions 
of R 2 ( t ) under both hypothesis are well-known, i.e., they are [1] 


and 


where 


P(R(t)|H ) = (2 o 2 )" 1 exp - M£)±!!L i r (t) ) 

2o o 


and 


Define 


P( R( t ) | Hq) = (2a 2 ) ^ exp - 

2o 


R(t) 


0 = 7 


m 


(l-|p|)t/T c 


p = s - s 


R(MT c )/2o 2 = Z 


(B.4a) 

(B.4b) 

(B.5a) 

(B.5b) 

(B.5c) 

(B.6) 


m 

17 


t=MT, 


- u-| p |) 2 (ir) n - 8 

0 


(B.7) 


then ( B .4a) and (B.4b) becomes 


P(Z|H 1 ) = e" (Z+e) I 0 (2/B?) 


P(Z|H 0 ) = e* 


(B.8a) 

(B.8b) 


It is now straightforward to obtain* [2] 


*P(y|H 0 ) is a version of Chi-square density function of 2N degree 
freedom and P(y|Hj) is the so-called noncentral Chi-square density with 
noncentrality parameter & and 2N degree of freedom. 
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MV = (N-l) ! 6 


and 

P(y l H l } = %^’ 1/2 e' (y+N6) I N _ 1 (2/N^‘) , 


where 

y = S./Zo 2 , 


N = Nj . 


Pp^ is then equal to 

p d ' I &) N ' 1/2e ' (y+N8>I N . 1 ( wl W) d y 
1 y 


- ,y ) (Generalized Q Function). 


Two possible false alarm probabilities must be distinguished 


Pp i - Pr[a false alarm occurs [H^ 


P F 0 ” Pr ^- a fa1se alarm occurs|H Q ] 


A . 

- one and only one samples in the section (M chips) 
under examination belongs to H n region. 


where 




H q - the whole section under examination belong to 
H q region. 

It is clear that 

P Fa * t h - 1 ) P FA‘ P f 

P F,o'sS P FA (8.12b) 

where 

Pf = P p {sample value y is the largest |y > threshold, ye H q} , 
and 

«» N-l 

p fa - i y hr 

* 1 - P(N,7) 

2 N-l 

= (1 + Y + Yp +•••+ ('N'_~i)l ) e Y (B.13) 

Eq. (B.10) and (B.13) are special cases of Eq. ( D .2b) and (D.2a) of 
Appendix 0 where computational aspect is discussed. However, (B.10) or 
(B.13) is much easier to evaluate, for example. Pari [3] has a very 
stable recursive method for computing QpjU.y). 

A simple lower bound for Pc i is 



" ~ L. 
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For large M, this becomes 


p * M n - p 

k F,1 31 FA F,0 


(B.14) 


Hence setting 


P ?k x = P F,0 = P F,1 


the resulting system performance, i.e., T , will be an upper bound, 


Finally, we like to mention that Pj^ in (B.10) is actually an upper 

1 


bound also. The "real" Pp_ is 


P Dj = Q N (B.Y)(1 - P F>1 ) 


However, Pp |_ « 1 is oftentime the case of interest, thus 


P D 1 = 


For large N and small 0, Pp^ and P Q ^ are related by [4] 


/N(i+2e)«f 1 (P Di ) = ✓Rfr” 1 (P FA ) - N0 


(B.15) 


(B.16) 


where 4 > (•) is the inverse complement error function. This relatin is 
obtained from Gaussian approximations of P(y | H q) and P(y|Hj). Hence a 
rough estimate of required N is 


✓TN e’WpA ) - /T?Zb ^(Pn )] 


(B.17) 

'dlinCo 
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■V * 
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ON THE DERIVATION OF THE RECEIVER OPERATING CHARACTERISTIC (ROC) 

OF THE SECOND STAGE TEST 

Consider the test shown in Fig. C.l, where we define 

MT c 

y(MT ) = / x 2 (t)dt , (C.l) 

c 0 

and x 2 (t) has a bandwidth of Wgp, the bandwidth of the BPF before the 
sque ing device. 

1 9 

If W << W BP* W BP MT c » 1 then x c (t) can be modeled as 

wideband process. Thus under Hq, it can be easily shown [1] that 

E[x 2 (t)|H Q ] = N 0 W Bp 

^ tn Q (C.2a) 

and 

Var[x 2 (t) |h q ] = (N 0 V.' Bp ) 2 



4 « 2 
" °0 * 

(C.2b) 

Hence 

MT 



C O 

E[y(MT c ) l H o3 = / ECx"(t)3dt = (N 0 W BP )MT C 

(C.3a) 


Var[y(MT c )|H 0 ] = N 2 W Bp MT c . 

(C.3b) 


• ] 



Under Hj, the mean and variance of y can be obtained in three 
steps. 

First, by [2] the power spectral d ?■ ;t v fu. ction of x(t)* (see 







Fig. C.l) is 


x(f) = (1-|P| ill ) 2 6(f) tilip 2 l sinc 2 (anp)6(f+ ^-) 


L+l 2 


ra= 
m*0 


+ P 2 V sinc 2 (7tmp/L)5(f + ^-y-) |P| < 1 


(C 


m*Q 

where sine x = sin x/x. 

Secondly, assuming Wgp = 2/T^ = 2R b and L » 1 then 


x(f) - s(l-|p|) 2 S(f) (PSD of x 2 (t)) (C 


Finally, we have 


ECy(MT c ) |H 1 ] = CS( 1- j P | ) 2 + N 0 W bp ]MT c ra 1 (C 

Var[y(MT c )|H 1 ] = [2S(1-|P|) 2 + N 0 W Bp ]N 0 MT c A ^ (C 

A Gaussian approximation via a central-limit theorem type argument 
(W BP MT C » 1) then yields 


*In the absence of noi se . 





P n = erfc{ 


°o erfc ^FA^ + m o " m l 


erfc'V ) - 2 ,£ ( l-| p | ) 2 /H 1 ^ 

= erfc{ — - -9— } 

/ (^ n ( 1 - l p |) 2 + 1 


(C.8a) 


0 


where 


6 T b 

n-7- 


(C .8b) 


In the above discussion, we have assumed that the second stage test 
(verification mode) use a circuit like that shown in Fig. C.l which may 
be taken from one arm of the tracking loop. On the other hand, we can 
utilize the first stage test circuit for the verification mode; in that 
case, the ROC is the same as eq. (B.10) and (B.13) derived in Appendix 
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APPENDIX D 

ON THE COMPUTATION OF P FAj (AND P^) 

D.I General Discussions 

From Appendix B, we know that the summation of N noncoherent 
integrations has the following two different distributions: 

p(y| H o> ■ ffPirr e ' y (oaa) 

and 

P(y| H l) = (|j£) N ‘ 1/2 e” (y+Ns) I N _ 1 (2/TW) (D.lb) 

where 

? E r 

B - M(l-|P|r TP . (D.lc) 

0 

and I i ( x) is the modified Bessel function of order N-l. Since all y 
correspond to code epoch in Hq region are independent (Appendix A), we 
have, for the SSP algorithm, 

Pp^ ^ Prfone of y in Hq region has the maximum value} 

= c 0 J p (y| H o){[/ Q P(x|H 1 )dx][/^ P(x|H Q )dx] 0 }dy 

" C N p (y| H 0 ) F M H i) FC ° 1 (y|H 0 )dy 

(D.2a) 


or, equivalently. 


= 1 - Pr{ the y from region has the maximum value} 


= 1 - / P(y|H,)[/ y P(x|H n )dx] C °dy 
0 0 

°° c 

=/ P(y|H,)[l - F °(y|H )]dy (D.2b) 

0 

where 

C Q ^ C-l . (D.2c; 


The above equation was shown by Lindsey [8] to have the following 
expansion: 


P FA‘ 


1 

Z 


= 9 n ” k = Q *»" * 0 


n=2 


(D.3) 


where ^ ( -y ) is the k-th generalized Laguerre polynomial of order 


N-l, and n is the k-th coefficient of x k in the expansion 


0 n 

( I 5r) n 

m=0 


nC, 


0 


I « kn * 

k=0 K ’ 


It is well known that L^"*'(-y) can be computed by [4, Ch. 22] 



7 


ir *'(-y) = y + N 


L k N1) ( _y ) = r (y +N+ 2 k - 2 ) L k-I 1 ) (-y) 


+ (N+k-2)L^2 1} (-y) 


With the initial condition: 


°0,n * 1 


, 1 < n < C. 


1/k! , 0 <; k < N-l 


a^p can be obtained via the recursion 


min(k.L-l) 

L — n — 


For the special case C n = 1, (D.2b) can be expressed by [1] 


-Ng/2 (N-l) N-l / kii.l 1 \ 1 


' k=i (k-i ) !2 


Unfortunately, both (D.3) and (D.6) become useless when N or is 
very large. For large N and C=2, Marcum [1] derives the following Gram 
Chari ier series expansion for Pp A : 


P fa (2) = \ Cl- < f,' 1 (T)] + C 3 «t, (2) (T) - C 4 4> (3) (T) - C 6 4- (5) (T)... (D.7) 


§. 


o 


2/In (i+e) 


• _ 1+20 
A O * 

* 8N(l+0)^ 


6 i6N(i+g) 3 ’ 


<t>(y) = 


J_ e -y 2 /2 

m 


> 



(y) denotes the i-th derivative w.r.t. y and 


r X (y) - — - \ e“ d /2 da 
/ZT -y 


ing the union bound along with (0. 


(0.8c) 

(D.8d) 

(D.8e) 


(D.8f) 


P FA < C 0 P FA(2) 


(0.9) 


In case of large Cq» more terms on tlie Gram-Charl ier series may be 
needed in order that (D.9) can render useful results. 

D.2 Numerical Algorithm 

We now proceed to describe a direct numerical evaluation algorithm 
for Pp A by use of (D.2b). 

C 

First let us define F(y) - 1 - F °(y|H 0 ). Since F(y) is a monotone 
decreasing function of y for any fix N and Cg» ^ ’ s eas ^ t0 see 


(D.10) 


P FA » / P(y|H 1 )F(y)dy 

■ J P(y|H,)F(y)dy + R(o) 

0 1 

where 

00 

R(«) * / P(y | Hj )F(y )dy 
a 

00 

< F(a) / P(y|H 1 )dy 

a 

< F(a) . (D.ll ) 

Pp A can thus be obtained by evaluating the 1-st term on the right 
hand side of (D.10) with a truncation error less than F(o). 

To have an accurate calculation of the integral 

/“ P(y |Hi)F(y)dy (D.12) 

one has to have algorithms to render precise values of P(y|H^) 
and F(y) in the first place. P(y|H^) can be obtained from standard routine 
[3] with minor modifications to avoid overfloating or underfloating 
problems. F(y) will be evaluated by the following algorithm* suggested 
by Chie [2]. 

Let Kq be the integer such that 


*Simple closed form approximations valid in some region can be found, 
fot example in [4], or [7]. 



then 


r~T = max 
V UUH-1 


k n = [x] - decimal part 



of x 


ORIGINAL PALi£ SS 
OF POOR QUALITY 

x > 0 
if N-l > x 


= N-l 


if N-l < x . 


Let [10~ z ,10 z ] be the range of real numbers the available computer can 
handle, then 

F(y|H 0 ) = 1 - e^yje^ (D.13) 


where 


Hence if 


then 


N-l k 


Vl (y)e ' y * kT )e ' y 


'0 


< (N-l) f-r e" y < 1 


. k o 


(N-l) jh- e' y < 10” z 
O’ 


F(y|H 0 ) = l . 


(D.14) 


Assuming (N-l)y ^e’^/Kg! > 10~ z henceforth, we shall describe the 
algorithm for two distinct cases. 


Case I. N-l < x 






2) Compute 


until 


k ORIGINAL FAClT 'S 

S, = [ t. OF POOR QUALITY 

* i=l 1 

(N-l-k)t k+1 /S k < 10" E or K = N-l, 


wher" E is a predetermined positive integer. 

3) Let K be the largest integer which statisfies the above 
condition and set F(y|Hg) = 1-S£. The resulting truncation error will 
be less than 10“^[l-F(y | Hq) ] . 


Case II. N-l > x 


Let 


2) Compute 


»i 4 fr *' y 

1 K q . 

V 1 

- i^iTyr e ' y 


s k - l S 

K i = l 1 


unti 1 

(K 0 - K )t K+ i /S K < 10 
or 


3) Let K be the largest integer which satisfies the above condition 


and let 



4) Compute 


l £-K 

S * = * S K + £ l K*1 


until 


or 


[N-l-£-(K 0 -K)]t £+1 /S A < 10 


-E, 


N-1-Kq < (£- K) 

The resulting truncation error, after setting F(ylHg) = l-S", where a 

X» 

is defined similar to K, will be less than 


-E -E 

(10 *+10 2 ) S £ . 


Observing that F(y) behaves Mke a step function (see Fig. 0.1), we 
further split (D.12) into two integrals: 

«i 

/ P(y | )F(y )dy + / P(y|H 1 )F(y)dy (D.15) 

0 ctj 

where the smaller one between Ffag) = 0.9 and 02 = N(l+e). Each 

integral will ir. turn be calculated according to an adaptive Gaussian 
quadratic rule [2]. 

D.3 Gaussian Approximation 

Note that both P(y|Hg) and P(y|Hj) converge to a Gaussian density 
as N ♦ « (see also Fig. D.l). This fact suggests that we substitute the 
following approximations 




Figure D.l. Behavior of F(y) and PCylH.). 






with 


P(y| H 0 ) = - _i_ — exp - (y-m 0 ) /2o 

/2ttoq 


p (yj H 1 ) * -zz - exp - (y-m 1 ) 2 /2oj 

/2 noj 

= N( 1+0 ) 

K 

N( 1+23 ) 

N , 



into (0.2b) and change the low limits to -» to get 

“2 “2 
(u+m) (u-m) 


p = 1 - * 

FA /27a 

00 

I [e 
0 

CM 
< O 
CM 

1 

c 0 

$ u (-u) + e 

where 

* A 
a = 

°i /a o s 

(l+2e) 1/2 


; ft 

V m o 

°o 

= /NB 


4>(u) 

A 1 U 

= ^=J 

#-**'*« . 


/ar -°° 


For Cn - 1, (0.17) becomes 





The Gaussian approximations (0.17) and (0.19) are compared with direct 
computation described before. Some selected numerical results are shown 
in Table D.la-D.lc. 


Finally, if we substitute the power series 


'1 


<fr(u ) * + 6 


-u fe /2 - 


/2iT n=0 


I P(n)u 


2n+l 


P(n) d [1.3*5 ... (2n+l) ] -1 


into (3.19), we find, a f ter some manipulations. 


i i ~2 2“2 

, _ 1 1 r f m C 0 > l 

FA " 7 'a>; P[ ' ( & 2 ‘ ai ’ 1 


where 


n=0 un /2U 


m/a 2 


f. I ] 

■(2n+l)'- y 


F(n) = 4n!/C 


n+1 


and 


r _ 1 + J_ 

7 2 ; 2 

2 t 2 

. p -y / 4 * -yt- •*— i 
°-n<>) ‘ TtST < 0 ' ‘ dt 


(D.20) 

( 0 . 21 ) 

( 0 . 22 ) 

(0.23) 

(0.24) 


is the Parabolic Cylinder Function. Unfortunately, this series 
expansion cannot tell us the behavior of the P FA explicitly. 



Table D.l. Gaussian Approximation, Pc A (G), for False Alarm Probability 
P fa (E). (a) C Nq = 10' * = 512. (b) C Nq = 10 3 , B = -6 dB. 

(c) N = 256, b = -6 oo 


Table D.la. 


B (dB) 

-6 

-5 

-4 


p fa( g > 

2.82(-2) 

1.69C-3) 

1.95(-5) 


p fa (e) 

3.50(-2) 

1.98(-3) 

1.75(-5) 



Table 

D.lb. 



N 

32 64 

128 

256 

512 

p fa (g) 

.924 .833 

.62/ 

.270 

.028 

p fa (e) 

.958 .891 

.703 

.329 

.035 


Table D.lc. 
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APPENDIX E 


ON THE PERFORMANCE OF S^A.B) AND S 2 (a,b,Bg,N t ) 

E.l Performance of Sj(A,B) 

Let {X ,n=l,2,3, . . .} be a sequence of i.i.d. R.V.'s. Define 

A k 

S k = £ X n , where Xg = Xg with probability one, and t 

A 0 = 0 

- inf{k: S k tf(A,B)}, where A > 0 >_ B. If f ( x | ) is the common 
probability density function of X n , n > 0 given that the hypothesis H. 
holds, then 


W ~ Pr £ X o = x 0’ S t ed i -1 Xge(A.B) 

A 

= Pr[Xg=x () ,S 1 ed i ] +/ P^ (z)f(z-Xg|H i )d 2 

B 

A A 

~ " i (Xg )+/ P i (z)f(z-x 0 |H. )dr (E.l) 

B 


where 



x > A} 

( E .2a 

x < B} 

(E.2b 


Similarly, if we define 


W = E l t l x o =x O* S t ed i^ 

(E.3) 

A 

E i ( x g)= 1 +/ E^ (z)f(z-x Q )dz 

(E.4) 


then 



UZm\ 


+i 


'tjCinC^c 


om 


Eq. (E.l) and (E.4)* were given by Anscombe [1] and Kemperman [2]. Two 
similar equations were obtained by Albert [3] for analyzing the so- 
called generalized sequential test. Kendall [4] then applied Albert's 
work to the special case where 


X i = aR i - b , 


(E.5) 


f(*|H 0 ) = i e -(x+b/a) 


x > -b 


(E.6a) 


f (x|H 1 ) = I exp-(^ +8)Ig(z/ f (x+b)) x > -b (E.6b) 


The first set of integral equations, i.e.. 


P n (x n ) = exp - j (A-x 0 +b) + / 


max(x 0 -b,B) 


_i -(y+b-x n /a) 

a dy 


B < Xg < A 


(E.7) 


. A -(y-x n +b/a) 

W * 1 * T / , „ E o<»> e d » 

max(Xg-b,B) 


B < Xg < A 


(E.8) 


were solved by Kendall [4], while the set corresponding to the Hj case 
has eluded all effort to derive a closed form solution. Nevertheless, 


‘Both belong to the class of Fredholm equations of the 2nd kind. 


oCinC^t 


om 


for the case we are interested in, i.e., Pfa 4 P 0 (°) << P 0 ~ p l(°)> or » 
A >> | 8 | , the so-called Wald's approximations are valid because of the 
fact that B < 0 < E(Xj) « A. Now, let us look at the solutions of 
(E.7) and (E.8) first. They are 

P FA P 0<°> 


G[-DB;Db] .-aD(A+b) 
‘Df A-B+b) :Dbl e 


(E.9) 


T 0 = E Q (Q) 


= e" aDB H[-DB;0b]+P n (0){l-e aD(A ' B+b) H[D(A-B+b);Db]}, B<0<A (E.10) 


where 


5(x;c) ^ 1 + l c < nc < x < (n+l)c ( E . 1 1 ) 

j=l J * 


-ax _ r ir 1 (jc-x) 1 


:;c) - (n+l)e - I I - c < nc < x < (n+l)c 

j=l i=0 ila 3 * 1 


( E . 12 ) 


a » e b/a . a -1 D 


(E.13) 


The above formulae become useless when (A-B)/b » 1. Good 
approximations utilizing Laplace transformation method were also derived 
by Kendall [Appendix, 4]. After some rearrangements of Kendall's 
approximation formula, we reach the following numerically stable 


equations: 


p = (aDb-l)e D8(T ~ a) -(l-YDb) 

FA (aDb-l)e D(h+b)(a " Y) -(l-yDb) 


(E.14) 


T Q = (aDb-1 ) 1 [ 2 (“nb-l) ' 1 * aDB l ’ (a-Y)(l-YDb) 


-DB(a-Y) 


p Ji r aD(h+b)+(ct0b) 2 /2(aDb-l)-l l 

K FA tx " 1 oDb-1 " (1-yOb) ( o-y) Jf 


(E.15) 


where a,y are two real roots of S = e SOb , y < e < a, and h ^ A-B. 

Eqs. (E.14) and (E.15) have a truncation error less than 

10-(3+[A-B/b]) f [ x ] being the integer part of x. Hence (E.9) and (E.10) 
A-B 

will be used if [-jj-] < N, where N! is the largest real number which 
computer can handle, otherwise use (E.14) and (E.15). 

On the other hand, Pg - Pj(0) can be obtained through the two-step 
algorithm: 

(1) Solve the transcendental equation 


E[e _X “] * 1 


(E.16) 


(2) Substitute the nonzero root* wj of (E.16) into 


W-iB 

p . ..L-e 

D u>. A u.B 
e 1 - e 1 


(E.17) 


Finally, Tj - Ej(0) is to be computed- by the use of Wald identity: 


*the existence and uniqueness of such a uj was proved by Wald [5]. 



* * * 
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E^O) = E(X t )/E(X 1 ) 


(E.18) 


E.2 Performance of S^a.b.Bn.M 


k . 

t = I R u 
i=1 u 


(E.19) 


R u - aR u- b 


lf S k-1 > B 0 


otherwise 


(E.20) 


where a > 0, b > 0, Bq < 0 and is defined in (8.1). 
Test S 2 (a,b,BQ,N t ) is defined by the decision rule: 


D 0 : lf S k,t 4 B 0 

D l,m : lf S k,m > B 0 


V 1 < z < M 

V 1 < k < N t , 


and = maxfS^} 


0: otherwise 


where 


(D.21) 


Dg - accept Hg ( see Appendix B), stop testing the present section 
and go to the next section. 

D l,m " re J ect H Q , 9° t0 second stage test. 

U ^ continue by taking another set of samples. 


In other words, S 2 (a,b,Bg,N t ) is a modified version of the so-called 
variable dwell-time ( VOT) system analyzed in [6], the modification being 
that Dg is not made until all ^ have crossed the boundary Bg and Dj j(T] 
is made at N t -th noncoherent integrations by choosing the maximum one. 


o 


-M 


‘cJCinC^c 
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Conceivably, the improvement made by this test will not be so 
impressive as that by the VDS over the fix dwell -time system. This is 
evident from Table E.l shown below. Other sequential procedures using 
CCD PNMF's are considered to be inferior to S 2 since each trajectroy 
S k ^ is independent of one another, no conclusion as to the 
categorization (H^ or Hq) of one trajectory can be drawn from the 
behavior of a group of different trajectories (or a single one). 

Table E.l. 


Mean 

Dwell -Time 




Detection 

False-Alarm 



Probability 

Probability 

Fix 

M=1 

p D 

P FA 

Dwel 1- 


Time 

M=m 

P D 

m p FA 

(FDT) 


VDT 

M=1 

P D 

P FA 

S 2 

M=m 

P D 

m p FA 


\ ' T 0 

The mean dwell time can be computed by modifying Wald's 
approximation as described in [6], or, if an exact answer is preferable, 
use the recursive algorithm in [7]. 

Let n be the random stopping time at which Dq is reached and 
F ( n | ) be the distribution function of n under H^, given that M being 
equal to one. After obtaining F(n|H^) from [6] or [7], we set* 


P D = Cl^NjHj)] 

P FA = (M-l)[l-F(N t |H 0 )3 


(E.22) 

(E.23) 


♦Assuming Pp A << Pq, 


aCinCom 
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* 
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APPENDIX F 

DERIVATIONS OF THE MEAN ACQUISITION TIME (T^) 

AND THE DISTRIBUTION FUNCTION OF THE ACQUISITION TIME 


The mean acquisition time for the SPS algorithm is 
straightforward. After the initialization period which takes MT C sec 
the acquisition process behaves like that of a two-dwell system with 
enlarged cells ^bins). Hence 


T ACQ = T d + (p Q “ 2^ T D 1 +P FA 1 T D 2 +P FA’ T p T c )N d N t 


(F-l) 


where 


T D = MT c* 

x D = Nx d = dwell time for the first stage test, 
t 0 ^ = dwell time for the 2nd stage test, 

T p T c = penalty time, 

p d = P D 1 * P D 2 * 

P FA = P FA 1 * p FA 2t 

N d • N t = total number of bins in the uncertainty region 

= (number of Doppler bins) • (number of time bins). 


Eq. (F-l) is actually an upper bound only since 1 + Cq is not a ^ wa ^ s 
equal to an integer multiple of M. In other words, one of the N t time 
bins must have less than M cells. Thus tq^N^N^ can be reduced to 
(l+C 0 )NT c« 

\n exact analysis expression for F(T^qq) cannot be found. However, 
if P FA^ << p FA^ t D 2 <<: T D|» then distribution function for single- 


csCinCoi 


^ a» W ■*- 





dwell system derived by Di Carlo and ,'eber can be used as an upper bound. 


For the SSP case. 


oo n 


Tflrn = T n + 1 1 CnT A +( n-k ) T ]Pr{correct detection at the i.-th try, 

u n =i k=l A p 

(n-k) false alarms} 


■ l .1 CnT A+ ( n-k )T D ( 1-P D -P FA ) 1 P E 


n=l k=l 


where 


t a = T D * (C N 0 +1)NT c * T D 2 • 


(F-2) 


T c ’ T p T c • 


After simplification, we obtain 


T * _ u _ _v_ T - T- 

(MO 7 (l-v) : 


y _ p _ l r A r u v i p u 'p u ■ 

t acq “ p d i ttv i t rrr? " tttt? j ~ thj ' » 


where 

u = P FA , 

v = i - P D - P FA . 


If P FA << Pq, then 


' ACQ 


V 

F" T [ 
0 1 


(F-3) 


which is eq. (4-5) in the main text. 
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APPENDIX G 

AN ADAPTIVE PN ACQUISITION SYSTEM 

The advantages of usiny CCD PN matched filters (PNMS) have been 
demonstrated in yreat detail [1]. It was shown that the presence of 
data modulation put a limit on the coherent integration period, M. 
According to Grieco [SJ the expected data transition loss (DTL), L, is 
given by 


I-M/4G p 

L = 

l-3M/4Gp 


(NRZ encoding) 
(Manchester encoding) 


(1) 


where 


code chip rate > M 
data rate 


Thus the TDRSS SSA service, having a large data rate operation range 
(0.1-300 i^b/sec), calls for an adaptive coherent integration length. 

Using the SPS algorithm, PN code can be acquired within 2.6 seconds 
with a probability of .92 [1J. The design system parameters for the 
above data are as follows. 


N d (Doppler bin) 

3 

Cy+1 (cells per Doppler bin) 

76UU 

6T c (step size) 

U.5 

Sampling Loss (dB ) 

2.S 

Doppler Loss (dB) 

0.4 
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Data Transition Loss (dB) 
Code Rate 
M (chips) 


U.2 

3 Mchi ps/sec 
512 


From (1), a 0.2 dB data transition loss is resulted from a data rate of 
1.05 Kb/sec. 

The RCA TC 1235A CCD correlator designed by Munroe of Lincoln 
Laboratory [2J has a programmable length of 64, 128 and 256 (chips). 

Two such correlators in tandem will allow a programmable length of 64n, 

1 <_ n <_ 8. For a code rate of 3 Mchips/sec, we propose that the break 
points be at data rate = 3 x 10^/64n Kb/sec. In other words, the 
correspondence between data rates and coherent correlation lengths is as 
shown in the following table. 


M 

512 

448 

384 

320 

256 

192 

128 

64 

(chips) 









Data 

0.1 

5.9 

6.7 

7.8 

9.4 

11.7 

15.7 

23.4 

Rate* 

5.9 

6.7 

7.8 

9.4 

11.7 

15.7 

23.4 

50.0 


(Kb/sec) 


*one channel 


These break points allow at most one possible data transition during a 
coherent integration period and the resulting system performance will 
not be worse than the design point, i.e., 90% acquisition time £2. 6 
seconds. To show the latter claim, let us consider the worst case when 
the data transition occurs at the middle of a coherent integration 
period. Since tne statistics under Hy is not altered by the presence of 
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data modulation we shall consider only the detection probability. For 
the above worst case, there is at most a 3 dB loss if we assume equally 
1 . k.e ly ±1 data modulation. Because then, there is one data transition 
for every two periods and whenever there is a data transition at the 
middle of that period the test statistic produced at the end of that 
period will be no worse than zero which accounts for the 3 dB loss. 
Hence, for example, when M switches from 512 to 448 at 5.9 Kb/sec a 3 dB 
worst data transition loss plus a 0.58 dB processing gain loss can be 
well compensated by a 7.5 dB C/Nq increase due to the data rate increase 
(from 1.05 K to 5.9K). With the same M, the OTL increase due to data 
rate increase can also be made up for by C/Ny incrase; see eq. (1). 

Other cases can be checked similarly. 

For data rate greater than 50 Kb/sec, we have to switch to 4- 
channel sequential detector which can be easily configured with devices 
borrowed from code tracking loop (a 2-channel double dither loop). 

Based on the numerical results presented in [3], we find, after 
necessary conversions, that the mean acquisition at 50 Kb/sec is less 
than 1.2 sec for the 4-channel sequential detector. The 90% acquisition 
time is expected to be less than 2.6 sec then. Two facts support the 
above estimation: (1) distribution of acquisition time for serial 

search algorithms can be well -approximated by that of a Gaussian random 
variable if Cq » 1 [4], (2) mean and variance of the acquisition time 
can be minimized by the same set of parameters [3]. 
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APPENDIX B 

PN CODE TRACKING LOOP 


B.l Introduction 

Code tracking is the second step of a two-step code synchronization 
procedure. It will be initialized right after code acquisition or 
coarse synchronization has been achieved to further reduce the code 
epoch offset between the received code and the incoming code and 
maintain this state for as long as possible. 

Several different types of code tracking loops have been 
investigated in the literature, and two of them, namely, the delay- 
locked loop (DLL) [1-2] and the tau-dither loop (TDL) [3], are widely 
used and thoroughly analyzed. Both of these configurations can be 
operated in a coherent or noncoherent mode. For our application, PN 
acquisition and tracking must be performed prior to carrier 
synchronization, and hence noncoherent tracking is of interest. 

One disadvantage of the DLL is the difficulty to maintain a close 
match between its two envelope detectors. Gain imbalance problem which 
results in tracking error can be eliminated by applying the TDL. On the 
other hand, the TDL suffers an approximate 3 dB loss in noise 
performance [3] due to the time-sharing of the early and late gates. 
Recently, Hopkins [4] proposed a new phase detector configuration called 
double-dither loop (DLL) which combines the desirable features of the 
DLL and TDL, i.e., the 3 dB dithering loss is recovered with no gain 
imbalance effects. By modifying Hopkins’ DLL the AIRS’code tracking 
subsystem employs a 2-channel DLL (2-DDL) in order to: 


• eliminate gain imbalance problem of the DLL 
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• pick up 3 dB dithering loss of the TDL 

• utilize the weak (I) channel signal energy not used by the one- 
channel DDL. 

B.2 Analytical Model and Timing Diagram for the 2-PDL 

The 2-DDL is illustrated in Fig. B.l. I E (I L ) and Q E (Q L ) denote 
the advanced (delayed) versions of the local I and Q channel PN codes 
respectively. The timing diagram of the dithering signals, to Sg is 
illustrated in Fig. B.2 where 6j and Gq denote the I channel and Q 
channel gains. Notice that at any given time if the input signal r(t) 
is cross-correlated with advanced (early) and delayed (late) versions of 
the local I channel PN code in the upper channel then in the lower 
channel it is cross-correlated with those of the Q channel PN code and 
vice /ersa. The dithering gain signals, S 4 , Sg ar. sign signals, S 3 , S 7 
are designed in such a way that proper gain and sign are assigned to the 
outputs of envelope detectors and a proper control signal is formed at 
any time. Related system parameters and definitions are listed in Table 


B.l. 

Following are assumptions made in the subsequent analysis of the 2- 

DDl: 

(1) [i| (loop bandwidth) « 1/T C (chip rate), i.e., code self noise 
can be neglected 

(2) B l *< 1/Tjj (bit rate) i.e., modulation self noise [ j can be 
neglected 

(3) 1/T d (dithering rate) « 2/T b (» bandwidth of the BPF). 

(4) Long PN code 

(5) Zero code Doppler 

(6) The time offset between early and late codes is one chip time. 
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Table B.l. System Parameters and Definitions. 


CODE RATE 
DATA RATE 

INPUT SIGNAL-TO-NOISE RATIO 
EARLY-LATE GATING 

LOOP BANDWIDTH 
Q: I POWER RATIO 
DITHERING RATE 
CORRELATION LOSS 
CODE PHASE JITTER 
MEAN SLIP TIME 
NORMALIZED MEAN SLIP TIME 


1/T b <=R b > 

E b /N 0 - ST b /N 0 

A = T c /2 
(Half a Chip) 

G 

1/T d 

c L 


log 10 (2B L T L > 
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The 2-DDL of Fig. B.l can be modeled as that shown in Figs. B.3 and 
B.4, given the above assumptions and following a similar argument of 
Hopkins [4], The corresponding dither signals which select appropriate 
outputs sequentially are shown in Fig. B.5. 

B.3 Code Phase Jitter Analysis 

In Fig. B.4, is the two-sided bandwidth of the bandpass filter 
of the i th bandpass detector BPD^) and B^ is that of the low pass filter 
of B p Dj. A correlator through which the received signal is cross- 
correlated with PN sequence and four parallel 

G-again BPD's in cascade is called a Ag-G detector (Fig. B.4). Let us 
denote the signal and noise components of the output of BPD^ in a Ag-G 
detector as G*y^ , and G*n^ ^ respectively. Then it can be shown [4] 
that 


\,i * k i n i (lt ‘ , A 0 ,i ) * v i 


E( V> = ° 


(B.l) 

(B.2) 


2 

°A i 
m B*' 


a Var( V> 

2 2 ^®i 

■ "i-hCttN 1 * 2 * „.i> 


(B.3) 


where 

all j are independent of each other 
k i = gain constant of BPD^ 

= dc offset of BPD.j 

n i = N 0 w i 
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<±r 


■v 

1 


0 


r -4 " - ' ■ ■ 'IJT' V <r "" 
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(B.4) 


OF PGO'< . J 

R ^ T + 2^^ ' f * *E 0r 

■ *B- I L or «L 

s is the total signal power 
T i * l/«l 

and 

T c T 

\[rt = E[PN(t)PN(t+ T ± ^-)] . (B.5) 

The output of the 2-DLL phase detector, omitting time and delay 
arguments, is 

e D = d lt G I Cy I L ,r y I E ,2 ] + 

+ d 2{ G I^I L ,2" y I E ,l ] * G Q^ y Q L ,4” y Q E ,3^ 

+ d 3 { G l ty l L ,3 _y I E ,4 1 + G Q^ y Q L ,l” y Q E ,2^ ^ 

+ d 4 ( G I^I L ,4“ y I E ,3 ] + G Q^ y Q L ,2” y Q E ,l^ ^ + n(t * 

(B.6) 

where n(t) is the totality of all noise terms. The phase discriminator 
characteristic (PDC) is defined as the expectation of e D (t,-r), which is 
denoted D(t), i .e., 

r 'c) = E(e 0 ) . (B.7) 

B> the independence of {n. }. ».q. (B.l), (B.2), (B.4), and the 



equation E(d^) = 1/4, we obtain 


0(t) = f (k 1 +k 2 +k 3 +k 4 ) (Gj+GqJQrI ( t - 2^Hc( t + ^J] • ( B - 8 ) 

average of the phase detector output is in proportion to the mean of all 
four gain constants. 

The normalized PDC or normalized S-curve g(x) is obtained 

T ) = O^fo T 

- (B.9) 

which is shown in Fig. B.6. 

The variance of e D can be derived from (B.l) to (B.l) and the 
identities: 


E(d i ) 

= 1/4 

E(df) 

» 1/4 

C(n 

m B 

i ,n C 1 ' 

»' '-p* • 

2 

V 

_ 2 

2 

V 

c vj t p J 

II 



It is found tnat 



Var(e D ) 


(B.10) 


( G I + 6 q) 


N 

l 

i=l 


i + 0 i i)] - 
» 1 * 


3E 2 (e D ) 


To assess the performance of the 2-DDL and to compare with other 
tracking loops, the following balanced conditions are assumed: 


B^ = B , for all i 
t = 0 


2 2 

Under these conditions E(e n ) = 0, all <£ • and o, - are equal, 
T T r V 1 V 1 

R (t - ~) = R c (t + ^ = 1/2, and (B.10) becomes 


*»r(e D ) - 4k 2 N 2 BH(l ♦ ^)(G 2 *G 2 Q ) 


(B.J1) 


Normalizing Var(e 0 ) with respect to D'(0) as in the case of E(e D ) yields 


2 A Var(e D } 
°D 


to* (o)] fc 

(G j+Gq) , , 

“ ~2 SNRq 1 + SNRj ] 


/ n 


( 8 . 12 ) 




r. 


c 

t 
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Substituting G = Gq/Gj into (B.12) gives 


2 

<\> 


2 (SNR 0 } t 2 + SNR, 1 


(1+G) 


(B.13) 


I 


It was shown [4] that the code phase jitter can be determined by 
replacing 2B with the equivalent noise bandwidth of the loop B^, if B L < 
2B. Hence the normalized phase jitter is 


°2-DDL = < SNR L> 2 + (SNI V *1 


(1+G) 


(B.14) 


where 


SNR 


2S 


L N q B l 


and 


SNR 


I n q w 


(B .15a ) 
(B .15b ) 


For comparison, we list the normalized phase jitters for TDL and 
DDL below, assuming only the stronger channel signal is tracked. 


<^ QL = (G.SNR l ) -1 [.905 + (.453 - (G-SNRj)" 1 ] (B.I6) 


o* DL = (G .SNR l ) _1 L.5+(G -SNRj )3 _1 . 


(B.17) 


When G = 1, it is noticed that 


_ 1 2 


5 2-DDL ' 2 °DDL 


I 2 

4 °TDL * 
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For high input SNR, the improvement for the 2-DDL over TOL is about 3 dB 
in terms of required SNR. Fig. B.7 and Fig. B.8 show the normalized 
phase jitters for three different loops in low SNR and higher SNR 
environments respectively. 

B.4 Mean Time to Lose Lock 

Given an initial code offset tg, the time to lose lock T L for a 
code tracking loop is defined as 


T = inf {t: D(T)=U,T=x n at t=0} (B.lBa) 

L tiO u 

The mean time to lose lock T L is then given by 


T l = E(T L | 1 =-^ at t=0) . (B .18b ) 


Evaluation of T L requires nonlinear analysis of the tracking 
system, and, for a first order loop, a general expression was obtained 
[5] in terms of a double integral for arbitrary phase detector 
configurations. In case tq=0 and the phase detector is characterized by 
(B.9), the double integral can be simplified into a product of two 
single integrals [6]. Since | tq| < .25 and in general tq is close to 
zero for our system, it is appropriate for use to use the result of [6] 
as an approximation. The expression for 7^ in this case is 


4a 2 (2B. 


7 U 
) 0 


3/2 


-[G( e)/ a ] 


3/2 


de}{J 


[G( e' )/a ] 


de' } (B.19) 


where 

G(e) = / g(t)d T , 


inK^om 
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g( e) is the normalized phase discriminator characteristic (PDC) shown in 
Fig. B. 6 and a- is the code phase jitter. The value 3/2 is derived 
from the fact that the late and early version of the local PN code is 
offset by one chip time. 

For an arbitrary initial condition tq, it can be shown that 

, 3/2 r( w 2 1^1 r , , w 2 

- T.l +-^f e ‘ G(e)/o f e G(e )/a d G 'd e (B.20) 

L T n =0 2 a“ -3/2 0 


The normalized mean time to lose lock, log^g (2BjT^), as a function of 
phase jitter a is dotted in Fig. B.9 and B.10 for selected values of a. 
B.5 Pull-In Timt and Lock Detector 
B.5.1 Pull-In Time 

The noiseless pull-in behaviors for various types of PDC have been 
studied by Spilker [1], Gill [2] and Nielson [7]. In the presence of 
noise, the pull-in time defined by 


T„ = inf It: I t| < t: T=x n at t=0} (B .21 ) 

K t>0 p u 


is a random variable, and so it is only meaningful to look at the 
expectation of Tp: 

T p = E(T p | at t=0) . (B.22) 

However, such an expectation does not exist because there is a nonzero 
probability that T p = «. The reason is as follows. 

The noise performance of the normalized code phase error 
e(t) - is governed by the differential equation: 
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d e ■ -4B L g(e)dt + dW Q 


(B.23) 


O 

where f dW Q is a standard Wiener process and a is the phase jitter. In 
deriving (B.23), we have made the assumption that the noise component 
out of the phase detector, n(t) of (B.6), is white Gaussian noise. 

For small e, g(e) *• e. Substituting g(e) by e in (B.23) yields 


de s -4B^ e dt + /3B^o dWQ 


(B.24) 


which shows that (e(t) } is an Ornstein-Uhlenbeck (0-U) process The 0-U 
process of (B.24) has a transition density function [8J 


where 


2 

°0 


-2oQt "Ogt 


P(t,x,y) = <>( 2 ^ (i-e ),xe ,y) 

♦(t.x.y) . e -(y-*) 2 /2t 

✓2 Tit 


(B.25a) 

(B.25b) 


^ = 4B L 


a Q = (8 B l ) 1/2 o . 


Eq. (B.25a) tells us that there is a nonzero probability from e = x = T g 
to y = 3/2 even with a PDC g(e) = t. However we have to be reminded of 
the fact that at e = 3/2, g( e) = 0, and so there is a nonzero 
probability that e will remain to be larger than 3/2 forever. 

In order to assess the noise pull-in performance with a suitable 
criterion we neglect the probability that | e| > 1/2 and apply the 



. r* » 
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approximation (B .24) for je(t), 0<t<T}, where T is a reasonably large 
number, say 10 seconds. That is, is now defined by 


/ tP(t,x,y)dt . 
0 


(B.26) 


B.5.2 Lock Detector 


The lock detector structure of the code tracking subsystem is shown 
in Fig. B.ll which is the same as that of a post-detection integrator 
used in the PN acquisition subsystem. The operation characteristic of 
such a detector has beer derived in Sec. 2, which we now rewrite in eq. 
(8.27) 


erfc { 


/ g p! ( 1 -M > 2 + 1 


} (B.2/> 


where 


FA = design false alarm probability 


E c 

^ ' N 0 T b 


E b T c 


V N 0 


G p = T b /T c £ processing gain 
M = integration time (chips) . 


)■ 

i 


Eq. (B.27) i? plotted in Fig. B.12 for P pA = 10" 8 , E C /N Q = -32 dB, l/\ 
5 3 M chips/sec and t = 1/4. Note P Q = .99 can be achieved within .085 
sec with the above specified worst case. 
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3.6 Band! imiting Effects of Pre-Despreading Filter and Bandpass 
Arm Filters 

B.6.1 Correlation Loss Due to RF Filtering 

The model under consideration is illustrated in Fig. B.13 where 


r(t) = /2S PN(t-r)d(t) cos( + n(t) 


= received waveform 


H( b>) = the transfer function of the RF bandpass filter which is 
centered at u^. 
r(t) = /2S PN(t--c)cos(a^t+(fi) 

= local estimate 

r(t) = the filtered version of r(t). 

Let the truncation function [7| (t) be defined by 

1 |t| <T 

0 otherwise 


fTj(t) = { 


and the cross correlation function between r(t)[7j(t) and r(t)[7|(t) by 


R(6,*) = lim|y/ E[r(t)r(t) fTj(t)]dt 
T ■*■<*> -T 


whc'e 
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Then it can be shown [10] that 


jc TB o 

R(0,0) = H / T c sinc^UT^Jda, 

* -nB 


H( uru^) = 


-1 mo 


, J aij < irB 


, otherwise 


In case the HU-^ n )has a nonlinear phase #(w) instead of constant phase 


/ST wB n o 

R(0,0) » — ^ f sine i uJ c /2) cos Sufdu, 

- wB 


where 


^ _ 1 d u>) 

^ do? orO 


The number R t (0,0) corresponds to a signal energy degradation in the 
best case where 5=^=0 and hence is called correlation loss. This 
correlation loss, as a function of bandwidth-time product is depicted in 
Fig. B.14 for the case a = 0. For further details see [9,10]. 

B.5.2 Bandlimitinq Effect in Tracking Loops 


The bandlimiting effect due to b*n Ipass arm filters in the tracking 
loop ir most easily accounted for by a signal energy degradation 
O' defined as 


oCinC^c 


^ otn 






JjnCi 


om 


where 


ff = E[(an)r] 

<T(t) = filtered version of d(t). 




It was shown [6] that ff can be expressed in the following way: 


O' = 


/ V f >|H (f)| 2 df 

— CD * 

/" s d ( f > df 


Hp(f ) = the transfer function of the BPF 
S d (f) = the spectral density of d(t). 


Simon [11] has examined the band limiting effects of the BPF for 
both DLL and TDL. For TDL, he found the optimal B,/R b , in the sense of 
minimizing jitter, lie between 2.5 and 3.0 for 0 <_ E b /N d 14 (dB) and 

B i T d = 4 * 

As for the choice of T d , Hartmann [3] suggested that B^/4 < 1/T d < 
B^/2, since, as mentioned before, the dithering rate 1/T d must be chosen 
small enough for the dithered signal to be passed by the IF BPF and, on 
the other hand, be large enough such that the sinusoidal component at 
l/2T d be out of the mainlobe of the loop filter. 
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APPENDIX C 

FREQUENCY ACQUISITION 


C.l Introduction 

The implementation and performance of the frequency acquisition 
subsystem are described in this Appendix. Supporting software and its 
set up can be found in the next two appendices. 

A carrier tracking loop does not acquire when the initial frequency 
error is comparable to the loop bandwidth. It is thus necessary to 
align the received and local carrier frei encies to within the loop 
bandwidth such that IF/or baseband process can take place. To avoid 
false lock in conventional frequency searching scheme and to make a 
rapid pull-in from a large initial frequency error possible we adopt a 
digital automatic frequency control (AFC) loop modified from Cahn's 
composite AFC/Costas loop [1]. The remainder of this appendix includes: 

• Analytic model and description of the AFC loop 

• Discriminator characteristics and pull-in range 

• Simulation model 

• Performance analysis 

• Performance results 

• Lock detector and its performance 

t Modification for DG-1 

C.2 Analytic Model an d ~>cription of the AFC Loop 

The AFC loop or frequency acquisition loop (FAL) for the AI..S has 
the configuration of Fig. C.l. The frequency synthesizer output is 
assumed to be of the form 
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The combination of mixer and low pass filter (LPF) serves as a frequency 
difference detector (FDD) which results in baseband outputs of 


I ( t ) = d j ( t ) cos[A<dt+$(t)] - /Gdq(t) sin[ Aait+(|>(t)] + n j ( t ) (C.la) 


Q ( t ) = dj(t)sin[Aa>t+$(t)] + v'Gdq(t )cos[ Aoit+ «j>(t ) 3 + nq(t), (C.lb) 

where dj(t), dq(t) are the ti/2S' valued data stream, assuming filtered 
by the FDD without distortion, of I and Q channels respecti vely. S' is 
the I channel signal power, Au is the frequency difference between the 
received and local carrier, G is the Q: I power ratio, and nj(t),nq(t) 
are independent narrow band Gaussian process with one-sided power 
spectral density Ng. 

The sampling and summation combination is the digital realization 
of an integrate-and-dump (I/D) filter. The relationship among the 
sampling rate (1/T $ ), the bit rate (1/T b ), the number of samples for 
each update (N) and the number of updates per bit time (M) can be 
summarized in the equation 



The control signal, e^, which is proportional to frequency error is 
derived by 
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e k = X k-l Q k " ! k^k-l 


(C.3) 


where I k ,Q k are the outputs of the digital I/D filters at time t= kTp = 
k(T b /M). The control signal, after being filtered by the loop filter, 
will produce a signal to adjust the local frequency synthesizer. 

C.3 Discriminator Characteristics and Pull-In Range 

The generation of the control signal (C.3) is motivated by the fact 
that an unnormalized frequency discriminator can be formed by [2] 


• t dQ dl 

9 = 1 dt - «3t 


(C.4) 


Approximating the time derivatives by finite differences over a time 
delay Tp, we have 


e = I(t-Tp)Q(t) - I(t)Q(t-Tp) 


(C.5) 


For I(t),Q(t) the same as those in (C.la) and (C.lb), ignoring 
noise and data modulation, we find 


0 = °a+G)S' sin(AaTp) ^ 2S[sin( AuiTp)] 


Since the I/D filter has a transfer function as 


(AM r 

H( oi) = sinc(-^-) 


(C.6) 


where 
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sinc(x) 


sin x 
x * 


it can be easily shown that the FAL discriminator characteristic (DC), 
which is defined as the mean of the control signal, is 

D( Aw) = 2S sinc^(-7j — ) sin(AoJp) . (C .7 ; 


D(aw) will be normalized with respect to O' (0) henceforth. 

Eq. (C.7) is shown in Fig. C.2, from which we see that the DC is 

Atiil* p 

negligibly small when | — ^ | > 2.0; see also Fig. C.6. The arrows on 

the curve indicate direction of frequency error variation with time in a 
noiseless environment. It is clear that the noiseless pull-in range is 

| AuiTp | < ir 

or (C.8) 

| Af T F j < 1/2 

where Af = Au/2ir. 

In the presence of thermal noise, the pull-in range is of course 
expected to be narrower and dependent on actual input signal to noise 
ratio (SNR j ) . Fig. C.4, C.5, C.6 are three typical frequency error 
trajectories for AfTp = .42, .525 and 1.05, respectively. 

C.4 Simulation Model 

According to the analytical block diagram shown in Fig. C.l we set 
up the following simulation model for the AFC loop; the corresponding 
system parameters and their definitions are summarizeo in Table C.l. 


1^1 i 

N Jl (S k,t +n ka ) 
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Fig. c .2. Phase Discriminator Characteristic for FAL (D(AfT r )/D' (0) ) 



Table C.l. FAL System Parameters and Definitions. 


Bit Signal to Noise Ratio 
Bit Rate 
C} : I Power Ratio 
Sampling Rate 
.oop Bandwidth 

Accumulation Numbe" for the I/D Filter 
UpJ'te Speed 

Control Signal at t * kTp 
Output of the I/D Filter 
at t 1 kTp 

Signal Component of I k (Q k ) 
t = (k-l)Tp+(t-l)T s 
.w se Component of I k (Q k ) 

Original Frequency (Phase) Offset 
Frequency Error at t * kTp 


E b /N 0 ■ ST b /N 0 

1/T b (R b i 

G (1 £ G £ 4) 
I/T s 

1 

N 

Tp (-. d /H) 

e k 

I k (I-channel) 
Q k (Q-channel) 


&( *U ( fy) 
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Ql. = 


- 7 (s Q +n Q ) 
N ' s k,t n k,r 


,1 *1 Q . i 

[fj 1 s k *) + n 

" t=l 


A Q Q 
" s k + n k 


(C.9b) 


s k,Jl = 6 I d I cos t\ + A«j c (A- 1 ) T s ]- G Q d Q sin[ 4| ( + A(ij^( *-l)T $ ] (C.lOa) 


SjJ ^ ijd, sin C^ + Au^(t-l)T s >Gpdg cost <j t +ac^(t-l)T s ] (C.lOb) 
lc k 


wi.ere n^.n^ are i.i.d. Gaussian random variables with variance 
Ng/(2sT f: ) > G? = 1//1+G and G* = /G/l+G. The difference equations 
governing the control signal e k , the local frequency and phase 
estimates <^,6^ ar.d the state variables, I k , and for a first 
order FAL are 


e k = 

^-l^k ' ! k Q k-l 

(C.13) 

“k+1 

= “k + ^k 

(C.14) 

\+l 

= K + *k T F 

(C.15) 

A vi 

= Aa b " “k+1 

(C.16) 

Vn 

Vl ‘ \+l 

(C.17) 

\+i 

= Aujj(k+l)Tp + • 

(C.18) 


C.5 Performance Analysis 

C.5.1 Frequency Jitter: Steady State Analysis 

From (C.14) and (C.16), we obtain 

♦Note the Gj and Gq defined here are square roots of those defined in 
Appendix 8. 


c -3 
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Aa k+1 = Aa k ~ He 


a<^ - -yECe^ * YCe k -E(e k )] (C .19) 


For BPSK signaling, (C.19) becomes 


n+1 = H " l0( H ) + n k 


where 


•YCe k -E(e k n 


(C.20) 


= ~^ n lLl n k“ n k n k-l) + ^ s k-l n k +n k-l s k' s k n k-r n k s k-l^ 


= -YC(hoise x noise) + (signal x noise)] 


(C.21) 


Assuming A&^Tp « 1, i.e., D(ao ) *» AtJp, we obtain 


H+l = H (1 * Y “ F } + n k‘ 


Letting a = 1-yTp < 1, the above equation becomes 


A Vl " aAta k + n k 


(C.22) 


Before we can proceed, a few observations and assumptions are needed. 
First, we notice that 


E(n k ) = 0 


(C.23a) 
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E( Vj ) * 0 


if i * j 


(C.23c) 


Hence {n k } is a sequence of uncorrelated and identical distributed 
random variables. Moreover, Au^ and n k are independent because of the 
independence between the initial frequency error, and rTg. Next, we 
observe that 

ECau^) = Eta^Ao^ + a^ + a* 2 n^ +...+ rT k _j) 

It 

a Au^j 

a 2 E[( a«<^_j) 2 ] + 2aE( Ao^_^n k _p + E(nj^) 

- ,2k. \ 2 . , l-a 2k 

- a (ao^j) + E(n k _ 1 ) ^ 2 . 

The second equation has used two observations mentioned above, i.e., 

E(A VlVl ) * E(Ai t - 1 ) E( V 1 ) " 0 VK 

and 

E(nf) • E(nf) - ... = E(n^. 1 ) . 

Now, if K >> 1 such that a k « 


and 

E[(a^) 2 ] « 


0 then 
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£( A(^} = 0, 

rf / *2-, _ E * n k-1* 

E£( Ao ^ ) 3 • 2 


(C.24) 


and the autocorrelation function is found to be 


E ( Aa»j A(dj ) 


a *£ ( Aa^ ) 

aW.l) 


i = j+t. 


(C.25) 


As has been expected, Au^ is asymptotically independent of a&^, and auj 
and auk are uncorrelated if j i - j I * t » 1. We shall establish (C.24) 
by another approach through which our simulation result on frequency 
jitter is to be justified. Let us form a new sequence {au£{ from the 
sequence { } by the sampling rule: 


Ao£ = AoJj 


i = kj + k. 


(C.26) 


where n and k Q is such that 
E( A<dj Aia; + j )/ ^(ao^)E( Ao^ + j) * 0 and a ^ • 0. 

For analytical simplicity, we assume {n k } are Gaussian so that the 
uncorrelatedness of {n^ } makes them i.i.d. The difference equation 
governing the new sequence |au£} is similar to (C.22): 


a ’K + n k 


(C.27 ) 
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where a* = a J and 
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n k = aJ lfl kj + ^ 2 "kj-l + ** ,+ an (k+l)j-2 + n (k+l)j-l 


* 

i 


i 


It is clear that and n£ are independent and n£ are i.i.d. Gaussian 

1-a 1 2 ~ 

with mean zero and variance equal to ~ * Var(rT-). The uncorrelated 

l - d ^ K 

sequence {&<*}, being driven by a white Gaussian sequence, is itself 
x k 

white Gaussian. This and the choice of kg (a «0) imply that {au£} is a 
stationary white Gaussian process. 

Since has finite mean, by the strong ergodic theorem [3], the 
sample mean 


N K= I +1 Aa k ~ m N (C -28a) 

converges to the E(a<^) with probability one. Furthermore, if {n£} is 
ergodic itself, then the sample variance 


(N-lf 1 


( LU«i> 

K=1 K 



(C.28b) 


converges to Var[aa£] with probability one [3], In summary, if {n^} 
is an ergodic, stationary, white Gaussian sequence then the sample mean 
and variance will converge to the assemble mean and variance with 
probability one. This will be our underlying assumption about {n^} in 
the subsequent discussions. Summarizing au£ from 1 to N and dividing by 
N, we obtain, from (C.27), 

1 N , N i N . 

N - 1 &u k+l = N~ I Aa k + ¥ I n k * 

K=1 " K=1 x n K=1 



-ion. 


For N » 1 and stationary {Au£} 


Or FGor\ Qi 


N N 

l Aa^ +1 “ l Ao£ 

K=1 ^ 1 K=1 ^ 


hence 


m N {1_a) " £ J x "k * 


N 


Now, substituting the identity 


N 


N 


I »Vl = a ‘ 2 I K 2 * 2a' l 4^ a J „'’ 2 

K=1 1 K=1 K K-l K=1 


and (C.29) into (C.28b), we find 


N 


(N-l)s N = (l-a'V l [ I n' 2 a 2a' J 4 ^n' - 1 ( J n' k 

K=1 K=1 


? \ L "k 

N(l-a')^ K=1 k 


Denoting the frequency jitter by o^., we have 


2 , . N 

a. c ~ I im — 5 - 
Ar at a ^ 


N-*-® 4 u 
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1 Am*"' 

lim — N - = 0 
N-h» " L 


Var(nj^) 

4„ 2 (l-a* 2 ) 


(C.30a) 


Var(n k ) 

4* 2 (l-a 2 ) 


4ir 2 (2 yTp-Y^Tp) ^ 2(SNRj) 2 SNR I ^ 


— 5— 1 [.5(SNR: 2 ) + snr: 1 ] 

4ir Tp(2- yTp) 1 1 


'l N rt 


Bp * 1/Tp . 


(C .30b) 


If yTp « 1 then 


JZ r_j_ + _!■■■■ 1 
8n 2 SNR I 2(SNRj) 2 


(C.31) 


From the simulation efficiency viewpoint, however, (C.26) is not a 
good sequence to work with, especially when a » 1, i.e., B L Tp « 1. We 
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would rather use the sequence {D kJ given by 


AUj - aAuij_j = n^_^ 


i = k + Kq, k > 1 


It can be easily proved that 


and 


E^Dj) =0 if i * j 


E <°?> * E ("l + k 0 » * 


Hence o^p can be obtained by first computing the sample variance of the 

white Gaussian sequence {D k } and then dividing it by 4n 2 (l-a 2 ). For 

2 

other modulation format. Op can be derived accordingly by substituting 
the corresponding D(ao^) and Var(n k ) into (C.28). We have found, 
however, that they are all equal except for SPQSK which renders 


A 


1 


V 8i 2 l a^SNRj) 2(a 2 SNR I ) 2 


] 


(C.32) 


where oj and <^> are functions of M( s Tp/T b ) and N(=T S /Tp). When both M 
and N are large, say 10, then « 1. 

So far, we have assumed a NRZ signal format and ignored the data 
transition loss caused by a reverse of the control signal at the bit 
junction whenever a data transition occurred. It can be shown that the 
average loss for a BPSK/NRZ signal is .75/M and 3(1+G 2 )/16M for 
UQPSK/NRZ. Because the Manchester ceding guarantees a sign transition 
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during a bit time the average loss fo* 8PSK is 3/2M and and 3(1+2G <: )/8M 

for UQPSK. It is the consideration of such a data transition loss 

whence a median (>_ 5) or large M requirement is derived. A more 

2 

accurate estimation of should replace SNRj defined in (C.28b) 

by SNRj - DgSNRj, where Dq is the data transition loss described above. 

C.5.2 Acquisition Time: Transient State Analysis 

C.5.2.1 Noiseless Behavior 

The noiseless acquisition time, T^gp, is defined here as 

T = inf {t: Aul=Q} , (C.33) 

nyi- t>Q x. 

i.e., the first time the frequency error becomes zero. 

There will be no attempt to address the nonlinear FAL analysis 
here. Throughout this subsection discussions are restricted to the 

linear operating region of a FAL, that is, Aug is always assumed to be 

such that AugTp < 1/6. In this region, the continuous-time version of 
the FAL shown in C.l can be equivalently modeled as that illustrated in 
Fig. C.3. For a first order loop, 

F(s) = sTJ * 

Hence the closed-loop transfer function becomes 

K K F( S ) 

H(S) = l+K r K 0 F(s) 


■ 1 


( C . 34a ) 



i / 
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Noise, 
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v„ - K r (<y%) ± 


F(s) 
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Compared with (C.22), we obtain 


and 


V 1 


K r = a 


K 0 S Y 


Defining the loop time constant as 


'JkIC"’ 1, t3 

r »F POOR QUALITY 


T L = Tj/K r K 0 = l/a Y 


(C.34b) 


and taking the inverse Laplace Transform of (C.34a), we obtain 


iw t 


a -1 ' t/T L 

±oL 1 ( H (s)) = Auge L 


(C.35) 


Therefore T^gp can be estimated via 


t = T. In — ^ 
L A<^ 


(C .36) 


so long as Aug •Aw^. > 0. 

C.5.2.2 Noise Performance 

Eq. (C.33) is used to define the so-called noiseless acquisition 
time because the behavior of {a<^, k=0,l,...} in a noiseless 
environment, being wel 1 -approximated by the deterministic equation, 
(C.35), of its continuous-time version, is a deterministic and 
decreasing function of k. As a noise performance measure, (C.35) is of 
course not appropriate any more, due to frequency jitter. Hence, we 
define the A-P frequency acquisition time as 
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T A Qp(6,P) = inf {|au^.| < 5 with probability P} (C.37a) 
t >0 


Usually, 6 << A«q and sTp is in the linear region of D(a<d), so Gaussian 
assumption is valid and (C.40a) is equivalent to 


t a QF U.P) • < «. erfc(^-) 


(C.37b) 


t>0 


aF 


Even with the assumption of 100% acquisition probability, T A qp(6,Pj 
does not necessarily exit, due to the fact that the inequality 




(C.38) 


cannot always be satisfied. A mo^e feasible performance measure is the 
mean acquisition time, T ft qp = E[T A qp( 6,P)], given (C.38) and 100% 
acquisition probability. An approximation of T A qp is dericed as 
follows. 

In the presence of noise, we have 


d(Ao^.) 55 -•f0( Ao^. )dt + y n(t)dt , 


and in the linear region of D(awj.) it becomes 


d(Au^) s -Y^pAo^dt + y n(t)dt 


where 


and 


n(t) 


t=kT c 


t = kTr 




= "vY 
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If we assume n(t) is an additive white Gaussian noise with a one-sided 
power spectral density N Q , where 


2Var(n ) N ? N> 

N 0 = T = 4T F + 2sT7 1 

Y Bp r F 


then { Aoij. , t >0 } becomes the well-known Crnstein-Uhlenbeck proress [4J 
which is also an approximation used to compute the mean pull-in time for 
the PN code tracking loop discussed in Sect*: a Y. The mean acquisition 
time, assuming a 100% pull-in probability, is thus to be evaluated by 


T AQF = / tP(t, A(^j,0)dt - / tP(t,A<^,0)dt (C.39) 


wher a 


T is a large number 


P(t,4 V y) - ♦(fj (l-e‘ 2ct },i(^e" at ,y), 


♦<t,*,y) - -L .-«»*«» /» , 


a - Y * T~ 


and 


a 3 Y^c • 


Eq. (C.39) is in fact a conditional mean only. The 'real' T^p does not 
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inborn 


exist, followed by a similar argument given in Sec. C.5.1. In case of 
SQPSK, is replaced by 


N 0 * 2T F t ^(SNR.) + 


2( o^SNRj )' 


] • 


(C.40) 


This complete our linear analysis of the FAL for DG2. 

C.6 Performance Results 

C.6.1 Simulation and The 'etical Predition 

Judging from (C.29), frequency jitter is controlled by the triple 
(y.Bp.SNRj) or equivalently by (B L ,Tp,E b /N 0 ,M) since 


Y = 4B L , 

B F = 1/Tp, (C -41) 

ST. T r i 

SNRj = |j — = |j— (jjj-) • (C .42) 

1 "o b 0 


The above fact suggests that a small B|_ and a large Tp (small M) be 
used. On the other hand, a larger Tp implies a smaller pull-in range 
and a small B L will delay the acquisition time; see eq. (C.8) and 
(C.35). One solution to this dilemma is to apply a multiple-stage 
strategy. More specifically, (B^Tp) will not remain constant but will 
vary with time, of course, we would like to have a large B L and small Tp 
initially, and B L (Tp) be switched to smaller (larger) values at some 
predetermined time. Such a multiple loop bandwidth and update speed 
system allows us to enjoy the advantages of both fast pull-in time and 
small frequency jitter in a wide dynamic range environment. 8efo» 
givim; our design proposal based on our worst-case assumption for DG-2 
we like to go through a verification-by-simulation test on the analysis 
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given in the above two subsections. 

To begin with, let us demonstrate three typical frequency error 
(aoj^) trajectories in Fig. C.4, C.5 and C.6 for AfgTp = .42, .525 and 
1.05 respectively. Notice that each trajectory converges to a stable 
operating point according to the directions of arrows shown in Fig. 

C.2. Nevertheless, this is not to say that given an initial operating 
point the dynamic behavior of Au^. will follow Fig. C.2 and converge to a 
predetermined stable operating point with probability one. Shown in 
Fig. C.7 is another trajectory for AfgTp = .525, which turns out to be a 
successful pull-in. These all have been predicted by eq. (C .37a) which 
tells us that the trajectory of Awt is drifted by a time-varying force 
-t0(awj.) and, at the same time, perturbed by the variance term yn(t). 

Roughly speaking, whenever the drift is overcome by the variance, 
which occurs with a nonzero probability in the presence of n(t), the 
trajectory will not follow the phase diagram shown in Fig. C.2 
anymore. Simulation results show that successful pull-in can be 
achieved with high probabilities for E b /Ng ^ 8 dB at 1/T b = 1 KHz if 

| Af 0 T p | < 1/4 (C.43) 


which is larger than the linear region of D(aw). 

As mentioned before, both frequency jitter and acquisition time are 
in proportion to B^. Fig. C.8 is the sample variance trajectory 
computed from Fig. C.4b. The horizontal straight line is the frequency 
jitter estimated by (C.31). The simulation result is very close to the 
theoretical computation. By eq. (C .34b ) and (C.35), it is predicted 
that at t * T l , the frequency error will be reduced to e“ A of the 
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Fig. C.8. Sample Deviation Trajectory Computed From Fig 
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original frequency offset Aug. Fig. C.9 and C.10 exhibit two different 
cases: in Fig. C.9, AfgTp = .42 while in Fig. C.10 AfgTp = .042. 

Obviously, the prediction based on (C.35) is too optimistic for the 
former but quite accurate for the latter. The reason for the 
discrepancy in Fig. C.9 is simple: the operating point for the former 

is in the nonlinear part of D(au) and the derivation of eq. (C.35) 
assumes a linear D(au), which implies a much larger D(au) at AfgTp = 
.42, and zero noise. 

C.6.2 Design Procedure for DG-2 

Our design proposal is based on the assumption that the available 
£ b /N 0 at 1/T b = 1 KHz is at least 8 dB and the largest frequency 
uncertainty is ±4.2 KHz. In order to operate in this worst-case 
condition, Tp must be such that jTpAfnj < .25, which implies 


T p < 6.0 x 10' 5 (sec) 


or 


(C.44) 


M > 4 KI T b - 17 

Our second requirement is that noiseless Aube reduced to within ±25 Hz 
in 0.2 second. From (C.36), the loop time constant is 


Auh -1 


Since = 1/4B L the required loop bandwidth is given by 




(C.45) 


Substituting t = 0.2, a«q = 4200 and Auj. = 25 into (C.45), we obtain 
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B l » 6.4 (Hz) 


The third design criterion is to have a frequency jitter less than 10 
Hz. By (C.31) this requirement is translated into 


B 


2 ( 110 ^) T f 


L C(SNRJ)" 1 +(2SNRJ 2 )* 1 ] 


(C.46) 


If M = 20 then we require 


B l » .012 (Hz) 


and if M=5 then B L » .36 (Hz). 

With Bj_ = 6.4 Hz, M * 17 the frequency jitter is about 186 Hz 
according to (C.31). Assuming 100% pull-in probability and a Gaussian 
distribution of af at t = 0.2 sec, the above jitter implies then that 
|af| < 211 Hz (one a point) with probability 0.68 and |af| < 30 Hz with 
probability 0.02 only. An equivalent statement is that |af| will remain 
large for a very long time; see Fig, C.ll and C.12. To get rid of this 
undesirable phenomena at least one intermediate stage is needed. If M 
and B[_ in this stage is such that the worst af, assuming to be the 3a 
point 583 Hz at t = 0.2 second, be reduced to 5 Hz (noiseless) when t = 
0.6 second. Again by (C.36) we obtain 


B L = 2.974 (Hz). 


Using M = 5, o^p is found to be 28 Hz. 

In addition to the above three considerations, the accumulation 
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number N must be large enough so that (1) the loss due to the data 
transition is negligible, and (2) the digital implementation of the 
I/O filter has a performance close to an ideal analog I/D low-pass 
filter. For this, we demand that 


N > 10. 


The above seemingly ad hoc desiyn approach can be stated in a more 
systematic way, i.e., as a mathematic programming problem, as follows. 
(1) Liven T b (R b ), E b /N Q and &f Q 
=£» decide initial stage M via 


M > 4 Klv 


(2) Given a 


a F 

decide final stage B^ and M by 


B 


2Uo /sp ) T p 


and 


L C(SNR I )' 1 +(2SNR^)' 1 ] 


SNRj = (^)/M 


(C-44) 


(C.46) 


(3) Given the number of intermediate stages, N $ 

-r-^-decide the corresponding Bj_ and M such that T A qp(6,P) 
defined by (C.37b) is minimized. 


Our approach for solving (?, described above is mathematically 
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unsound and clearly incomplete. An estimate using Gaussian and 100% 
acquisition assumption shows |Af| < 46 Hz with probability .9 at t = 0.6 
sec. However, simulation results (Fig. C. 13-14) are better than what 
has been predicted perhaps because of the discrete nature of the AGC 
loop. | 

For high E b /Ng, two-stage loop is often good enough and the design ■ 

procedure can be summarized as that in Table C.2. 

C.7 Lock Detector and Its Performance 

Lock detector (LD) is an auxiliary circuit used to provide a 
recognition signal when the frequency error is reduced to a preset 
level. Conventional LD algorithms for Costas or squaring loops do not 

i 

work for FAL because of either large pre-sampling bandwidth 
2 2 

(hence I k +Q k fails) or inappropriate indication signal 
2 2 

(since I R -Q k a cos <^). 

A good LD for FAL has to take into account the fact that the AGC 
may not be able to respond appropriately when the frequency error is 
large. So the design philosophy is to avoid comparing with a signal- 
power-dependent threshold. This consideration lead to the following LD 
algorithm: Define 


C k ' Vk-l + Q k Q k-l 


(C.47) 


Accumulate C k and e k for Mg times, i.e.. 


k 0 +M 0 

‘i - i ft I 

k= V i 


(C.48a) 
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Table C.2. FAL (2-stage) Design Procedure for High E b /Ng. 


SYSTEM PARAMETER 
First Stage Tp (M) 

First Stage B L 

Second Stage Tp(M) 
and B l 


DETERMINED BY RELATED DESIGN 

PARAMETER 


T f ' «KT • 
H>4|rf 0 |/R b 

M > 5 


l Af ol " 

Pull-In Range in Hz 
Data Transition Loss 


B l = ( 4t ) ~ 1 



= 2(i I a AF ) 2 T F 
L [(SNRj)' 1 +2(SNR|)" 2 ] 

SR "1 - (")' 1d o 


t - Noiseless Pull- 
In Time (sec.) 

= Pull-In State 

= Design Frequency 

Error at t 

°aF = Frequency 

Jitter (Hz) 

E 

Tj- = Bit Signal-to- 
N 0 Noise Ratio 


Accumulation Number 


> 10 


D n = Data Transition 






Then declare acquisition if 


X i > x 0 • 


GK!G:;- 

OF PC- 


In (C.48), k Q is chosen to be ’arge enough that e^ and are in steady 
state. The relation between (C^.e^) and (C k ,e k ) is the same as that 
between ag£ and Au^ ; see eq. (C.26). It can be shown that 


AwT, 


o UUJI r 

E{C k ) = 2S sinc^^-^Jcosl Au^Tp). 


(C.49) 


Thus for AioTp « 1, 


E(C k ) = 2S cos Ao^Tp 

Let C k be written as 


C k = E ^ C kl Au k^ + *- C k‘ E ^ C lkl 


A 


Z + C, 


k* 


It can be shown that and C” k are uncorrelated and so 


E(C. ) = E(C. ) + E«T. ) 



h 1 


. * 


! 


[ 
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2 2 
-Af/2o 


aF 


d(Af); 


aF 


moreover. 


El# 


E(C*) + E(Cf ) 


= ^=— / E 2 (C k |Af)e 

/o n ^ 


/ 2 ito aF 0 


Af 2 /2o 2 At : B F 

tf d(Af) ♦ B v (^-2) 


- a. 


(C.50) 


where 


.A o/°aF^2 


«v ' 2 hf ) 


and E(C k |Af) is given by (C.49); replacing a<d by 2uAf and x by 1. We 
notice in steady state |AfTp| « 1; thus C^- 1 and |^ k | » 1. This can be 
justified as follows. 

The assumption that | AfTp | « 1 in steady state with a high 


probability is derived from the small frequency jitter and large Bp 

2 


(small Tp) assumptions. Under these conditions, i.e., a^p/Bp « 1, Bp 
« 1, 


,tt 2 


E(C‘) - /— §— / cos (2itAfT F )e 




d( Af ) 




and ^ « 1, SySp « 1, The right hand side of (C.52) is apparently 

O 

dominated by the first term, EfCjp; so C k * C k and 
E(|C k |) - E(|CJ) 


r-R — « -Af / 2a c 

-✓-= 2 “ / cos(2nAfTp)e ^dUf) 

0 



OKiv- 

OF FC ■ 


As for e k , the same assumption leads to 


e k - 2nTp Af + e k 


(C.51) 


where e k is defined ,n a way similar to C^. Based on (C.7) and the fact 
that Af being a Gaussian variable, e k written as 


e k + e k 


has a mean E(e k ) 


E(e k ) * E(e k ) + E(e k ) - 0 


and variance Var(e k ) 


Var(eJ = E(e*) + E(ef) 



'jyifr ■ 


cjCinC^t 


om 




) — 9 ~ " ? i -Af^/2 a^r ®F 

■ /Tl^kl^ dUf) + Mlf •) 

"°aF 0 L 


Hence 


(x.52) 


, , 0 -x -*W 0 - x -x^/2ay 

E(|eJ) - / ~~ e ° dx + / "ZZ e d 


•» /2 7f a n 


0 /2 na n 


.. /ii # 


(C.53a) 


Var ( l e k l ) = % - f ^ = C 1 - \) i 


2, l A “2 

J °0 “ °0 


(C.53b) 


Notice that by taking absolute value of e k , the variance is reduced 

by a factor of (1-2/n)”* » 2.75, Var(e k ) is of the order of (STc/N n )"* 

2.7SST . 

while Var ( | e k | ) = 0(( — ^ -)" 1 ). Similar reduction between Var(C k ) and 

Var(jC k |) is also obtained. Define 


- £ (| c icl) 


o « o -&f^/2a^p 

— / sine ( wAfT F )cos(2irAfT F )e ^d(Af) (C.54a) 

✓T o r r 


a 2 2 2 
°1 = °1 “ U 1 * 


oCinC^oi 


a* ^ *»=• 
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For Mg >> 1, Xj and Xg can be approximated by Gaussian random 

variables, using a central limit theorem type argument, with means Mg^, 

2 2 

Mggg and variances Mg a,, MgOg respectively. 

The detection probability of the ID is then given by 


\ ' Pr lY x o) 


■ Pr(x r x 0 >0| 


= erfc( 




/,-tr 2 , 


l 


✓ (of+cJJ/Mg 

> 1/2 (C.55) 

The inequality, Pg^ > 1/2 is derived from the fact that > wg which is 
obvious from observing (C.53a) and (C.54a). Of course, Pg^ is still a 
function of SNRj, though a very insensitive one (Fig. C.15), but the AGC 
sensitivity problem has been eliminated. 

C.8 Modification for DG-1 

0.8.1 Modification in Loop Structure 

PN sequences are used in DG-1. In the receiver, code 
synchronization must be estaolished before frequency acquisition 
subsystem is initialized. For mode 1 and 2, the I and Q channel 
modulating PN sequences are derived from the same generator but the Q 
channel sequence is offset by x + 1/2 chips relative to the I channel 
sequence where x >_ 20,000. It is appropriate to model them as two 
independent PN sequences as far as performance is concerned since Tp/T c 
« 20,000 in all practical cases. To utilize signal energies from both 
independent channels, like the 2-DDl case we 'enhance' the FAL of Fig. 
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C.l to that of Fig. C.16. Statistically, the control signal from the 
upper pair of channels is the same jS that from the lower ones, except 

y 

for a scale factor. It is thus proper enough to discuss the upper ha'.v 
only. Unless otherwise specified, it is henceforth understood that the 
signal and noise under consideration are from the upper (I) channels. 

To simplify our notation, the argument for distinguishing upper and 
lower channels will be suppressed if there is no danger of confusion. 

For this 2-channel FAL, (C.IOa) and (C . 10b) should be modified to become 


S k,»“ G I d I k cosC «k * n ( ‘- 1)T s ] * "It 


4 C ‘I *1 

* h,i n k,i 


(C . 56a ) 


where 


S k,t = G l“l k 5lnC ^ " + "k , : 


t $ * 

n k ,£ 


(C.56b) 


"k.t = "n pN G Q sin[ VH (t-l}T s J 


Vl = n PN G Q cosC VH lt - 1)T s 3 

T s T s 

n pN »}-/ d(t)PN 1 (t)PN 2 (t)dt = j- dPN 1 (t)PN 2 (t)dt 


(mode 1 and 2 or mode 3, Bi-^) 


(C .57a) 


1 ^ 

4 - / d(t)PN, (t)dt 

s 0 1 


1 5 

4 - f dPN, (t)dt 
's 0 1 * 


(mode 3,NRZ) 
(C.57b) 


- 226 - 


n. r I 

oLinL^om 

S> 
















d(t) is a +1 valued data stream, d is a +1 binominal random variable, 

T 

and PN^(t), PN2 (t+ ) are tw0 independent PN sequences offset by half 

a chip time (T c ). 

C.8.2 Noise Model for DG-1 

Mathematically, eq. (C.57a) and(C.57b) are equivalent to 


OK1G''aA'i- ■ ■ - 

OF POOR C- 


'PN 


2T s /T c b. 


l T- 


i=l 


1 VT ‘ 
■f l b i 

L i=l 1 


(C.58a) 


(C.58b) 


where b^ are independent equally likely +1 valued binominal random 
variables and T s /T c - L is assumed to be an integer. We can rewirte 
(C .9a ) and (C .9b) for DG-1 as 


1 N 


(s "k,t +n l t tn k,t> 


and 


" S k + N ^ n k,ji +n lc, ^ 


- S^n 1 
k n k 


(C.59a) 


Wu + i 1 ‘"2.. 

t-l 


£*1 


'k,t k,4' 


^ S? + n? 


(C . 59b ) 
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It is well known that n PN approaches to a Gaussian variable as T s /T c 
goes to infinite. For finite T s /T c , tan be shown [Appendix I, 5] 
that the rms approximation error e, de ned by 


l [P h (k)-P n (lc)]‘ 


where 


P b (k) * Pr(n pN =K) = ( k+l )2" l , -L < K < L 

~2~ 


V k) - \ [ er U^)-erf(^)] 

9 c /2L V2L 


is less than 3% if L 12 for mode 3 or if L ^ 6 for mode 1 and 2. 

Based upon this, n* and nj^ will be approximat’d by two independent zero 
mean Gaussian variables with variance equal to 


Var(nJ) = j 


1 

>L'il+( 


♦sB 


(mode 3,NRZ) 


(otherwise) 


(C.60a) 


Var(njJ) = Var(nJ) 


(C.60b) 


L* = Ml > 6 (mode 1 and 2 or, mode 3, Bi-$) 


> 12 (mode 3,NRZ) 


(C.70) 


C.8.3 Frequency Jitter 

Forming the control signal e k by adding the error signals from the 


/ 
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upper channels e£ and that from the lower ones e£, i.e.. 


e = e u + e L 
e k e k e k 


the frequency jitter is found to be 


. '■■■ - • 
OF POOR Q- 


. i 


V 


t 

8 


, 1 r 1+G 2 1 . 1 1 

n 2 1 2 1 a 2 L* 2 (l + G) 2 Ur ^ (SNRj) 2 J 


1+G 


+ T* 


1 


a 0 L* (1+G) 


2 (SNR) 


I 


(C.71) 


where 


a Q = 4 for mode 1 and 2 or, mode 3 with Bi- 


? for mode 3, NRZ. 


For L* » 1, 


2 s ir 1 + _J_ 1 

^ 8 it 2 2(SNR x ) 2 SNR I 


(C.72) 


If T” 1 = 3 M bits/sec, H = 20, 1/Tg = 1 k bits/sec, then L' * 1.5 x 10 2 , 
and if M=5, L' = 6.0 x 10 2 . 

In both cases, L‘ » 1 and 1/L ' « 1/SNR j so the design parameters 
in C.6.2 remains good for DG-1 if the design criterions are the same; 
see Fig. C.17 for a comparison. Suppose only the lower channels, i.e., 

the stronger Q channel signal is tracked, then (C.71) becomes 

2 _ ^F fo -l r 1 . 1+G ,2, 1 , 1+G , , r 7 o\ 

°aF 8i 2 t 2 t ajjL r GTsM^T J a^lT G(SNR“) (C * 73) 
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(dB) 


C.17. Frequency Jitter for DG-1, Mode 1 ( ) and DG-2 
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APPENDIX D , 

INTRODUCTION TO THE APPLICATION OF AIRS DEMONSTRATION PROGRAM (APD) 

I 

D.l Purpose * 

The AIRS demonstration program (APD) described herein is a FORTRAN j 

language software package capable of analyzing and simulating the 
performance of various subsystems of the AIRS. It is an interactive j 

software package developed by LinCom for NASA/ Goddard Space Flight 
Center to run on the Perkin-Elmer computer, and will be used as a ^ 

software analysis tool for the AIRS. This appendix and the one j 

following discuss the organization and applications of the APD. 

Definitions of input and output parameters are documented, and, for 
explanatory purposes, a sample run is included at the end. 

D.2 APD Design Philosophy and Architecture 

The functional guidelines which were followed in the course of 
designing, developing the APD are stated below: 

(1) APD will allow a great deal of flexibility for any system 
change and expansion. 

(2) It will serve as a engineering analysis tool in both system 
design stage and system operation stage. 

(3) It should be easy to understand, use and maintain. 

(4) The outputs should contain all the information necessary for 
system performance assessments. 

(5) The powerful capability of the existing software package-CLASS 
should be utilized. 

(6) The computation time should be kept as a minimum. 

As a consequence of the above requirements, the APD is organized in 
a logical and structural modular form as shown in Fig. D.l, basic 
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Figure D.l, Current Architecture of the AIRS Demonstration 
Program (ADP). 
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routines referring to those used in more than one configurations within 
one subsystem and/or in several subsystems. Detailed architectures for 
the PN acquisition system are shown in Fig. D.2 and Fig. D.3. The PLL 
and bit sync simulation package are discussed in Appendix E. 

In order that the PAD be able to communicate with and thus utilize 
the CLASS an interface routine is included. 

D.3 Properties and Parameters 
D.3.1 The Acquisition Subsystem 

There are fourteen system parameters of the Acquisition Subsystem 
to >e specified: 

(1) C/N q (dB-Hz) 

(2) Data format (NRZ or bi -phase) 

(3) Data rate (kHz) 

(4) Code rate (M chips/sec) 

(5) I and Q channel power ratio G, 0 < G < 1 

(6) Designed false alarm rate 

(7) Search step size (chip) 

(8) Doppler loss (dB) 

(9) Data Transition loss (dB) 

(10) Number of Doppler bins 

(11) Number of chips per Doppler bin 

(12) Length of the CCD PNMF (chips) 

(13) Number of the noncoherent integrations (1st stage) 

(14) Integration time (K chips) fo the 2nd stage test, three options 
to be chosen from: 

(1) Receiver configuration (SSP or SPS) 

(2) Optimization with respect to parameter 13 above 
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(3) Performance criterion - minimizing the me-.n acquisition 
time or, minimizing the time for 90% acquisition 
probability. 

Note that the minimization is carried out by choosing the optimal 
detection (& false alarm) probabilities of both stages. For the second 
criterion, a 90% detection probability per cell is assumed. Besides the 
evaluation of the acquisition time the APO also offer operation 
characteristic analysis for various detector structures. 

D.3.2 Code Tracking Subsystem 

Though a 2-DDL is proposed and found to be superior to other 
existing code tracking loops, the performances of DDL and TDL are also 
presented. The effects of the pre-despreading RF filtering and IF 
filtering in both I & Q arms of the tracking loop are considered and 
modeled as signal energy degradations. Lock detector is part of the 
code tracking system and its performance evaluation is included in the 
task of APD. System parameters pertaining to the code tracking are 
summarized as follows. 

( I ) - ( 6) : same as (l)-(S) & (7) of the acquisition subsystem. 

(7) Loop bandwidth (Hz) 

(8) Cut-off frequency (normalized by data rate) of the arm BPF 

(9) Chip time - bandwidth product of the RF filter 

(10) Phase nonlinearity of the RF filter 

(II) Designed fa’se alarm rate for the lock detector (LD) 

(12) Integration time (chips) for LD. 

Two options dc exist, namely, 

(1) Loop configuration (2-DDL, DDL or TDL) 

(2) Choice of the IF filter (ideal BPF or 3-pole Butterworth) . 
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Code phcse jitter, normalized mean slip t me and mean pull-in time for 
the ?nd order loop are calculated, where a 2 dB degradation from the 1st 
order loop in jitter square is added to the 2nd order loop. One-delta 
loops, i.e., one chip delay between early and late gates, and zero 
original code phase offset are assumed in all computations. 

0.3. 3 The Frequency Acquisition Subsystem 

A digital 2-channel, multiple-stage, 1st order automatic frequency 
loop is investigated. The maximum number of stages in a AGC loop is 
limited to be three. Besides ( 1 )- (3) and (5) in D.3.1, there are eleven 
parameters to be determined, namely, 

(5) Original frequency error (Hz) 

(6) Number of samples per update (N) 

(7) Number of upuates per 'ic time for the 1st, 2nd & 3rd stage 

(8) Loop bandwidth (Hz) for the 1st, 2nd & 3rd stages 

(9) Switch time (sec) from the 1st (2nd) to the 2nd (3rd) stage 

(10) Seed (a double-precision leal) for random number generator. 

To control the frequency error trajectory to be plotted, three 
parameters are added: 

(11) Number of points (at most 100) in a curve 

(12) Sampling interval, i.e., number of updates between 2 points 

(13) Plot-starting time (sec). 

By using proper values in (11)-(13), we can examine closely both 
"global" and "local" behaviors of the frequency acquisition process. 
Parameter (10) is used to generate different and independent 
trajectories by entering different values. 

0.4 Accessing and Running APD 

To access the APD simply type the system command 
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*RERUNRFI/G AIRS 

and APD will respond with a series of options-and-parameters sections to 
guide the user, in a page by page procedure to obtain the numerical 
resulcs he or she needs. To terminal APD, a user needs only to choose 
the “ENDING (APD, OR PN ACQUISITION, etc.) ANALYSIS*' command in the 
control routine. The application of the APD in designing an optimal 
system f . s a standard recursive process as shown in Fig. D.4. All 
the related system parameters are determined in at most three sections, 
where each section is designed in such a way that it will fit one 
screen. T facilitate the usage of APD, all system parameters are 
defaulted to some values and users can have options to change or not 
change any of them for their design purpose. For a better understanding 
of the APD we include below a sample run of it. 

i 
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OF POOR QUALITY 


AIRS DEMONSTRATION PROGRAM 
(AGP) 

VERS I ON J..0 
LIMCOM CORE. 


■ HE CURRENT DATA CHARACTERISTICS ARE : 

’ . PN USE = 0 

0 => I ?<© ( DG1-1 ,2) 

1 => I ONLY (DG1-3) 

2 => NONE (DG2) 

SYMBOL FORMAT (OR PIT FORMAT IF UNCODED) 

0 -> B I -PHASE 

1 => NRZ 

2 => UNKNOWN 

MODULATION FORMAT = 3 

=> BPSK ( DG2 ) 

1 => SCPSK f DG2) 
i => BPSK *2 (001 ) 

4 **> UNKNOWN < DG2 ' 

CODE RATE = 1 

0 => UNCODED 

1 => RATE 1/2 

2 *> PATE 1/2 

3 *> UNKNOWN 

fc. I BIT RATE (KHZ > * 1 

7. Q BIT RATE (KHZ) = 1 

3. Q/I POWER RATIO (BETWEEN 1 & 4) * 4 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 
ENTER 0 IF NO PARAMETER CHANGE IS TO BE MADE 

C MTER NEW VALUE FOR PARAMETER 1 

i 

* 

^NTER THE PARAMETER NUMBER TO BE REDEFINED 
E M TSR 0 IF NO PARAMETER CHANGE IS tO PE MADE 

'WE CURRENT DATA CHARACTERISTICS ARE : 



ORIGINAL PAGE IS 
OF POOR QUALITY 


D * 

Oif 


PN UPS ~ 1 

0 ! i-C i DO 1-1 t 2 ) 

! => { :. , UL V (DG1-3' 

2 *»> NOME (DG2) 

SYMBOL FORMAT (OR BIT FORMAT IF UNCOE'ED ) 

0 => ? I -PHASE 

1 => NRZ 

2 => UNKNOWN 

MODULATION FORMAT * 3 

O -> BPSK < DG2) 

! => SQPSK ( DG2 > 

3 => BPSK*2 <DG1) 

4 *> UNKNOWN ( DG2 > 

CODS RATE » 1 

0 **> UNCODSD 

1 *> HATE !/•:> 

2 *> RATE 1/3 


3 => UNKNOWN 

.. I PIT RATE (KHZ ) = 1 

' 0 BIT RATE (KHZ ) » 1 

8. O/I POWER RATIO (BETWEEN 1 &4) * 4 

ENTER 1 FOR PN ACQUISITION SUBSYSTEM ANALYSIS 

2 FOR PN TRACKING SUBSYSTEM ANALYSIS 

3 FOR FREQUENCY ACQUISITION LOOP ANALYSIS 8* SIMULATIONS A8tS> 

4 FOR BIT SYNC A?<S 

5 FOR PLL A&S 


m "V i 


♦ 

1 
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A'l C'L 1 1 S IT I ON SUBSYSTEM ANALYSIS 


ORIGINAL PAGE l« 
OF POOR QUALITY 


SN T ER 1 FOR ACQUISITION TIME EVALUATION 

2 FOR RCVR. OPERATION CHARACTERISTIC INVESTIGATION 


THE FOLLOWING IS A LIST OF PARAMETERS TO BE DETERMINED 
INPUi C/NO(DB> = 36.5 

4. PM C ODE RATE (M-CH I PS/SEC. > «= 3.00 

5 . ACQUISITION ALGORITHM = ;> 

1 => SPS( SERIAL SE ARCH) 

J => SSP( PARALLEL SEARCH) 

6. DESIGN FALSE ALARM RATE = 0. 10D-04 
SEARCH STEP SIZE (CHIP) = .50 

•*:. DOPPLER LOSS(DB) * .40 
•. DATA 1'RANSM I SS I ON LOSS(DB) = .20 
10. n OF DOPPLER BINS * 3 
t!. n OF CHIPS PER DOPPLER BIN * .3SD+04 
16. # OF NONCCH. 1NTEG. FOR THE CCD DETECTOR * 500 

1 / . LENGTH OF CCD PNMF(CHIPS) = 511 

IS. INTEGRATION TIME(K-CHIPS) FOR THE 2ND STAGE TEST = 100. 

IS. AUTO-OPTIMIZATION ON THE 2ND STAGE TEST = 2 

1 => YES 

2 => NO 

20. PERFORMANCE CRITERION = 1 

1 => MIN. THE MEAN ACOUSITION TIME 

2 *> MIN. THE TIME FOR 907. ACQUISITION 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 

ENTER 0 IF NO PARAMETER CHANGE IS TO BE MADE 
« 

ENTER NEW VALUE FOR PARAMETER IS 
4 SO 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 
ENTER O IF NO PARAMETER CHANGE IS TO BE MADE 
19 

ENTER NEW VALUE FOR PARAMETER 19 
• « 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 
ENTER 0 IF NO PARAMETER CHANGE IS TO BE MACE 
16 

ENTER NEW VALUE FOR PARAMETER 16 

A f ’O 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 
ENTER 0 IF r.""t PARAMETER CHANGE IS TO BE MADE 
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rMC CURRENT ACQ. SUBSYSTEM PARAMETERS ARE ! 


INPUT C/NO( OR ) = 36.5 

*N CODE RATE! M-CH IPS /SEC. ) = 3.00 
ACQUISITION ALGORITHM = 2 
1 => SPS( SERIAL St ARCH) 

3 => SSP < PARALLEL SEARCH) 

6. DESIGN FALSE ALARM PATE « O. J0D-04 
•. SEARCH STEP SIZE (CHIP) = .50 
. DOPPLER LOSS ( DB ) = .40 

DATA TRANSMISSION LOSS (DP) » .20 
’O, 4 OF DOPPLER BINS - 3 
' 3 . 4 OF CHIPS PER DOPPLER BIN * . 333+04 
6 . n OF MONCOH . INTEG. FOR THE CCD DETECTOR -*= 600 

: /. LENGTH of CCD PNMFI CHIPS) « 5 u 

■3. INTEGRA F I ON TIME(K-CHIPS) FOR THE 2ND STAGE TEST = +50. 

AUTO-OPTIMIZATION ON THE 2ND STAGE TEST = 1 

1 => YES 

2 => NO 

•0. PERFORMANCE CRITERION = 1 

t *> MIN. THE MEAN ACOUSTTION TIME 

2 => MIN. THE TIME FOR 907. ACQUISITION 

-NfER t FOR PRINTING OUT RESULTING PERFORMANCE 
2 FOR TABULATED PERFORMANCE VS. PARAMETER 
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ORIGINAL PAG£ il* 
OF POOR QUALITY 


'cJlindom 


PERFORMANCE WITH THE ABOVE SPECIFICATION CAN BE SUMMARIZED 
FOLLOWS. 


PD 

3= 

0. 59529D+00 

WHEN 

T 

at 

0. 24308D+01 

SEC. 

PD 

a 

0. 83621D+00 

WHEN 

T 

a 

0.48*150+01 

SEC. 

PD 

at 

0. 9337 10+00 

WHEN 

T 

a 

0. 72922D+01 

SEC. 

-n 

a 

0. c '7.:-:l7D+00 

WHEN 

T 

a 

0. 9722SD+01 

SEC. 

-D 

a 

0. 939) 4D+00 

WHEN 

r 

= 

0. 121540+02 

SEC. 

:-D 

a 

0. 995MD+00 

WHEN 

T 

a 

0. 14584D+02 

SEC. 


a: 

0. 9PR22D+00 

WHEN 

T 

s 

0. 17015D+02 

SEC. 

PD 

a 

0. 99928D+00 

WHEN 

T 

=r 

0. 19446D+02 

SEC. 

pO 

at 

0.99971D+00 

WHEN 

T 

a 

0. 2187&B+Q2 

SEC. 

PD 

a 

0. W9S8P+00 

WHEN 

T 

a 

0.24307D+02 

SEC. 

->D 

a 

0. 9999550+00 

WHEN 

T 

a 

0.26738D+02 

SEC. 

“D 

a 

0. 99998 D+00 

WHEN 

T 

a 

0.29168D+02 

sr\ 

PD 

a 

0. 999990+00 

WHEN 

T 

a 

0.315990+02 

SEC. 

PD 

a 

0. 1 OQOOD+O J 

WHEN 

T 

a 

0. 24030D+02 

SEC. 

PD 

a 

0. 10000D+O1 

WHEN 

T 

a 

O. 364*0D+02 

SEC. 


1EAN ACQ. TIME* 5 0.408340+01 SEC. 

VETEC. PROB. = 0. 595290+00 
WELL TIME ON THE 2ND STAGE TEST = 0. 15067D+00 SEC. 
-ALEE- ALARM PROB. < 1ST STAGE)* 0.4047:0+00 
FALSE-ALARM PROB. (2ND STAGE)* 0 . 247091 - 0 '' 
FALSE-ALARM PROB.* 0. 100OUD-04 
DETEC. PROB. (1ST STAGE)* 0.59529P+00 
DETEC. PROB. (2ND STAGE)* 0.100000+01 
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ORIGINAL PAG£ !1P 
OF POOR QUA 


- M ' qq ; cp e:MDfNO ANALYSIS. 2 FOR CONTINUATION 


, *#**#*->#*****#*##*#*##*##*#****■*#■*■«■•*•»■•*#■*•*■*■#■■*■«••«•■¥ ********* 

'••!;=■ CURRENT COTA CHARACTER! ST I CS ARE : 


1 pAJ t I-- C' — 1 

0 => 1*0 (DO 1-1, 2) 

1 «> i ONLY (001-3) 

2 => NONE ( DG2 ) 

. . SYMBOL FORMAT ( OR BIT FORMAT IF UNCODED) » 1 

A — U f _C'i_i^ger 

1 -'.: t)R2 
> -;■• unknown 

:-'. MODULATION FORMAT = :? 

!.» •=.:■ yRSK (DG2) 

1 -> SQPSK < DG2 ) 

:j> => APSK»2 < DO 1 > 

4 => UNKNOWN ( DO 2 ) 

: . CODE RATE « 1 

0 => UNCODED 

1 --> RATE 1/2 

2 => RATE 1/3 

3 =‘> UNKNOWN 

I BIT RATE (KHZ) = 1 

/. Q BIT RATE (KHZ) * 1 

0. Q/I P^'WER RATIO (BETWEEN 1 ft 4) = 4 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 
ENTER 0 IF NO PARAMETER CHANGE IS TO BE MADE 


ENTER 


1 FOR PN ACQUISITION SUBSYSTEM ANALYSIS 

2 FOR PN TRACKING SUBSYSTEM ANALYSIS 

3 FOR FREQUENCY ACQUISITION LOOP ANALYSIS ?< SIMULATIONS! A&S) 

4 FOR BIT SYNC A&S 
■5 FOR PLL. A?tS 



- 248 - 


r 


- * 


s 


D 

B 



AM 


4l 


-« 


0 

0 

0 


0 



ORIGINAL PAGE » 


PM TRACK I MO SUBSYSTEM ANALYSIS P00R QUALITY 

THE FOLLOWING IS A LIST OF PARAMETERS TO BE DETERMINED 
1. INPUT C/NO(D8) « 36.5 

DATA FORMAT (1 FOR NRZ. 2 FOR BI-PHASE) * 1 
.. DATA PAT£<KH7) * 1.00 
*. PN COOS RATE ( M—CH I PS /SEC. ) * 3.00 
SEARCH STEP SIZE(CHIP) « .50 
LOOP BANDWIDTH! HZ ) * 0.5 

.:6. CHOICE OF LOOPS!# OF CHANNELS) * 4 

1 ■»> TAU-DITHER LOOP 

2 *> DOUBLE & DITHER LOOP (DDL) 

4 «> 2-CHANNEL DDL 

??. CHOICE OF BPF « 2 

1 *> IDEAL BPF 

2 => 3-POLE BUTTERWORTH 

IS. CUT-OFF FREQ. (DIVIDED BY DATA RATE) OF THE BPF » 2.0 
29. FRED. (DOPPLER) OFFSET (HZ) * 0.0 

:?1. PHASE NONLINEARITY OF THE PRE-DESPREADING FILTER » 0. 10 
22 . TIME (CODE CHIP) -BW ( PREDESPE AD I NO FILTER) PRODUCT - 2.0 

32 . INTEGRATION TIME(CHIPS) FOR LD * 0. 1000+06 
34. DESIGNED FALSE ALARM RATE FOR LD » 0. 100D-07 
: t . ISO POWER RATIO (G.OCGOl) * 0.250 
ENTER the parameter number to be redefined 
f N r ER 0 IF NO PARAMETER CHANGE If .0 BE MADE 


0 

0 

0 

0 

0 

B 

I 


M 

ENTER l FOR PRINTING OUT RESULTING PERFORMANCE 
2 FOR TABULATED PERFORMANCE VS. PARAMETER 


ocucoRMANCE WITH THE ABOVE SPECIFICATION CAN BE SUMMARIZED 

. -'OLLUWSi 

SIGNAL DEGRADATION DUE TO BPF » -0.94 OB 
UQRELAT I ON LOSS * -1.37 DB 


CODE PHASE OUTER (7.) * 1.069 St 1.199 (2ND ORDER LOOP) | 

NORMALIZED MEAN SLIP TIME ® 0.340D+02 & 0.340D+02 (2ND ORDER ! 
MEAN POLL-IN TIME » 0.9S3D-01 St 0.111D+00 (2ND ORDER LOOP) 
DETECTION PROBABILITY FOR LD * 0.973E+00 
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ORIGINAL PAGE W 
OP POOR QUALITY 


ENTER 1 FOR ENDING ANALYSIS. 2 FOR CONTINUATION 


THE CURRENT DATA CHARACTERISTICS ARE : 

1. PN USE = i 

0 => T?<0 (DG 1-1,2) 

1 => I ONLY (DG1-3) 

2 => NONE < 002 ) 

I.-. SYMBOL FORMAT < OR BIT FORMAT IF UNCODED) * 

0 »> 3 I -PHASE 

1 *> NRZ 

2 => UNKNOWN 

MODULATION FORMAT -» 3 

0 =•> 9PSK (DG2) 

1 *> SQPSK (DG2) 

3 => BPSK*2 (DOT ) 

4 => UNKNOWN <D02) 

4. CODE RATE = 1 

0 *> UNCODED 

1 => RATE 1/2 

2 «> RATE 1/3 

3 *> UNKNOWN 

6. I BIT RAFE'KHZ) « 1 

7. Q BIT RATE (KHZ) = 1 

8. 0/1 POWER RATIO (BETWEEN l ?< 4) = 4 

ENTER THE PARAMETER NUMBER TO BE REDEFINED 
ENTER O IF NO PARAMETER CHANGE IS TO BE MADE 


ENTER 1 FOR PN ACQUISITION SUBSYSTEM ANALYSIS 
2 FOR PN TRACK I N3 SUBSYSTEM ANALYSIS 

.3 FOR FREQUENCY ACQUISITION LOOP ANALYSIS & S I MUL AT 1 ONS ( AfcS ) 

4 FOR BIT SYNC A&S 

5 FOR PLL AfcS 


AnCoi 




imi 1M I a * < 


ro o if only 2 stages are used) 


» OF UPDATE / SIT ftng FOR THE FINAL (3RD) STAGE 
(SET TO O IF ONE-STAGE FAL IS USED) 

LOOP BN (HZ) FOR THE 3RD. STAGE a 0.340 


SNITCH TIM? (SEC. ) FROM 2ND TO- FINAL STAGE * 0,400 

It OF ROT NTS << 100) IN A CURVE * 50. 

SAMPLING INTERVAL! I* OF UPDATES BETWEEN TWO POINTS) 


ENTER 0 IF NO PARAMETER CHANGE IS TO BE MADE 
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ORIGINAL PAGS flf 
OF POOR QUALITY 


■!£ CURRENT FAL. SUBSYSTEM PARAMETERS ARE i 

INPUT C/NO<D B> »» 33.0 

FREQUENCY ERROR (HZ) * 0.420D+04 

# i..F SAMPLES PER UPDATE * 10 

# OF UPDATE / BIT TIME (1ST STAGE) « 17 

LOOP BN (HZ.* ST STAGE) * 6,400 

4 OF UPDATES /BIT TIME (2ND STAGE) = 5 

LOOP BW(HZ> FOR THE 2ND< INTERMIDI ATE ) STAGE « 2.974 

(SET TO O IF ONLY 2 STAGES ARE USED) 

SWITCH TIME (SEC.) FROM 1ST TO 2ND STAGE = 0.200 
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0ATA-AI0E0 LOOP SIMULATION 


E.l Simulation Software Organization 

The data-aided loop simulation software is organized into eight 
main subroutines as shown in Figure E.l. These subroutines are modified 
slightly to accomodate different types of modulations. The main program 
handles the input question and answer session whose purpose is to define 
the simulation parameters. The main program also controls the flow of 
the program execution. The main loop in Figure E.l starts with SAMPLE, 
which generates noise-corrupted signal samples. The sampling instants 
are determined by the timing error of the loop at that point in time. 

The phase error of the loop is used to rotate these samples. The number 
of samples generated correspond to a symbol time. The data symbol as 
well as the noise samples are generated using pseudo random 
techniques. Therefore, the Dal is a strict Monte Carlo simulation. 

The genet .ted samples are processed by DSP to detect the symbol and 
to generate the dynamic bit sync error and the dynamic carrier recovery 
loop error. The dynamic errors are filtered by FIL to control the 
carrier and bit synth synthesizers modeled by SYNTH. The SYNTH routine 
computes the current phase estimate and the sampling clock estimate. 
These estimates are used in the DOPPLR routine, which models the Doppler 
profile, to generate the phase and timing error for the next call of 
SAMPLE. This completes a simulation cycle. 

The five routines SAMPLE. ?«, FIL, SYNTH and u^PLR are designed 
to duplicate the AIRS functv.'S c' a one-to-one basis. Hence, the 
software can be used to realistically model the AIRS operation and 
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predict its performance. 

Aside from these routines, there are three additional ones that are 
used to perform housekeeping functions. Subroutine DETECT keeps track 
of the SNR estimate and the BER. STAT2 keeps track of the means and 
variances of the phase and timing jitter. Finally, PPLOT plots the 
timing and phase error transients during acquisition. 

E.2 Sampled Signal Model 

The incoming signal is represented by 


s(t) = /2P^ dj(t) cos »^t + ^?P^ d 2 (t) sin ^t 

where d^(t) and d 2 (t) represent unity amplitude (+1) baseband symbol 
sequences with data rates 1/Tj and 1/T 2 (t) respectively. This 
discussion assumes that NRZ formats are being used and can be easily 
modified for Manchester format. The received signal is a noise- 
corrupted version of s(t) given by 


r(t) = s(t) + n(t) 


where n(t) is an additive white Gaussian noise (AWGN) with one-sided 
spectral density Nq watts/Hz. Figure E.2 shows the front end processing 
of the AIRS receiver used to convert the received signal to digital 
samples. The received signal is first converted to its I,Q baseband 
components via the coherent demodulation process. The phase tracking 
error of the local carrier recovery is modeled by the phase error +. 

The I-Q ba.-oand components, ignoring the double frequency terms, is 
given by 
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r j(t) * ^ dj(t) cos ^ d 2 (t) sin ♦♦ n^t) 

«*g(t) * djCt) sin ♦♦ d 2 (t) cos ♦+ n Q (t) 

where n { (t) and nq(t) are independent baseband AMGN processes with one- 
sided spectral density Nq watt/Hz. The baseband processes are lowpass 

A A 

filtered and sampled according to the local clocks t fc and t £ . The 

A A 

nominal clock rates of t^ and t £ are integer multiples of the bit rate 
given by Nj/Tj and N 2 /T 2 . 

In Figure E.2, the vectors _x(t) and jrJ(t) denote the two components 
of the lowpass filter outputs for rj(t) and rq(t). Similarly each of 
the vector samples x^ and has 2 components. 

For convenience, the lowpass filter is modeled in the simulation by 
a moving window integrator as shown in Figure E.3. The gain in front of 
the integrator is selected so that the signal samples are nominally 
of unity amplitude level. The lowpass filter for rq(t) is analogous. 
When the bit sync is perfect, the local sampling clock coincides with 
the transmit data clock (defined by the transition insta *') every 

A 

samples. If there is a bit sync error e s t-t, then some the samples 
will not be of unity amplitude as shown in the last example in Figure 
E.3. In this example, the local clock leads the transmit clock so that 
e < 0. 

The components of the vector samples and are given by 
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where 


0R13:xA‘ * 
OF POD A : 


X lk “ ^\k C0S ** *V P 1 ^Zk sin ** "] 


*2k ‘ ^lk s,n ♦ * J 2k cos ♦ + "( 


*u ' -*Fi s u sin ♦ * J zt cos ♦ * n < 


* 2 * * *17*2 "u cos ♦ * d 2 t *i" ♦ * "a 


H t*N/T. 

ff ik ’ T7 /- V t)at 

n 9 t +n,/t, 

d ii * T I* d.(t)dt 

1* >2 t 1 

M* 

for i=l ,2- The noise samples are Gaussian with zero means and variances 
given by 


E[n-f k ] ■ ECS*,] ■ (N 1 /2)(P 1 T 1 /N 0 )- 


Etfijj,] - E[^ t ] • (N 2 /2)(P 2 T 2 /N 0 ) _1 


The sample pair (n^.n^) is independent as well as The 

sample pairs and (n^ »^q t ) are correlated. Notice that the 

samples for the filtered data waveforms are normalized to unity as shown 
in Figure E.3. 

E.3 Loop Parameters 

In the DAL simulation, the loops are characterized by its order and 
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its one-sided loop noise bandwidth. In a digital receiver, the loop 
action can be described by a difference equation* 


Vl * (E-l) 


Presently, we have ignored the Doppler offset as it has no bearing on 
the loop bandwidth characterization. In this equation, F(z) describes 
the action of the loop filter and is given in z-operator form as 


G 2 G 3 

F(z) * 6 ,+—^- + ITT 

1 1 -z 1 ( 1 -z *) z 


For a first-order loop, G 2 * G 3 * 0. For a second-order loop, G 3 = 0. 

A third-order loop »as all nonzero gains. In the simulation, the 
filtering is performed in the subroutine FIL. The error e k is generated 
in the subroutine DSP and is normalized so that ECe^] ~ 4 ^ while 
tracking. 

To determine the loop bandwidth, noticed that (E-l) ' 

approximated by a differential equation 


; * -(Gj/T ♦ G 2 /T 2 s + G 3 /T 3 s 2 W (E-2) 


where T is the time between update of the loop phase error, typically 
equal to a symbol time, and the correspondence (1-z -1 ) = sT is used. 

The loop bandwidth for the analog loop described by (E-2) is well known 
and is summarized in Table E-I. Loop toandwidths used in the simulation 


*We concentrate presently on the carrier-recovery loop. The bit sync 
will be treated later on. 
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Table E.l. Filter Parameters in Terms of One-Sided Loop Noise 
Bandwidth B, . (For a Second-Order Loop, the Damping 
Ratio c = l 0.707 is Used. For a Third-Order Loop, 
the Wiener Optimum Filter is Used.) 


LOOP-ORDER 

FILTER PARAMETERS 

1 

Gj = 4B l T 

2 

Si = § b l t 
G 2 4 (f 8 l T) 2 

3 

G 1 ■ 3 b L t 
g 2 = 4 B l T) 2 

G 3 = \ <4 B L T > 2 
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is defined in this manner. 

For the bit sync, the determination of the noise-bandwidth is 
similar except that (E-l) is replaced by the difference equation 
governing the normalized timing error samples 


■ \- F < z > e k 


(E-3) 


The DSP routine normalizes e k so that ECe^] a for small A^. This is 
defined for 50% transition density with NRZ data. 

E.4 Effect of Data Modulation 

For a single channel operation, i.e., data appearing only on the I 
or Q channel, the determination of equation (E-l) is trivial. In that 
case, T is either the symbol time for NRZ data or half the symbol time 
for Manchester-coded data. For DG2, the situation is a little 
complicated. If the data rates are equal and coherent, then T is the 
symbol time and e^ Is the average of the error signal provided by the I 
and Q channel. For unequal data rates, the loop error signals are 
produced at irregular intervals. Figure E.4 illustrates this point. 

The DAL processes the error signal samples as they become available. 
Therefore, Instead of (E-l) we have two simultaneous difference 
equations, namely. 


♦ktl ■ - F < z > e Ik 

(E-4a) 

■ V F(Z)e Qt 

(E-4b) 


1 

I 


P 1 


I 

; j 

i 

1 ! 


where k and i denotes the Indexes of two different sequences of epochs 
in time. According to the optimum MAP estimation, the ratio of 
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ECej k J/ECepj^] must be equal to Pj/Pq * I:Q channel power ratio. The 
subroutine DSP Is set up to normalized e^ and eq K so that ECe^] ~ 
^/[^(Pq/PjHRq/Rj)] and E[e Q4 ] ~ (Pq/Pj) ♦ t /Cl+(P Q /P I )(R Q / R I )3 for 
small small errors. (We use Rq/Rj to denote the symbol rate ratio of Q 
to I channel.) 

We can combine (E-4b) with (E-4a) by Increasing the update rate of 
(E-4b). Then since In one I-symbol time T, there are only Rq/Rj parts 
of error signal on the average from the Q channel, we have, 
approximately. 


\+l 


♦k - F < z >[ e Ik + e Q 4 ] 


(E-5) 


It Is easy to check that the average value of the bracketed term In 
(E-5) Is proportional to +. Hence, we can use Table I to determine the 
loop parameters. 

E.5 Sample Runs 

Figures £.5 shows a typical run using the command RUN BPSK. 
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ENTER EB/N0(DB>, SAMPLE/BIT, DATA RATE (BPS) 


II > 5 , 8 , 2192 

*• ENTER DOPPLER(HZ), DOPPLER RATE (HZ/S) 

> 0, 765 

f ENTER LOOP ORDER AND BN FOR PLL 
I. >3,70 

ENTER LOOP ORDER AND BN FOR BIT SYNC 
r >2,34 

OPTION* ENTER l(ACQ), 2 (TRACK), OR 3 (BOTH) 
1 ■ > 1 


original 

OF POOR QUALITY 


ENTER #ACQ SAMPLES ( (*=100) , INC STEP (BIT), #TRI ALS 
> 100,3, 10 

ENTER 0 ( RANDOM DATA) , 1 (SQ) , 2 ( ALL1 ) , 3 ( ALL0) 

>0 

ENTER INIT TIMING, PHASE ERROR (BIT, DEG) , SEED 
>0,0,9 


PHASE ERROR (DEG) AND TIMING ERROR (KBIT) 
VS TIME (BIT) 


32.55 +- 


27.50 ♦- 


22.44 +- 


17.39 +- 


12.33 +- 


7.28 +- 


ppp 

PP! PP 

- — PP — + P 

P ! P 

-pp + PP— 

P ! P 


-p — + — + — 

P !T T TT TTTTTTT 


pppppp 

— PP — +P-PPPPP — 
P ! P 


2.23 +-PT — TTT — TTTTTTT-TTTT7TTTTTTTT-T — T +-TTTTTTP — TTTTTTTTT -TTPPPP-+ 

PPTTTT !T ! ! PPPTTT ! TTT PP 

-2.83 + + pp + p + + 


-12.94 + 


-17.99 + 

0 . 000 


59. 989 


119.979 


179. 968 


239. 957 


299.947 


PPTTTT 

! T 

; 

! PPPTTT 
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Figure E.5. A Typical Run by Using RUN BPSK. 
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APPENDIX F 

AIRS AUTONOMOUS DOPPLER COMPENSATION 
Summary 

Techniques for autonomous Doppler compensation for AIRS are 
proposed and their performance are evaluated. A first-order 
compensation technique is shown to be capable of reducing the maximum 
Doppler offset error to less than 800 H z during a powerflight if the 
Doppler compensation is applied once every second. The duration of 
maximum Doppler rate Is also reduced from 9 sec to 1 sec. The 
performance can be further Improved by reducing the time between Doppler 
compensations, Tg. An approximate relation is that the duration of 
maximum Doppler r*te is Tg, and the maximum Doppler error is the amount 
of Doppler accumulated over Tg due to the Doppler rate. 



1.0 Introduction 

Figure 1 shows a functional diagram for the AIRS Doppler 
compensation scheme. During acquisition the ADPE Doppler predictions 
are used to control the downconverter synthesizer frequency, nominally 
at 335 MH Z . After the AIRS has acquired and the carrier recovery loop 
is in lock, the AIRS takes over the Doppler compensation. The AIRS 
accomplishes this as follows. First of all, the frequency of the 
recovered carrier is generated by the AIRS tracking unit. This 
recovered carrier frequency, nominally at 35 MHz» kpown * The actual 


r 


n 

LJ 
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Doppler on the received SSA signal is the difference between the sum of 
the recovered carrier frequency and the frequency setting of the 
downconverter frequency synthesizer at the particular time instant, and 
the nominal frequency at 370 MH Z . This Doppler frequency is sampled 
periodically and the samples are used by the AIRS to extrapolate the 
Doppler frequency to the next sample point. The AIRS then uses this to 
control the downconverter synthesizer until the next sample is 
obtained. T .is function is performed by the downconverter synthesizer 
control unit and Figure 2 illustrates this technique. 

2.0 Extrapolation Techniques 

The problem of performing AIRS compensation can simply be stated as 
follows. Given a iro .u Doppler profile f(t) between 0 and t, what is 
the best way to extrapolate, i.e. , determine an estimate f^t+x) for 
f(t+x) at a later time t+r? Because of practical constraints, we only 
wish to measure the Doppler and perform extrapolation at regular time 
intervals separated by Tg secs. The problem then reduces to: given fg, 

f^ where f^ is f(kTg), how can one predict f^+j? 

2.1 First and Second Order Extrapolation 

In what follows, we consider two ways to perform extrapolation. 

The first method is exact if the Doppler profile is linear but there is 
a constant error if the profile is quadratic. The second method is 
exact if the profile is quadratic but there is a constant error term if 
this .rofile is cubic. In either case, the constant error term is 
proportional to Tn. the interval between Doppler measurement updates. 


-r 9- 




2.1.1 First-Order Extrapolation 


The first order extrapolation is given by 



z V f k -i 


to 


The associated error is 




( 2 ) 


Note that the error term (2) is the second increment of the Doppler 
profile. Note also that two measured Doppler samples are needed for the 
extrapolation. 


2.1.2 Second-Order Extrapolation 

The second-order extrapolation is given by 

r 


f k+r 3 V 3f k-l +f k-2 


and the associated error is 


(3) 


e k+l =f k+r 3f k +3f k-r f k-2 


(4) 


Notice that the error term is the third increment of the Doppler profile 
and three measurements are needed for the extrapolation. 


3.0 Performance Examples 

Tc see how the extrapolation .echniques perfor.n, we select an 
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example fron .he WDD dynamics profile us 2 d to test the cycle slipping 
performance during an uncompensated powered flight. This is shown in 
Figure 3 where a 765 H z /S Doppler rate (AIRS spec maximum) is applied at 
t-4 secs. Since the ADPE may be off by 9 secs, the current system 
Doppler error associated with the ADPE controlled downconversion is the 
difference shown. Figure 4 shows the extrapolated Doppler used by the 
AIRS downconverter when the first and second order extrapolation 

techniques are used. For simplicity, the frequency of the downconverter 

• 

is linearly interpolated between Doppler updates (2 seconds in this 
case). The Doppler seen by the AIRS carrier tracking unit is the 
difference between the actual profile and the downconverter frequency 
profile and is shown in Figure 5. Note the marked improvement. 

The doppler compensation error can be further reduced by decreasing 
Tq. Figure 6 shows that the peak error is roughly linearly proportional 
to the Doppler update period Tq. 
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Figure 1. Functional Diagram for the AIRS Autonomous Doppler Compensation Scheme 





ORIGINAL PAGt i*» 
OF POOR QUALITY 


Actual 

Doppler 















DOPPLER Ef 










)R. HZ 







JunCc 


om 


APPENDIX 6 

TECHNOLOGY ISSUES INVOLVING CCD AND SAW DEVICES FOR AIRS 


SUMMARY 

It appears that a dual 256-tap CCD correlator will be available 
next summer from Fairchild and Ford. The device is based on a mask 
developed by S. C. Munroe of the MIT Lincoln Laboratory. The device 
appears to be directly compatible with the AIRS acquisition 
requirements. Because of the size of the bulky substrate required, SAW 
correlators are deemed to be unsuitable for the AIRS application. 

CCD Devices 

The attached paper describes work performed by MIT Lincoln 
Laboratory on the development of a 256- tap CCD correlator prior to 
1982. Subsequently, they have developed a mask for an improved version 
of the TC1235A device to be transferred to Fairchild [Contact: David Wen 
(415) 858-6166] and to Ford Aero [Contact: John Roschen (714) 720-6151] 

for commercial (military?) production. The device should be available 
by next summer. 

o 

The new chip has an area of 50,000 min and is based on 4 pm 
technology. A lot of the external functions shown in Figure 3 of the 
paper will be incorporate^ into the new chip. The new chip will be a 
dual correlator accepting both the I and Q signals and can correlate 256 
samples (512 ISO) simultaneously. The inputs to the new chip are the 
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reference PN code, the associated clock and a latch signal. The 
reference PN code is clocked sequentially into the chip and the latch 
signal controls the timing which determines when the PN chips will be 
loaded in parallel to form the fixed binary weights for the 
correlation. While the correlator chip is performing the correlation, 
it can accept the next set of sequential PN code reference input. Since 
the latching of code weights takes place within 400 ns, there is no 
problem in updating different sections of the PN codes in sequence. 

This simplifies the hardware requirements since we do not need another 
device for setting up the next set of weights, fhe output of the chip 
interfaces directly with an A/D converter. 

For the AIRS application, it is anticipated that up to 1,000 PN 
chips will need to be correlated. This means that 8 of these devices 
must be used in parallel, each matched to an appropriate portion (128 
chips at 2 taps/chip) of the total section of 1,000 chip code. Since 
the A/D operates at a relatively slow speed of 3 MHz, accumulating the 8 
outputs does not appear to be a major hardware problem. 

SAW Devices 

A programmable SAW tapped delay line consists of a single 
piezoelectric substrate with an interdigital transducer at one end which 
launches a SAW toward an array of interdigital tap transducers. Each 
tap transducer is wire bonded to a substrate containing integrated 
circuit switching chips which combine the taps with their respective 
polarities chosen to match a selected code sequence. Typically the SAW 
travels at a velocity of 3x10^ cm/sec which corresponds to a delay of 
3.3 us per cm of substrate. Correlating a 512-chip sequence for the 3 
Mcps T0RSS chip rate (171 ps) requires a 52 cm length of substrate. 
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THE USE OF CCD CORRELATORS JN THE SEER COMM 
SPREAD SPECTRUM KECEt VEKSt (U) 

, S. C. Munroe __ 

( ( n? - #/z - sr-,'^-3 

M.I.T. Lincoln Laborator, 

P.O. Box 73 

Lexington, Massachusetts 02173 


CU) ABSTRACT 

(U) The SEEK COMM radio uccs • pair of 
Charge-Coupled Device (CCD) correlators to detect 
the dehopped , pseudo-noise nodulated signals in 
each of three receiver channels. Developed Jointly 
by H.I.T. Lincoln Laboratory and RCA, the TC1235A 
analog-binary correlator had specification goals 
which were strongly driven by SEEK COMM systen 
requirements. With minor exceptions our experience 
with the device performance. Interface 
requirements, and the total analog-binary 
correlator board have been very favorable. In 
terms of performance and practicality the TC1235A 
Is a major advance over earlier correlators. The 
correlator boards have been designed to permit 
plug-in replacement of the CCD correlators with o 
single, one-time gain trim. To date, however, not 
one TC1235A has been replaced due to failure. The 
minor deficiencies in the device have been 
sufficiently characterized, and are well enough ‘ 
under' -ood, to be either minimized or eliminated in 
the next generation chip. 

TABLE 1 


(U) INTRODUCTION 

(U) CCD correlators were chosen and developed for 
uee in the SEEK COMM receivers primarily because e j 
the advantage analog-binary correlators have over 
all digital correlators in a jamming environment. 
Rather than dwell on device design and developsee; 
which have been covered previously, [1,2,3], the 
emphasis here will be on CCD correlator system 
requirements and device performance, Interface 
requirements, board level performance, and 
experience gained. 
jr • 

(U) SYSTEM REQUIREMENTS AND DEVICE PERFORMANCE 
(U) The design of the SEEK COMM receivers placed • 
.mber of requirements on the correlators. Many t i 
these specifications were translated directly lnte 
design goals for the RCA TC123SA CCD correlator. 
Table 1 summarizes both the system requirements aid 
the actual device performance. The following 
comments apply to that table. 


i : 
\ 1 

i 
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(U) CCD CORRELATOR SPECIFICATIONS 
System Requirements 

Device Length: 512, 256 points; programmable 
CCD Clock Rate (> 2X chipping late): d MIU 
Program (Reference) Register Lcffd (Lite: ? MHz 
Reference Latch Update Time: «3ai>s (symbol 1-ngth) 
Dynamic Range*: > 40 d3 (full length mode) 

Signal-to-Noise Ratio: ' 35 dS 

Distortion: “Low” 

Code Incensltlvlty (to code balance, distribution) 

Device Interchangeability 
Low Power Dissipation: { 2U 


TC1235A Performance 
512, 256, 128 

> 10 MHz 

> 10 MHz 
■ 5 pa 

> 50 dB 

> 35 dB 

> 1Z 

No code-dependent bias, some 
devices sensitive to distribution 
of I’s and 0's. 

Yea - with one-time 
gain trim. 

900 mW 


Reasonable off-chip support circuitry 
‘Dynamic range defined as 20 log 


Yes 


Signal irput at max, output 


UNCLASSIFIED 


10 Signal input at min. detectable output 


IThls work was sponsored by the Department of the Air Force. The views and conclusion!, contained in this 
document are those of the contractor and should not be Interpreted as necessarily -epresentlng the official 
policies, either expressed or Implied, of the United States Government. 

DISTRIBUTION LIMITED TO US GOVERNMENT AGENCIES . Approved for presentation in closed session and 
publication. Acquests lor additiooal copies of this document should be obtained from HQ ESD/TOFL (Llncols 
Laboratory), Hanscom AFB, Bedford, Massachusetts. 
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IJI Tfc, nrrd In alr lalit lnifrmoHul.it Jon between 
,,,**1 and jammer >e the distortion 
•a«< Ideal Ion. 

j-j t>i« code-lnscnsltlvlty issue haa two aspects. 

first Is the sensitivity of the correlator to 
Iwawlanrrs In the numbers of I's and O's In the 
«fe“code 1*1. In the TCI235A this problem 
bvtn designed out*''. The second aspect Is 
M»elti«ity to the distribution of I's and O's In 
|Wm refer enre code. The aource of this subtle and 
saeapected problem, which reached unacceptable 
I avals In only a ‘ew correlators while In the 
fali~i*ngth mode, was traced to differences 
latMin the first and second halves of the 
Saw tees. These nonunlf oral ties can occur because 
I ka chip Is folded at midpoint and js r » 
reslarratlon errors parallel tc the CCD channel 
adact the first and second halves differentially. 
|y eliminating several extremely asymmetrical 
rmdas and Including the r k ,rrelator In an auto-zero 
lamp the problem has been nypassed in the SEEK 
°°*i radio. With the understanding of the 
machantsa gained from the 7C1235A It should be 
paarlble to iiuuce code seu>itlvlty to negligible 
le> ala in the next generation CCD correlators. 

CC) figure 1 Is an oscillogram showing correlator 
performance while single sampling a 3V p _ p , 8 MHz, 
255-blt, cyclic M-sequence As expected, a single 
terra latlon spike occurs every 32 ps. While there 
la aoae nonuniformity suf’-l *>r>sed upon the 
roaatant clock feedthrough In ; e baseline, lta 
value la consistent with tL. . ice being 
•1 ••etched to the cyclic code oy 1 bit (256 va 255 
bite). In terms of peak-to-sidelobe ratio, 

I — refore, the TC123jA approaches the theoretical 
value for cyclic M- equences. 


(U) Figure 2 la a photomicrograph of the TC1235A 
mounted and bonded In a 64-pln ceramic DIP. Aa 
can be aeen not all pin* are used, and some of 
rhe bonded pins are redundant. The uhlp Is 
folded at .he midpoint, thereby producing mirror 
Image symmetry about a cental horizontal a^ls. 
Although the chip la large It was fabricated with 
8 pm design rules. If fabricated with today's 4 
pa rules the area would more than halve, making It 
comparable In area to a 64 Kbit DRAM. Design 
changes contemplated in aoae circuits for the next 
generation will shrink power, and pin count 

still further while enhancing performance. 
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TC1235A CCD CORRELATOR 
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*• L (U) TC1235A CCD correlator. 
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Fig. 2. (U) TC1235A CCD correlator chip. 


(U) SEEK COMM A/D SLICE 

(U) The analog-binary correlator L-*ru In the 
SEEK COMM radio la often referred to as the A/D 
slice because all signals routed to the board are 
analog, while all those leaving sire digital. 

Figure 3 Is a block diagram of the slice. Due to 
the I and Q nature of the signals two correlators 
are required per receiver channel. For diagnostic 
reasons analog switches have been placed toe. ore the 
correlators so that u tent Input can be routed to 
either or both device*. Thla feature facilitates 
rapid Isolation of fi.t.s. 

(U) The CCD off-chlp s .pport circuit-/ consists of 
the following. A single 20V p _ p CCD shift riglste- 
transfer clock and 8V p _ p input sampling (strobe) 
pulse are generated In are circuit. The load logic 
circuitry controls loading of the reference code 
and updating of the latches. The on-chlp program 
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Fig. 3. (U) TC1235A analog-binary correlator board block diagram. 
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register clock driver is used here, and both the 
input reference code and clock input to the on-chip 
driver are TTL compatible- Length control switches 
enable the correlator to be either 512 or 256 
stages long for sync or data Bodes, respectively. 

A very Important support circuit consists of two 
feedback loops per correlator which coopletely 
stabilize the devices against threshold voltage and 
teaperaturc variations. These coapensation loops 
enable plug in replacement of the correlators with 
only a one-time gain trim performed during board 
calibration. If an emitter resistor can be 
Justified to stabilize a single bipolar transistor 
against wide variations in 0 then surely a dual op 
aap and a few passive cooponents can be Justified 
to stabilize an 1C with nearly 10.000 transistors 
and many orders of aagnltude more signal processing 
capability. Incorporation of this automatic 
compensation circuitry, which consumes little power 
or area, was a major step toward the goal of a 
"rweahless" A/D slice. When the CCD correlators 
are In production status, so that the 'abrlcatlon 
conditions are more uniform, the 110Z maximum gain 
variations seen now should be reduced considerably. 
Once the maximum correlator gain variations are 
less than about ±5X then even the one-time gain 
trim can be eliminated, therby resulting in a 
completely adjustment-free board. 

(U) The post correlator processing circuitry begins 
with the transformers, which provide scaling and 
level shifting of the differential correlator 
outputs with high common mode rejection. Cain 
adjustment is effected by trimming potentiometers 
on the secondaries of the transformers. Fast, 
pirclsloo rectifiers convert the bipolar 


differential Inputs Into easily sucoed unipolar 
current source outputs. The combined currents arc 
fed to a fast op aap for voltage conversion and 
scaling, with the result routed to an C-bit flash 
A/P converter after modification by a 
not linearity. The purpose of the nonlinearity is 
to expand the low signal end (by 2X) and compress 
the high end (by 2X), thereby achieving 9-blt 
resolution. Finally, an auto-aero loop forces the 
A/D output to be a digital zero for the zero 
analog signal condition. Our experience has bees 
that an auto-zero loop would be necessary even If 
the correlators were perfect because of the 
temperature-dependent input offsets of both the 
fast op amp and the A/D converter. However, the 
auto-zero loop could be simpler, and the fuoctlcm 
implemented less frequently, if the correlators 
were outside the loop. By reducing the 
sensitivity of the next generation of CCD 
correlators to asymmetrical codes it should be 
possible to do this. 

(U) Figure A is, a photograph of the SEEK CKSI 
analog-binary correlator board. In the upper lef: 
corner are the l.t), and test inputs, together wiil 
two analog switch IC's. The two CCD correlators 
are In the upper right quadrant. Nearby are the 
feedback compensation loops, length, and program 
register control circuits. In the lower center w 
the transfer and strobe clock circuitry. The A JZ 
converter and support circuitry is in the lower, 
left hand corner, while the digital timing cktlz 
Is In the lower right. For convenience the 
rretltlers, fast op aap, nonlinearity, and 
auto-zero circuitry have been placed on a dauff-v 
board mounted behind the correlators. 
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Pig- 4. (V) TC1235A analog-binary 
correlator board. 


Fig. 5. (U) TC1235A analor-binary board. 


CC) Th« A/D * 1 lee boards are completely 
latarchangeablc. * n <* dissipate about 12W typical, 
••tk i be eicept loo ot power and dynamic ran;c the 
•Hclllutlogr for the board are almost Identical 
Iboae Hated la Table 1 for the TC1235A Itself. 
*** “west case dynamic range measured for some 
r ^* !• about 35 41, down from about 50 41 
****** fot the device alone tested with cyclic 
•!<>»((•. The typical dynamic range la about 
41. with highs near 50 dB. One factor 
•■mirtbutlng to the dynamic raoge reduction on the 
*4 la the more severe test conditions of 
•avrlodle codes which are not M— sequences- This 
a tuatlon Illuminates defects such as sensitivity 
I* **p°*trlcal codes. Another factor Is 
c ^* r In the post correlator processing 

^ttultry. l n any case, elimination of several 
J*"*ly »|rsattrlt«l codea rendered hoard-level 
101 1 ™**" *° r * than adequate for the radio. 

I, board and CCD correlator reliability have 

Figure 5 la an os-Ul-'g .an 
U , ** board performance after conversion of the 
jg * ou tput to analog fora by a D/A converter. 
f #r *^" P ***^ u lth Fig. I. which shows performance 
c °rtelator alone, it la clear chat other 
4 lff * * c *l e change the only significant 
•““wolf'** **' C w> »»eforms la the Increased 

Jhj # . or "fty In the baseline of the board output. 
m>g|| la due primarily to the Intentional 

!„, ,******* on the board which amplifies low 
*!gaals by a factor of two. 


(U) SUMMARY 

(U) The TC1235A Is an analog-binary CCD correlator 
developed Jointly by RCA and MIT Lincoln 
Laborat .y for use In the SEEK COMM receivers. It 
fa a third generation device which exhibits a 
quantum leap In performance over earlier 
correlators- When used with proper support 
Circuitry the TC1235A displays reliability and 
capabl lty more than adequate for use In spread 
spectrum receivers. Its few flaws, which are 
relatively minor, have been characterited and are 
understood well enough to he significantly reduced 
or eliminated In the next generation. 
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OPTIMAL GAIN CONTROL FOR ADC 


H.l Introduction 

Gain control is required to match the dynamic range of the input 
signal and the analog-to~digital converter (ADC) to minimize distortions 
introduced by the A/D system. Here, a techniuqe for generating the 
required feedback to adjust this gain is described. The technique also 
provides a May to estimate the received signal level. 

H.2 System Model 

Figure H.la shows a system model for analysis. The incoming signal 
is modeled by 


s(t) = d(t) + n(t) 

where d(t) is a t/S NRZ waveform with bit time T and n(t) is white with 
density Nq/2. The purpose of the low pass filter is for antialiasing 
and is assumed to be ideal with bandwidth B Hz. It is also assumed that 
BT » 1 so that the signal is undistorted at the filter output. The 
gain in front of the ADC is adjusted so as to minimize the probability 
of error at the output of the digital integrate and dump (summer). 

Figure H.lb shows the model used for the ADC. The output of the 
gain controlled amplifier is sampled and then clipped to the range of 
(-V/2.V/2) where V is the dynamic range of the ADC. The clipped signals 
is then quantized to 2 n levels where n is the number of bits used by the 
ADC. The normalized converter characteristics is 
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for kA< y < (k-1) a 


where a= 2~ n is the normalized quantization step size of the ADC. 

Assuming the IAD clock for the summer is coherent with the data 
transition, i.e., perfect bit sync available, the input for each I&D 
period is a sequence of signal samples of the same sign. Without loss 
of generality we assume that a +JS signal is sent... Then, the input to 
the clipper is a Gaussian sample with mean = G«£ and variance = 
G^NqB. The probability of the sample x k at each of the levels of the 
ADC output can be computed easily. In particular, the probability of 
assuming the lowest (1) P^, the highest (2 n ) level P H , and the sum of 
the two, P Qut are of interest as we shall see shortly. They are given 
by 


P H = \ erfc C(Y - a)/y-1]^p 
P L = |-erfc[(^ -A)/tfl]/p 


P out * P H + P L 


where y 3 G/S/V, p = S/2o and erfc(x) = /2/» / e dx. After computing 

x 

the probability of each ADC output levels, the probability distribution 
of the summer output can be easily determined, for example see [1], The 
probability of making an error, i.e. P^ can then be obtained. 
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H.3 Dependence of Pp on A/D System Parameters 

The resultant probability of error is a function of the ADC input 
signal-to-noise ratio, the integration length N, the number of bits n 
used, and the setting of the gain G. Figures H.2-H.5 show the 
performance of the digital I&D system compared to an ideal analog I&D 
system. The degradation, in terms of signal-to-noise ratio loss at a 
fixed E b /NQ, is plotted vs the normalized gain y/ p as a function of 
different combinations of E^/Nq, n and N. 

In a practical system implementation, one must somehow use 
information from the ADC outputs to adjust the g$jn G. The last set of 
curves in Figures H.2-H.5 unfortunately are related to the system 
parameters in a complicated way so that they can be of limited use for 
error signal generation. Intuitively, it is reasonable to assume that 
the AOC outputs, regardless of the A/D parameters, must be 
“statistically preserving" of its samples, in order to optimize the 
system performance. Figure H.6 shows the performance degradation of the 
digital I&D system in a form motivated by this intuitive reason! n The 

degradation is plotted against the out of range probability P 0ut for a 
variety of input signal conditions and integration length with n fixed 
at 6 bits. From this figure, it is obvious that if P ou j. is restricted 
to 2’ 5 = 1/32 or less, the degradation is less than 0.02 dB. The 
approach suggested by Figure H.6 to adjust for optimal gain is to keep 
track of the occurrences of the highest and lowest samples and adjust 
the gain G so that P out - 1/32. 

H.4 Analysis of Proposed AGC/Signal Strength Indicator 

Figure H.7 shows the out of range probability as a function of the 

Q 

normalized gain y//p = ^ ^ owith varying system parameters E b /N 0 , N 
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Figure H.3. Digital Integrate-and-Dump ADC Loading Sensitivity. 
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Figure H.4. A Digital Integrate-and-Dump ADC Loading Sensitivity 
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and n. It can be observed that P out can serve as an error feedback to 
adjust the gain. For a fixed P out no system parameters, it is shown in 
Figure I!. 8 that the mean output of the ADC m is a good estimate of 
signal y = G/Sf/V. If one measures the mean output of the A/D, then one 
can estimate the input signal amplitude /S via 


. m 

Since V is a known hardware parameter and G is the steady state gain of 
the AGC system. Figure H.9 shows the algorithms^required to implement 
this A'lC/signal strength indicator system. 


Reference 

[1] C. M. Chie, "Performance Analysis of Digital Intpgrate-and-Dump 
Filters, 11 IEEE Trans. Comm. Tech. Aug. 1982, pp. 1979-1983. 



I 

















m 


APPENOIX I 

BASEBAND VS IF I-Q SAMPLING 


Figure 1.1 shows two approaches to obtain inphase and quadrature 
(I4Q) samples of an IF signal s(t) = fi . P o(t) sin *^t. The first 
approach is straightforward, it uses two A/D converters to sample the 
demodulated baseband outputs after the mixing operation. The second 
approach samples at 4 x the IF frequency and yields 2 pairs of I-Q 
samples every IF cycle as shown in the example in Figure 1.2. Notice 
that the I and Q samples have to be inverted alternatively because of 
the sign change of the IF carrier every half a cycle. 

Obviously, the IF approach saves hardware, namely, 2 mixers and 1 
A/D converters. However, there are disadvantages. First of all, the 
sampling speed is currently a technology constraint. Given a fixed 
sampling speed the baseband technique offers two times more samples, 
hence two times more resolution in time. Secondly, in the IF sampling 
scheme, the samples must be taken precisely at multiples of */2 of the 
IF carrier. However, for data detection purpose, the system performance 
can be enhanced if the samples are taken coherent with the data clock as 
done in the baseband scheme. This advantage becomes more apparent as 
the number of samples /symbol decreases. 

Let us consider an example. A state-of-the-art commercial flash 
A/D converter has a conversion speed of about 60 MHz. For 6 Msps 
operation there are 10 samples/symbol. With the I-Q sampling, there are 
only 5 samples per channel. Even if the bit sync can operate perfectly, 
the inherent time quantization error has a peak-to-peak value of 1/5 = 
20% since the sampling instants are forced to be coherent with the IF 
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Figure I . Z IF I-Q Sampling 
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carrier instead of the bit sync. This also affects the integrate-and- 
dump detector performance. On the other hand, since the sampling clocks 
in the baseband technique are controlled and adjusted* by the bit syncs, 
the quantization error can be made to be very small. 

Of course, at low data rate where there are many samples/symbol, 
the IF technique becomes attractive. 


*To the resolution of the synthesizer providing the clock, which can be 
made very fine. 
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MANCHESTER SYMBOL AMBIGUITY RESOLUTION 

When the baseband data is Manchester-coded, there is a need to 
resolve the ambiguity associated with the symbol alignment. Figure 0.1 
illustrates this problem. Depending on the bit alignment, the symbol 
stream can be interpreted either as all ones or zeros in this example. 
If the data symbols are correctly aligned, there is always a transition 
in the middle of the symbol. If the data symbol is purely random, then 
the occurrence of a transition in the middle of the misalignment symbol 
is only 50%. 

By keeping track of the number of transitions using the two 
different symbol alignments, one can resolve this ambiguity 
statistically. Let us assume that we keep track of M * 2 N symbols and 
see which symbol alignment yields the most mid-symbol transitions. Let 
p c be the probability of error in making a correct decision on the half 
symbol so that 


= j erfc 


A 


2 E s /N 0 


A transition is declared if the first half of the detected symbol is 
different from the second half. Given that the correct symbol alignment 
is used, we can easily determine the following probabilities 


Pl ■ U-P c ) 2 ♦ p l 

<*1 = ‘-pi 
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Figure J.l . Manchester Symbol Ambiguity 
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of detecting the transition (p^) and not detecting the transition (q^) 
at the middle of the symbol. The probability of detecting a transition 
at the wrong symbol alignment is 


p 2 ~ P T *Pi + 


0 


where pj is the probability of having a transition. Notice that pj is 
related to the data transition density of the symbol pattern, p t via pj 
= l-p t . We also use q 2 = l-p 2 to denote the probability of not 
detecting a transition. , . 

After observing for M = 2 N symbols, the probability distribution of 
the number of observed transitions k, at the correct alignment is 


And for the incorrect alignment. 


1 


P 2<*> - (>2^‘ 

The probability of picking the incorrect alignment based on the number 
of observed transition in both cases is the error of the ambiguity 
resolution process and is given by 


p£ = Prob(k<t) 

M * 

- I l PiOOP 2 U) 
*=l k-l 1 c 
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Figure J.2 shows the probability of incorrectly aligning the Manchester 
symbol. This is computed for the worst case data transition density p t 
= 0.125 per the TDRSS requirement. As an example, for a worst case 
E b /NQ = 5 dB, P c < C.05 so that the ambiguity can be resolved with an 
error probability less than 10"^ if 1,000 bits are used. 






Figure J.2. Manchester Symbol Ambiguity Resolution Performance 
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LOCK INDICATION 


Since the AIRS operation requires switching to different 
configurations at various stages of the receiver operation, means must 
be provided to monitor and verify whether what stage is the receiver 
operating at. This is the task of the lock indicators. In AIRS, lock 
indication is derived from comparing two signals, (selected to be 
insensitive to other system conditions), of which one signal will be 
high relative to the other when the system under .eonsiderai ton is in 
lock, and vice versa. This avoids, for example, the need to establish a 
fixed absolute voltage reference, which may be difficult to maintain 
during acquisition. 

As an example, the lock indicator for the carrier loop is obtained 
by comparing the in-phase output (proportional to |cos <j>| ) and the 
quadrature output (proportional to (sin <j>|) of the A/D subsystem. Since 
the gain of the receiver is identical for both signals, the relative 
magnitude of the two signals is independent of the receiver power level 
or whether the AIRS is in the acquisition or the tracking mode. 

Figure K.l shows the signals to be compared for lock indication. 

The signals used for the FLL is proportional to |s1n A<uTp| and |cos 
AoTpl where Tp is the update speed. Before the FLL acquires, both 
signals are of comparable value. After the AIRS has acquired, j si n 
AuiTp| is very small compared to |cos Aalp|. 

For the bit sync, the signi. a are comparable if the bit sync is not 
locked. Once locked, the full-bit ACM output is much larger than the 
• mid-bit output. 
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K.l. Lock Indicator Signals for AIRS Subsystems. 
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The lock Indicator signal for the AGC is a little different. The 
signal measures the probability that a A/D sample is out of range, if 
the AGC is set incorrectly, then this probabil. j will be different from 
the design value (e.g. 1/32). If the gain is too high then the 
probability will be greater than 1/32, and vice versa. When the system 
is in lock, then the measured value will be within a prescribed range 
around the design probability. 

Analysis 

The performance of these lock indicators can best be determined via 
analysis than via a Monte Carlo simulation since .the latter is very time 
consuming in the range of interest. Tne analysis is not difficult and 
is straightforward. Figure K.2 shows the analytical approach applicable 
to AIRS. These procedures can be followed during the final design phase 
of the AIRS. Because M is usually large, the law of large number 
applies and \ | A | and \ J B | can be replaced by Gaussian random 
variables. In that case, the analysis can be simplified further. As a 
matter of fact, standard results can be used and are found in radar 
literature under the topic of noncoherent integration. 




- 309 - 









i r 


t 

K 

I 

I 


oCindt 


om 

APPENDIX L 
MEAN SLIP TINE 

The mean time to first slip t for a CW PLL is given by [1] the 
approximation 


2 .cos * ss +26t ss 


T = 


4B l cos cosh *8 


where a= linearized loop signal-to-noise ratio, 4 ^ $ is the steady state 
phase error, and 6= a sin If we approximate the BPSK and QPSK 

loop S-curves by sin 2 4 and sin 44 respectively, then with a simple 
substitution, the mean time to first slip can be shown to be 

^cos 

*e ' 

4B l (cos N4» s$ ) cosh *^N 


t - 


where 3 = ^sin N<» ss and N = 2 (BPSK) or 4 (QPSK). For simplicity we 

can assume that is zero since it is negligible with a third-order 

loop for the TDRSS Doppler profile. It is instructive to express a in 

2 

terms of the steady state rms phase jitter o^. Then we have an 
approximate formula for the mean time to first slip 


t = 


, 2/N 2 o 2 

— e ♦ 
4b l 


where N is the number of phases in the modulation scheme and a. is in 

♦ 

radians. 

As an example, let Bj_ » 50 Hz. Then to achieve a 90 minute mean 
time to first slip, we require 
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N < 2/lr»(4B L x/ w) 


or o < 22.7°/N. Hence BPSK requires a. < 11.3° and QPSK requires 

? ♦ 

o < 5.7° to meet the slip time requirement. 
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